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The principal problem of theoretical chemistry is the elucidation of the 
causes of molecular stability. The first successful explanation of this 
stability was the first-order perturbational atomic orbital treatment of 
Heitler and London. 

This was followed by the development of a quite different method by 
Hund, Mulliken and J. Lennard-Jones, the method of molecular orbitals. 
It should perhaps be emphasized at this point that this is also only a first 
approximation. 

The principal explanation of chemical valency forces, according to the 
atomic orbital treatment of Heitler and London, arises from the indis- 
tinguishability of similar particles in quantal space. The eigenfunctions 
describing the molecular system are composed of atomic orbitals which . 
retain their original spherical symmetry. 

On the other hand, molecular orbital theory ascribes chemical valency 
forces to an extreme distortion of the individual one-electron eigenfunctions, 
the product of which is the molecular eigenfunction. 

Depending on our choice of approximate molecular wave functions we 
get two different interpretations of chemical valence theory. 

The atomic orbital method is more in accord with classical valence 
theory, and each molecule does not pose a special problem, as it does in 
molecular orbital theory. The method also has the advantage that it 


appears to be precisely what it is, namely, a mere first approximation, 
The chief disadvantage consists in its complexity. 
Using molecular orbital theory each molecule is a more or less special 


problem, not only quantitatively, but qualitatively. The general pro- 
cedure for such calculations requires that we assign various orbitals to the 
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electrons in such a way that each orbital belongs to some irreducible 
representation of the point symmetry group of the molecule. The net 
result is a high degree of orthogonality between the various orbitals and 
a resultant simplification calcu- 
lationally. At the same time, 
such orthogonality places an 
unnecessary restriction on the 
orbitals. Theoretically such or- 
thogonality is only required 
between eigenfunctions repre- 
senting different eigenvalues of 
the entire molecular system. 
While both of the preced- 
ing methods have yielded qual- 
itatively information of great 
value, they have not been very 
successful in calculating the 
things of primary interest to 
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dipole moments. It seems imperative, therefore, that we pass from 
first to second approximations, using a self-consistent field method. 

As a first approximation to the solution of the Fock-Hartree equation 
we shall use an eigenfunction capable of assuming molecular or atomic 
orbital forms as special cases. 

We shall apply this method to the hydrogen molecule and then generalize 
the results to more complicated cases. 

The Self-Consistent Field Applied to the Hydrogen Molecule.-Using the 
method of Fock! we find the following set of equations for the hydrogen 
molecule : 


L TSS op cast yf (vs(2))* 7 
<r 4 / gi(2) 1 ¢(2)dr14 ¢(1) + f a dra ei(l) + 
* gi(2)0)(2 
if gilZ)e an ¢(1) + | / e(2)He(2)dr et) |= e¢i(1); (1) 
2 * 


e, is the energy of the 7th electron in the coulombic and exchange field of 
the other electron. 

The correct antisymmetric eigenfunction for the ground state of the 
hydrogen molecule has the form 


I 
Y= { oi(1)¢)(2) + ¢:(2)¢,(1)} f a(1) B(2) — a(2)B8(1)}. (2) 
V2 + 2S? 
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As a first approximation we use the eigenfunction 


V { (yall) + Awe(1))(Aga(2) + on(2)) + (¢a(2) + 


1 
V2 425 


Ago(2))(Aga(l) + go(1))} (3) 
where ¢, and ¢ are hydrogenic 
ls orbitals. This eigenfunction 
was first used by Coulson and 
Fischer.’ 

We can then take advantage 
of a transformation to elliptic co- 
ordinates 
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where 
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mation. 

In the simple transformation, apart from the sign, the coefficients of 
uw and » are equal. By destroying this equality we replace the atomic 
orbitals by elliptical eigenfunctions 


e7 7? 1» e - (ap + Br) ero e ~ ap + Bry (5) 


and the eigenfunction of equation (3) becomes 
l 


V ie . le ~ ar ai + a) [(e~ + de®) (Xe ~ Bre + e”*) + 
NV 2 + 28? 


(e~ Met MPm)(ne-O" + €")]} (6) 
where N is the normalizing factor. 

Applying the self-consistent field method we find that a = 0.835, B = 
0.775, X = 0.137. Such an eigenfunction gives a binding energy of 4.20 
e.v. or 89.0 per cent of that observed experimentally. A comparison of 
the resultant accuracy is given in table 1. 

The first approximation to the solution of the Fock-Hartree equation 
has been called an antisymmetric molecular orbital. For reasons which 
will appear la‘er we shall refer to this first approximation, the eigenfunction 
of equation (6) and subsequent generalizations as semi-localized orbitals. 
Briefly the reason for this is that it is not yet clear that in complicated 
molecules even the first approximation related to the eigenfunction of 
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equation (4) will involve the atomic orbitals of all the atoms in the molecule. 
Secondly, the formalism of the method of semi-localized orbitals is pre- 
cisely that of atomic orbitals and the final molecular eigenfunction has the 
symmetry properties of molecular orbitals only for the extreme case of 
exact coincidence with molec- 
ular orbitals. 

The charge distributions of 
the one-electron eigenfunctions 
of atomic orbital, molecular 
orbital and semi-localized orbital 
methods are given in figures 1, 
2 and 3. 

Extension of the Method to 
Complex Molecules.—The_ ex- 
tension of the preceding method 
to more complex molecules is 
made through the use of the 
formalism previously developed 
for the atomic orbital method. 
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Ga > e Me™ 
~ yu, t+ bv 


¥h = ¢€ e (8) 


or even more simply by the analogous combination of atomic orbitals. 


TABLE 1 


HypDROGEN MOLECULE 


BINDING ENERGY . TOTAL ENERGY 
METHOD g.V KCAI E.V KCAL 


A.O. Variation" 3.76 79 66 30.965 96.1 
S.C.F.M.O.° 3.608 76.27 30.808 96.79 
SCF Semi-localized 4.202 8902 31.407 98 38 
Orbital 

* Wang, S., Phys. Rev., 31, 579 (1928) 

® Coulson, C. A., Proc. Camb. Phil. Soc., 34, 204 (1938). 


As an example of the simplest sort of eigenfunction, we will report 
briefly on the application of the method to diatomic hydrogen fluoride. 
Here we have utilized the eigenfunctions 


(Ga + Aig); (Ax¢a + ¢p)- (9) 


Using a complete Hamiltonian, neglecting multiple exchange integrals, 
and omitting all canonical structures except the one corresponding to 
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classical valence theory, binding energy and dipole moment calculations 
were made. Since the molecular orbital and atomic orbital calculations 
are special cases of semi-localized orbitals, these calculations were made to 


serve as comparisons. The results are given in table 2 and the calculation 
will be reported in more detail in a subsequent publication. 


TABLE 2 
DIATOMIC HyDROGEN FLUORIDE 


BINDING ENERGY DIPOLE MOMENT 
METHOD IN ALU IN Db. 


Atomic Orbital —0. 9688 0.86 
Molecular Orbital 0.1352 0.67 
Semi-localized 0. 1648 1.56 
Observed 0.2440 1.91 


In this case the semi-localized orbitals are much closer to the molecular 
than the atomic orbital treatments. The contrast in dipole moments 
between all three methods is very noteworthy. ‘The final form for the 
three treatments is given below. 

Atomic Orbital: 

6° 
%. = = ¢ ** Fr, COS Gy 6 = 2.60% = 
Vr 

Molecular Orbital: 

Ga + 0.571 >. 

Semi-localized Orbitals: 


(ga + 0.14¢,); (1.04¢, + ¢). 


In this calculation the inner electrons are assumed to be unchanged. 
This assumption is borne out by the small amount of distortion which the 
fluorine atomic orbital undergoes. This distortion is assisted by a large 
exchange force. In the case of the other non-bonding electrons the dis- 
tortion is opposed by slightly weaker exchange forces. If the valence 
electron orbital of fluorine is only slightly distorted, it seems likely that the 
non-valence electron orbitals are practically unchanged. 

Calculations with the semi-localized orbitals have now been made for 
two extreme cases of chemical binding. In both cases considerable im- 
provement in numerical results has been obtained. Several general con- 
clusions seem to be in order. Indiscriminate use of either molecular or 
atomic orbital methods may lead to nonsensical results, for example, the 
atomic orbital calculation of the binding energy of hydrogen fluoride. 
In any domain the semi-localized orbital method is necessarily superior. 
In the two cases considered the molecular eigenfunction belongs to an 
irreducible representation of the molecular point symmetry group. The 
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individual electronic eigenfunctions do not by themselves belong to any 
such irreducible representation. Hence the molecular orbital formalism 
is not correct for either molecule considered. It seems likely that this 
formalism may hold only for the trivial case of a completely ionic bond. 
The basis for such a conclusion is that for the highly ionic bond of hydrogen 
fluoride the identity of semi-localized and molecular orbitals is far from 
being complete. 

In any case the exploitation of the method of semi-localized orbitals 
appears to be a profitable undertaking. 

* This paper was prepared in connection with a conference on ‘‘Quantum-Mechanical 
Methods in Valence Theory” held at Shelter Island, September 7-10, 1951, under the 
auspices of the NATIONAL ACADEMY OF SCIENCES (see ‘‘Summary of the Conference” in 
a subsequent issue of these PROCEEDINGS). 

! Fock, V., Z. Physik, 61, 126 (1930). 

2 Coulson, C. A., and Fischer, I., Phil. Mag., 40, 393 (1949). 


SUPERPOSITION OF CONFIGURATIONS: THE HELIUM ATOM* 


t 


By G. RUSSELL TAYLOR' AND ROBERT G. PARR 


DEPARTMENT OF CHEMISTRY, CARNEGIE INSTITUTE OF TECHNOLOGY, PITTSBURGH, 
PENNSYLVANIA 


Communicated by R. S. Mulliken, January 7, 1952 


Introduction.A very nearly exact approximation to the ground state 
wave function for the helium atom was obtained by Hylleraas in 1930.' 
His approximation is so accurate that calculation of the ground state 
wave function and energy for helium no longer remains an interesting 
problem per se, but the extension of his method beyond two-electron 
problems is quite hopeless.2, The work reported in the present paper 
therefore was undertaken as a step in the development of a systematic, 
practical procedure for obtaining an accurate wave function for any atom 
or molecule. 

The quantitative shortcomings of the methods that have been most 
commonly applied to many-electron systems heretofore, the Hartree-Fock 
method for atoms and molecules and the valence bond and molecular 
orbital methods for molecules, are becoming more apparent with each 
new application. There has resulted an increased use of wave functions 
involving configuration interaction.* 4 Wave functions of any desired 
accuracy in principle can be obtained by the superposition of a sufficiently 
large number of configurations,® but if the method is to be of practical 
value the number of configurations required should be small; i.e., the 
configurations should be chosen with discretion. Present knowledge re- 
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garding configuration interaction effects is very meager, however, so that 
inspiration now serves in place of discretion. Prospects of progress have 
been termed melancholy,® but melancholy or not there appears to be no 
other road. At the beginning lies the simplest non-trivial example, the 
helium atom. 

In the problem of the helium atom the number of simplifying assump- 
tions is a minimum, the required integrals can be evaluated to any desired 
accuracy without prohibitive effort, and accurate experimental data are 
available for evaluation of results. Furthermore the relative simplicity 
of the system permits the variational determination of the best effective 
charges to be used in one-electron functions, which is not feasible for most 
more complex systems. Simplicity, of course, has as its concomitant 
disadvantage the loss of generality. Nevertheless, conclusions reached for 
the helium atom can be used as a guide in further exploratory investigations 
of successively more complex systems, which might have as their goal the 
formulation of a set (probably a large set, accompanied by massive tables 
of definite integrals) of empirical rules whereby an accurate configuration 
interaction wave function for any atom or molecule would be obtainable 
with a minimum of trial and error effort. 

Approximate Ground State Wave Functions and Energies..-We seek a 
normalized approximate wave function for the ground state of helium which 
has the form 


v= Da”, (1) 
where the V, are normalized antisymmetrized product wave functions 
(AP's) for various configurations. The best set of coefficients a, is de- 
termined by minimization of the energy E = Jf ¥*//¥dv, where H is the 
Hamiltonian operator for the helium atom neglecting spin-orbit interaction 
and relativistic effects and assuming infinite mass for the nucleus, namely, 

H = —-Vi-V:—4 ry — 4/r2 + 2/ria, (2) 
where the unit of energy is the atomic unit e?/2a) = 13.603 e.v. The 
best coefficients satisfy the equations 

> aH — SE) = 0, co (3) 
and the corresponding best energy £ is the lowest root of the secular equa- 
tion 
lH.; — Syé| = 0, 


where 


Hy, = SVHVde and 8, = S¥iWade. (5) 


The one-electron functions used in the construction of the AP’s are 
taken as products of the form 
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(i) = (nlm)7"(ny), (6) 


where (nim,)“"' is the hydrogen-like atomic orbital of total quantum 
number n,, azimuthal quantum number /;, magnetic quantum number 
m, and effective charge Z,,, and (n,) is a spin function (either @ or 8). 

A configuration of an atom is defined by the specification of the total 
and azimuthal quantum numbers for each electron in the system. Be- 
longing to a given configuration there are in general many different AP’s 
corresponding to the many possible assignments of magnetic and spin 
quantum numbers. We could in equation (1) include all the AP’s for 
each configuration which we wish to superpose and proceed to solve 
equation (4) by brute force. Since we are neglecting spin-orbit inter- 
actions, however, we know in advance that the secular equation may be 
factored into a product of independent equations, each factor of which 
involves interactions only between those LS vector-coupled functions 
belonging to the several configurations which have the same eigenvalues 
for L*, S*, L,, S,.7 Thus in the LS coupling scheme the |S ground state 
wave function for helium, W, is a superposition of '\S wave functions, 
?,, belonging to the several configurations: 


Vv = Yb. (7) 
The #, for a given configuration are linear combinations of the W, for that 
configuration; the procedure for obtaining the former from the latter 1s 
well known.’ 

Let us designate the |S vector coupled function belonging to the con- 
figuration (1s)7(1s)* by &(1s“, 1s*"), that belonging to the configuration 
(2p)@(2p)”" by &(2p%, 2p%"), ete. The approximate wave functions 
for the ground state of the helium atom which we have investigated may 
then be written as follows: 

%, = y By oe ae mngP( MS”, Ns”), (S) 
n= 1 mm = 1 
Vy = ByP(1s”, 1x7) + Byb(2p%, 2p") + Byb(3d™, 3d”) + 
Byo(4f“, 4%), (9) 


bo = CyP(1s“', 17") + Cxb(2p*, 267) + Cb(3d™, 3d”) + 
Cy(4f, 4f%). (10) 


Table 1 gives numerical results of variational calculations using these 
functions. For a given function, one assumes values for the effective 
charges, computes the matrix elements //;,, and Sj, (this requires only 
well-known procedures’), solves equation (4) for the energy /, and repeats 
with new values for the effective charges until a minimum £: is obtained. 
Discussion..-A single AP (e.g., function A-1 of Table 1) can be a solution 
of the wave equation for an atom only if the 7,, terms in the Hamiltonian 
operator are neglected or replaced by a spherically symmetric potential 
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function. The principal difficulty with such wave functions is that they 
allow the electrons to move independently of one another; hence, in 
general, the functions are not small, as they should be, when the electrons 
are close together. One says that the wave function does not correlate 
the motion of the electrons, and the difference between the correct ground- 


rABLE 1 


VARIATION FUNCTIONS FOR THE NORMAL HELIv! 


CONFIGURATIONS 
FUNCTION® SUPERPOSED PARAMETER VALUES 


A-l = (1s)? /, 6875 5.69531 0.1118 
b . 

(1s)?, (1s2s) Z, 757 5.7085 0981 

(1s)?, (1s2s), (1s8s) 2 = 1.777 5.7117 0049 

(1s)?, (1s2s), (2s)? /; 78 5.7095 0971 

(1s)?, (1s2s), (1s3s), Z TtH4 5.7148 O918 


(2s)?, (283s), (33)? 


(1s)?, (2p)? Z, = 1.6896, Z. = 3Z, 5.73870 O679 
(1s)?, (2p), (3d)? Z; and Z; as above, Z 5.742867 0637 
62, 
(1s)?, (2p)?, (3d)?, 2Z,, Z2, Z3 as above, 2, = 743770 0628 
(4f)? 122, 4 
(Isls’) Z, = 1.19, Z;’ = 1.8352, 5.751316 0552 


(1Ists’), (2p)? Z, and Z,' as above, Z: = 5.790369 0162 
1.162, 
(Isls’), (2p), (8d)? Zi, 2\', 22 aS above, Z; = 5.794033 0125 
9.52, 
C- (1sis’), (2p)?, (3d)?, 21, 21’, Za, Zz, as above, 5.79486 0.0117 
(4f)? Z;, = 16.02; 
Hartree Fock@ ‘ 5.734 0.073 


“ Functions labeled A, B, C are of the types ¥4, Va, We of equations (8), (9) and (10), 
respectively 

” Kellner, G. W., Z. Physik, 44, 91 (1927) 

° Eckart, C., Phys. Rev., 36, 878 (1930) 

4 Bethe, H., Handbuch der Physik, Vol. XX1V/1, p. 324 

* One atomic unit = e?/2a) = 13.603 e.v. More strictly the values of Feaic. should 
be interpreted as in units of ¢e?/2an., where dye/ay = 0.999863, which corrects for the 
finite mass of the nucleus 
’ Kexpu. = —5.80656e?/2a)—-see Paschen, F., Sitsber. preuss. Akad. Wiss., 1929, p. 662 
* These values are calculated assuming those in the previous column are in units of 
2ay. Strictly the correction indicated in note e should be applied, which increases 
all Feate.-Eexpu. values by 0.0008 


e? 


state energy of a system and that calculated with the best function of 
this type (the Hartree-Fock function) is called the correlation energy. 
Thus the correlation energy for the helium atom is 0.073 atomic unit or 
0.99 e.-v. (cf. Table 1). 
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There are essentially two ways (other than explicit introduction of 
ry. in the wave function) in which electron correlation can be brought into 
an analytical wave function for helium: either the electrons may be 
separated into concentric spherical shells or the wave function may be 
given angular dependence (e.g., one electron may be concentrated along 
the x axis and one along the y axis). In Table 1, all functions in set A 
and function C-1 fall in the first category, all functions in set B in the second 
category, and functions C-2, C-3 and C-4 in both categories. Set A 
separates the electrons into the concentric so-called K, L and M shells, 
as the same effective charge is used in each of the hydrogen-like orbitals 
ls, 2s, 3s. Essentially the same effect is achieved in C-1 through variation 


of the charges in Is orbitals. It might appear as though such an effect 


TABLE 2 
ErrecrivE CHARGES: Brest VALUES VERSUS VALUES DETERMINED BY CRITERION OF 
MAXIMUM OVERLAP” 
Z\¢ Z;' 22 Za Zs 
Function B-4° Best“ 1. 6896 32) 62, 12Z, 
Max. overlap 1. 6875 (10/3)Z, 72, 122, 
Function C-4‘ Best 1.19 1.8352, 4.162, 9.52; 16.02, 
Max. overlap 1.19 1.8352, 4.572, 9 62; 16.62; 
® Criterion of maximum overlap: to determine the effective charge for a given high 
energy configuration wave function, maximize the three-dimensional overlap of the 
radial part of this function with the wave function for the lowest energy configuration. 
® See equation (9) and table 1 
“ See equation (10) and table 1 
“ Variation of the effective charges was not carried out as precisely with the B func- 
tions as with the C functions. 
* Best and maximum overlap Z; values in the case of the B function differ because 
the former is an adjusted best value for the four-term function B-4 while the latter is 
the best value for the one-term function B-1. Adjustment does not give appreciable 


improvement in the case of the C function 


is also operative in the B set of functions, but there it is minor. For 
although hydrogen-like orbitals for K, L, M and N shells are used, the 
effective charges are varied, and the best values as given in Table | result 
in almost maximum overlap of the radial portions of the higher energy 
configuration wave functions with the Is* wave function (see below). 
Both effects are operative in the functions C-2, C-83 and C-4. Here 
(1s*', 1s“") separates the electrons into two shells, and the angular correla- 
tion is obtained from 2p*, 3d*° and 4f° configurations, where again the 
radial portions of the higher energy functions have almost maximum 
overlap with the lower energy function. 

In Table 2 the variationally determined effective charges for the higher 
energy configurations in the type B and C functions are compared with 
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values determined by a maximum overlap criterion. The agreement is 
striking, particularly in view of the fact that the calculated energies are 
not very sensitive to small changes in the parameters near their best values.* 

The angular and radial correlation effects are almost independent of 
each other. Thus, the improvement resulting from superposition of 
2p*, 3d? and 4/* configurations is about the same in both B and C sets of 
functions, even though the radial correlation obtained with B-1 and C-1 
are very different (zero in the former case). Successive improvements 
from introduction of angular factors in the two cases are, in atomic units, 
0.0434, 0.0042, 0.0009 for set B and 0.0390, 0.0037 and 0.0008 for set C. 
On this basis the angular correlation appears to account for about the 
fraction 0.045/0.073 or 62% of the total correlation energy, radial correla- 
tion for about 38°. The best function in Table 1, C-4, contains the 
fraction 0.0643/0.073 or SS% of the total correlation energy, about 64% 
of the radial correlation energy, and about 97°) of the angular correlation 
energy. For this function, C):C2:C3:C4:: 1000: —60: — 12: —4. 

Summary. The wave function for the ground state of the helium atom 
is approximated by a superposition of vector-coupled 'S functions for the 
configurations (1s, 1s“), (2p7)2, (3d”)?, (4d%*)?, where the one-electron 
functions are hydrogen-like orbitals with the indicated charges. The 
best calculated energy, 5.79486 e?/2a) (experimental value: 5.80656 
e?/2a 9), obtained with Z,; = 1.19, Z,;' = 1.835Z,, Z. = 4.162Z,, Z; = 
9.52, Z,; = 16.02, contains S89, of the total correlation energy (Hartree- 
Fock energy-experimental energy). 

Comparison with similar calculations putting Z; = Z,’ suggests that 
total electronic correlation energy in an atom may be considered as a 
sum of radial and angular correlation energies. Radtal correlation energy 
is tentatively defined as the difference between the Hartree-Fock energy 
and the energy calculated using the best superposition of all configurations 
obtainable from a complete set of purely radial one-electron functions 
(e.g., all ms hydrogen-like functions), and angular correlation energy is 
defined as the difference between total and radial correlation energies. 
The calculations then indicate that the correlation energy in the helium 
atom is about 62% angular and 38% radial, and that the best wave func- 
tion obtained contains 97% of the angular and 64° of the radial correlation. 

In attempting to regain radial correlation (separate electrons into 
concentric shells) weakly overlapping shells are found to effect little 
improvement; hence the use of the complete subset of ms functions with 
only one charge would be a slowly convergent process. Similarly in re 


gaining angular correlation, maximizing the overlap of the radial factors 
of the high-energy angle-dependent terms in the wave function with the 
purely radial first term is found to yield almost optimum values for the 


charges. 
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* This paper was prepared in connection with a conference on ‘‘Quantum-Mechanical 
Methods in Valence Theory” held at Shelter Island, September 7-10, 1951, under the 
auspices of the NATIONAL ACADEMY OF SCIENCES (see “Summary of the Conference’’ 
in a subsequent issue of these PROCEEDINGS). The paper is based on part of a thesis 
submitted by G. Russell Taylor in partial fulfillment of the requirements for the degree 
of Doctor of Science, Carnegie Institute of Technology, 1951. The work was assisted 
in part by the Office of Naval Research, U. S. Navy, through Contract Nonr-493(00) 
with Carnegie Institute of Technology, in part by a grant from the Research Corpora- 
tion. G. F. Hadley aided in the computations 

t DuPont Postgraduate Fellow, 1950-1951; Research Corporation Fellow, 1951. 
Present address: Mellon Institute of Industrial Research, Pittsburgh, Pennsylvania 

1 Hylleraas, E., Z. Physik, 65, 209 (1980) 

2 James, H. M., and Coolidge, A. S., Phys. Rev., 49, 688 (1936). 

3 For complex organic molecules, see, for example, Parr, Craig and Ross, J. Chem 
Phys., 18, 1561 (1950) or Coulson, Craig and Jacobs, Proc. Roy. Soc., 206A, 297 (1951). 

‘For important developments in theory, clear statements of the general method of 
attack and a calculation on the beryllium atom, see Boys, S. F., Proc. Roy. Soc., 200A, 
542 (1950); 201A, 125 (1950). 

® This has long been known and is implicit in statements of many authors, yet 
correct statements still crop up in the literature. Thus the LCAO molecular orbital 
method including configuration interaction is identical with the valence bond method 
including ionic structures, and as the number of AO’s is increased both methods converge 
to the true wave function 

§ Coulson, Craig and Jacobs, op. cit. 

7 Condon, E. U., and Shortley, G. H., Theory of Atomic Spectra, Cambridge University 
Press, Cambridge, 1935; Ufford, C. W., and Shortley, G. H., Phys. Rev., 42, 167 (1932). 

* Only after the present work was complete, indeed only after this paragraph was 
written, was it noticed that S. F. Boys (op. cit., paper I1) had earlier also come to the 
conclusion that maximum overlap of radial functions gives maximum gain of angular 
correlation energy. 


A COMPARATIVE SURVEY OF APPROXIMATE GROUND STATE 
WAVE FUNCTIONS OF HELIUM ATOM AND HYDROGEN 
MOLECULE* 


By R.S. MULLIKEN 
DEPARTMENT OF Puysics, UNIVERSITY OF CHICAGO, CHICAGO 37, ILLINOIS 
Communicated January 7, 1952 


For both atoms and molecules, there is a great need to find the simplest 
and most rapidly convergent succession of steps to improve on the familiar 
approximate wave functions of electron configuration or SCF (self-consist- 
ent-field) type. This may be done, in principle, by using wave functions 
depending explicitly on interelectronic distances 7,; (but this is cumber- 


. . , . . . . . - 
some in practice), or by forming linear combinations of wave functions of 
the electron configuration type (configuration interaction). In following 
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the latter path, Boys has indicated! that it is best to use generalized forms 
of electron configuration type, rather than forms which could, by them- 
selves, be SCF functions. In the present Note, a comparative survey 1s 
made in the light of this idea of various approximate forms which have 
been studied for the wave functions of the ground states of the He atom 
and the Hz molecule.' The survey also throws some light on the directions 
in which further studies might be aimed. 

Using the 'S, helium atom ground state as an example, assumption of 
an open-shell electron configuration (s‘s") permits better energy mini- 
mization than the usual closed-shell configuration (ls)*. Leaving aside 
spin and normalization factors (as will be done everywhere below), the 
corresponding wave functions are: 


For (1s)*: LsP ps: for (s’s"): gs’ SVs" \? 4 gs Dee CD (1) 


where Is, s’, s” all refer to AO’s (atomic orbitals) of s symmetry. The 
best (1s)? function is that using Hartree SCF ls AO’'s. The best (s’s”) 
function would involve two independent but not orthogonal forms s’ and 
s”, to be determined by a generalization of the SCF procedure. The 
superscripts | and 2 in Eq. (1) refer to the two electrons. //ere and below, 
symbols of the types («)* and (wy) are used as conventional abbreviations 
for electron configurations and also especially for (wo-electron wave func- 
tions having forms such as are illustrated by Eq. (1). 


The simplest approximations to s’ and s” are ls hydrogenic AO’s, Is 
t . ; eae 
and Is'', as used by Eckart and by Hylleraas, with effective Z values 
Z’ and Z” determined variationally. Somewhat better would be: 


s' =~ Ist + Alstt. 5” = istt + nist: (2) 


a little further improvement would be obtained by replacing the variation 
parameter \ in s” by an independent parameter. With s’ and s” chosen 
as in equations (2), the function (s’s”) can be written: 
(s’s”) = (Ististl) + [2n/(1 + v2] [Cst)? + (ist). (3) 

This formulation is of interest not particularly per se, but in its analogy 
to certain approximate expressions for the wave function of H, (see below). 

If we regard (ls)? as the ‘“‘true’’ ground-state electron configuration for 
helium, then the use of a better function of the type (s’s") (with s’ and 
s” two s functions) is a first step in taking “configuration interaction” or 
“correlation’’ into account. It corresponds to partly taking into account 
interaction with excited states of the types (ms, ns) and (ms)*. To take 
these fully into account, one might set up a series of mutually orthogonal 
forms of the type (s’s”) and seek to determine these and their coefficients.” 
After that, further improvements, converging toward the exact wave 
function, would be sought in the addition of 'S, forms of the types (p’p”), 
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(d'd"), etc. compare the paper of G. R. Taylor and R. G. Parr! elsewhere 
in this issue for part of these. 

If we regard (s’s”) itself as a generalized SCF form, its use may be said 
to reduce considerably the magnitude of the correlation energy or con- 
figuration interaction. The improved correlation in (s’s”) is accomplished 
by letting the two electrons move in considerably different AO’s, thus 
increasing their mean distance. This evidently more than outweighs the 
unfavorable exchange repulsion effect involved in the symmetric form 
which (s’s”) possesses (cf. equation (1)). 

Using the MO method, a completely analogous treatment can be given 
for the '}> +, ground state of the hydrogen molecule. The ordinary SCF 
closed-shell configuration is (1¢,)’, where lo, denotes a BMO (best MO). 
But the open-shell configuration (0,’ ¢,”), equivalent to partly considering 
configuration interaction with excited states of types (0,0,) and (a,)’, 
must give a considerably better wave function.* The next step would 
be to admix (la,)*, or better (o,'0,”"), and so on. (It is of interest that 
the best Io, has nearly the form of a united-atom hydrogenic 2po, although 
commonly it is approximated by the LCAO form (1s, — 1s)). Finally, 
other '3> +, MO-product forms of the types (x,,’1,"), (1,’r,”), (6,'6,”), and 
so on, should be added in due proportions.‘ 

Using the AO method, the ')>*, ground state of Hy has an electron 
configuration (0¢) of open-shell type; in rough approximation, 0 and @ 
are Z = | AO’s 1s, and 1s, of the two H atoms (Heitler-London method). 
It is seen that use of either of the configurations (¢,’0,") or (0%) for He 
may be regarded as a way of improving on (lo,)* by the same type of shell- 
opening as was used in equation (1) for He. And just as in the latter 
case, the possibility of improved correlation here is realized by letting 
the two electrons move in two considerably different orbitals, with a net 
gain in stability over even BMO (lo,)? if the two orbitals are chosen in 
the best way, in spite of the per se unfavorable exchange repulsion effect 
involved in the symmetric form (cf. equation (1) for (s’s”)) of the wave 
function.» The (0%) configuration for H», however, differs from (s’s”) 
for He and (¢,'0,") for Hz in that, regarded as MO’s, the orbitals o and @ 
do not belong to irreducible representations of the molecular symmetry 
group; they are nevertheless expressible as simple linear combinations of 
proper MO’s of the two species o, and oy. 

For the best (¢%) functions, PAO’s (polarized AO’s) may be used for 
the axially-symmetric orbitals o and %, with @ always the mirror image of 
o on reflection in the mid-plane perpendicular to the molecular axis. To 
approximate o and @ analytically, a possible first step is to add hydrogenic 
2po, and 2po, (with positive ends directed toward the center of the mole- 
cule) to ls, and 1s), respectively, with separate adjustments of effective 
Z values for ls and 2pe so as to minimize the energy. More generally, 
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o = Is.” + a2pe,'*” + B3do,'2” + .... (4) 


Using the expansion (4) or some equivalent form, one might proceed to 
determine BPAO’s o and @ using the variation method. This should 
lead to a wave function distinctly better than the well-known Rosen 
function. Next, additional terms of the same kind as (0%) but orthogonal 
to the latter might be added, analogous to what was suggested for He. 
To improve the wave function further, the next step would be to add 
functions of the type (wr), '}>+,. These were (implicitly) included by 
Hirschfelder and Linnett in their recent paper on H.? 

An interesting and instructive alternative to equation (4), although not 
leading so far toward BPAO’s, is® 


a = Isa + AIH; & = I? + Alsy. (5) 
The approximations of equations (5) for o and @ are formally analogous 


to those of equations (2) for s’ and s” of He. From equations (5), using 
for (o&) the analog of equation (1) for (s’s”), there follows: 


(08) = (1sq%15)%) + [24/(1 + A*)][(1547)? + (15p7)?]. (6) 


This expression, which is formally analogous to equation (3) for He, is 
seen to be the familiar Heitler-London-Wang-plus-ionic approximation 
for H». From what has been said above, it is clear that this cannot be 
as good as (0%) with BPAO’s for o and @. As is well known, equation 
(6) can also be recast in the form 


(1o,)*tcaccz) + Ylou)*Leaoz)» (7) 


which is simply the (Z-adjusted) LCAO approximation to (1¢,)* improved 
by taking into account configuration interaction with Z-adjusted LCAO- 
approximate (l¢,)?.° 

It will be noted that the PAO’s o and @ of equations (5) both reduce 
to the familiar LCAO approximation to the lo, MO if X = 1. In fact, 
Coulson and Fischer,’ arriving at these PAO’s from the LCAO-MO point 
of view, call them unsymmetrical MO’s! 

An interesting and characteristic feature for all the open-shell sym- 
metric forms (s’s”), (o,’0,"), (o&) is that when the two orbitals used are 
chosen in the best way, so as to give minimum energy, they are not orthog- 
onal. For example, as Slater has shown,” in the case of (0%) fot H2, 
if Lowdin’s OAO’s (orthogonal AO’s) are used for o and @, the energy is 
much worse than if Is, and 1s, are used, or than for (lo,)?. In fact, the 
change from natural AO’s 1s, and 1s, to BPAO’s for o and @ is in the direc- 
tion of increased non-orthogonality (increased overlap). In all of these 
cases, however, one could make use of orthogonal orbitals in various ways; 
for example, one could write 
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ma . ” a= ” . 
iC i Sty S — iC; diy 
_ i ae S oe ao P 
wa ili Tgiy Gy — iGi Fyiy 
= re ” ‘ > “ae , = 
¢ 2 >? Gi + 5 jCj Oxuj; F = > wr Cn = wre Cujs (S) 


where in each case the forms s,, 9,;, o,; are members of a convenient com- 
plete orthogonal set of orbitals of the indicated symmetry, with the c’’s 


” 


suitably different from the c”’s. Or, for example, supposing s’ known, 


* + 


s” may be written as s’ + as*, where s* is orthogonal to s’; @ and the 
particular form of s* (and of s’) must be determined by variation methods. 
The viewpoints presented above can be summarized and somewhat 


extended by means of the following table. 


* This paper was prepared in connection with a conference on ‘Quantum-Mechanical 
Methods in Valence Theory’ held at Shelter Island, September 7-10, 1951, under the 
auspices of the NATIONAL ACADEMY OF SCIENCES (see “Summary of the Conference” 
in a subsequent issue of these PROCEEDINGS) 

' This survey Note was inspired by pre-conference reading, including especially 
S. F. Boys, Proc. Roy. Soc., 201A, 125 (1950), Kotani’s and Taylor and Parr’s abstracts 
in the conference program, and references given below. Kotani and especially Lennard- 
Jones at the conference presented similar viewpoints as to Hy, im connection with reports 
on extensive and interesting computations. See also T. Inui, Proc. Phys.-Math. Soc. 
Japan, Sec. 3, 20, 770 (1938) 

2 The form (s’s") may be rewritten in terms of two mutually orthogonal AO’s s,; and 
s, in place of the non-orthogonal s’ and s”, as follows. Let 


so = a's, + b's; 5" = a”sy + b" So, (A) 


” 


where, if s’, 5”, s;, S) are now all taken as normalized, 


a’? + bb’? = landa”? + 5”? = 1 


7 


Then, for normalized (s‘s"), 
(s’s”) = 2a’a”N (s,)? + 2b’b"N (s.)? + (a’b” + a”b')N (s\83), (B) 
N = [201 + Q2)|'2; Q= JS s's" dr, (C) 


N being the factor which normalizes the form (s’s”) of (1), and the forms (s,;)%, (s,)? 
and (s,s.) being defined in the same manner as in (1). Equation (B) takes various 
forms depending on how a’ ora” is chosen in (A). Of interest are the choices a’ = 1 
(hence b’ = 0) and a” = 1 (hence b” = 0); or especially the choice a’b" + a"b’ = 0. 
The latter gives 

(s’s”) = N{(1 + OV(s.)? — (1 — O(s,)7}, (D) 


with 
s, = [211 — Q)]~'/%(s’ — 5”); 52 = [201 + Q)|~'/Xs’ + 5”) (E) 

Because of this equivalence of the forms (s‘s”) and a(s,;)? + 8(s,)*, the same final results 
should be obtained if one minimizes the energy by varying the latter form based on two 
orthogonal AO’s, or the former based on two non-orthogonal AO’s 

3 Hylleraas, if I remember correctly, once tried such a function, with an H,* BMO 
(or nearly so) for ag’, and a similar but lower-Z MO for o,”. 

‘The forms of (r’r”) and the like are somewhat more complicated than that of 
(o’0”). To obtain a form of '2*y type, one must take 
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(w'n”) = [f(1)g(2) + f(2)g(1)] cos (¢1 — ¢2), 


where x’ = fe'?; x” = ge'”. 

5 Usually the stability in the AO method of the symmetric form (1sq 1s»), as com- 
pared with Is, (1) 1s» (2) alone, is said to be due to the occurrence of a negative exchange 
integral instead of the usual positive integral found for atoms. A detailed analysis 
reveals, however, that the apparent negative exchange energy in (1sq 1s») is the sum of 
a positive exchange energy as for atoms, plus other numerically larger negative terms 
of the same nature as the resonance terms which make H,* stable as compared with 
H plus H* (cf. R. S. Mulliken, Chem. Rev., 9, 354 (1931)). 

* An interesting short-cut toward BPAO’s is that proposed by Gurnee and Magee 
(J. Chem. Phys., 18, 142 (1950)), who took o as 1s”) centered about a point somewhat 
displaced from nucleus a toward nucleus bd. 

7 Hirschfelder, J. O., and Linnett, J. W., J. Chem. Phys., 18, 130 (1950). 

* Coulson, C. A., and Fischer, I., Phil. Mag. [7], 40, 386 (1949); also Slater, J. C., 
Phys. Rev., 82, 538 (1951). 

® The amount of configuration interaction, i.e., the value of y, is of course considerably 
greater here than if BMO’s had been used for lo, and lo,; the same sort of thing is 
always expected when wave functions are not based on BAO’s or BMO’s. 

Slater, J. C., J. Chem. Phys., 19, 220 (1951) 


THE RADIATION SENSITIVITY OF CATALASE AS A:\- FUNCTION 
OF TEMPERATURE 
By R. B. SETLOw 
Bropuysics Division, YALE UNiversity,* NEw HAVEN, CONNECTICUT 
Communicated by G. E. Hutchinson; read before the Academy, November 5, 1951 


Introduction.—_A complete solution to the problem of the direct action 
of radiation on large molecules necessitates the knowledge of three things: 
(1) the amount and form of energy lost in passing through the material be- 
ing investigated, (2) the manner in which this energy is dissipated in the 
material and (3) the response of the system to the energy put into it. If 
one uses a type of radiation whose physical effects are well known, then the 
sensitivity of large molecules to this radiation can give information about 
the structure of the molecule. If the radiation employed is a beam of 
fast deuterons, we can be reasonably certain of the first item.' The man- 
ner in which the energy released in a molecule is dissipated is not well 
known, but an extension of the present work should shed additional light 
on this problem. We will be primarily concerned here with the response 
of an aggregate of catalase molecules to the energy released by the passage 
of deuterons through them. Charged particle irradiation of large mole- 
cules with specific functions usually gives the result that one primary ioni- 
zation inactivates the entire molecule.” * However, previous work on the 
irradiation of catalase by deuterons‘ indicated that an ionization within 
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the catalase molecule did not destroy all its catalytic activity. Some vi- 
ruses show results qualitatively similar to those obtained for catalase, and 
the usual reason given for the failure of an ionization to inactivate them is 
that the released energy only destroys a part of the structure, and the rest 


is able to function by itself.' In the case of T-1 bacteriophage it has been 


shown that the radiation sensitivity increases as the bombardment tem- 
perature increases.° The present work reports somewhat similar results 
for the irradiation of catalase by deuterons. ‘The results will be analyzed 
in terms of structure of the catalase molecule. 

Experimental Procedure.—In order to eliminate any indirect effects due 
to a suspending medium, all the irradiations were carried out upon dry 
catalase. The catalase used was derived from beef red cells. + Samples for 
bombardment were prepared by pipetting 0.05 ml., containing amounts 
up to 100 ug. solids, onto small 
glass coverslips. The coverslips 
were dried in a vacuum desiccator 
and then placed around the periph 
ery of a circular copper disk 
which was put im an insulated 
bombardment chamber of the cyclo Beam ood 
tron in which the pressure was 
<10°-4 mm. Hg. The disk could Shutter 
be rotated from outside the chamber 





so as to expose different samples ~— 
to the deuteron beam. The cur Tinch 
rent carried by the deuteron beam FIGURE 1 

was measured by a galvanometer. ‘tis saeiiaeiiiiiek haieiiads wiibioies 
The area of the beam was measured vision for raising the temperature during 
from the blackening of photographic — irradiation 

film placed in the beam, and the 

time of exposure was controlled by an electrically operated shutter. 
Further details may be found elsewhere.* 

For irradiations above room temperature, the bombardment disk was 
heated electrically, as shown in figure 1. ‘The temperature during the bom- 
bardment was measured by a calibrated chromel-alumel thermocouple. 
The temperature was kept within +2°C. of the average temperature. For 
bombardments below room temperature a device similar to that shown in 
figure | was used, except that the electrical heater was absent. The in- 
side of the bombardment disk could then be filled with liquids at different 
temperatures, such as liquid air or a water-ice mixture. 

After bombardment the samples were removed from the cyclotron, the 
catalase dissolved off the coverslips in water and assayed for enzymatic 
activity by comparing the amount of hydrogen peroxide decomposed in a 
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known time with a previously constructed calibration curve. The activity 
of control samples, which underwent the same treatment except for irra- 
diation, was also determined. The hydrogen peroxide was determined either 
by the amount of I, it formed from KI‘ or from its absorption at 247my 

as found in a Beckman spectrophotometer.° 
If the samples are deposited on quartz disks, it is possible to determine 
how much material is left on the coverslip after immersion in water by 
measuring the absorption of the disk in the ultra-violet region of the spec- 
trum. Data on irradiated and con- 





100 scons. seca al trol samples indicated that at least 
90) per cent of the material was dis- 
solved off. 

The energy of the deuterons used 
was 3.8 m. e. v. It is possible to 
calculate from published data’ that 
these deuterons, in passing through 
protein, will lose on the average 235 
e. v./100A. It is known that fast 
charged particles can lose their 


$s Ss 


Activity 


$ 


Relative 


energy in three ways: by ioniza- 
tion, excitation and nuclear recoils. 
For 3.8-m. e. v. deuterons only the 
first two are of significance, and 
because of the low quantum yield 
for excitation,® we shall neglect the 
1 : ; : ——J second and consider that the impor- 
2 4 6 8 +1 i mg 
Devterons per severe centimeter tant energy loss item is that due to 
ionization. It is known from analy- 
sis of energy loss by charged par- 
ticles in gases that the average 











FIGURE 2 


A semilog plot of the relative activity of 


catalase vs. the number of incident deuter- AGS , , 
energy lost per ionization is about 


110e. v. We shall assume the same 
value holds for solids.'- This primary 
ionization can of course give rise to energetic electrons which can pro- 


ons/em.? for several irradiation tempera 
tures 


duce ionizations of their own. 

Results. Some of the experimental results are shown in figure 2. The 
relative activity of catalase is plotted on a logarithmic scale against the 
number of deuterons per cm.*. It is seen that the room temperature curve 
is a straight line. This indicates that the relative activity may be given 
by an expression of the form . =e >”, where Ao is the activity of the 


41) 


controls, A the activity of the samples irradiated by a dose of D deuterons 


per cm.’ and S is a proportionality constant which we shall call the inac- 
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tivation cross-section. This expression of course simply indicates that the 
loss of activity under irradiation is a random affair, just as is radioactive 
decay. If we want to be literal about it, we can give the proportionality 
constant S physical significance. S has the units of an area, and thus it 
can represent the cross-section of a volume which is inactivated by the 
passage of a deuteron through this volume. For example, from the data 
taken at room temperature, one can compute that the cross-section cor- 
rane ate ; asd 
responds to 1420 A*. This is seen by finding the dose for which yey ” ie 
0 

= 0.368. The cross-section is then simply the reciprocal of this dose. 
Since catalase has a molecular weight of 225,000,° and thus a cross-sectional 
area of about 6000 A’, the inactivation cross-section is obviously too small. 
But a deuteron can pass through a catalase molecule without losing any 
energy at all, since the distribution of primary ionizations along a deuteron 
path is a random one. The chance of such a miss may be found by using 
bombarding particles of different velocities.?. Previous work with different 
energy deuterons and with 2-m. e. v. electrons,’ which because of their 
very sparse ionization give an inactivation volume instead of cross-section, 
indicates that a cross-section of 1420 A® corresponds to a volume whose 
molecular weight (assuming a density of 1.33 g./em.*) is about 110,000. 
Thus if we assume catalase has a molecular weight of 225,000, then a pri- 
mary ionization at room temperature knocks out one-half of its activity. 
It is seen in figure 2 that irradiation at liquid-air temperature yields a 
significantly smaller cross-section than at room temperature. 


Data taken at 43°C. give a more complex curve than at room tempera- 


ture. This curve, however, can be resolved into at least two components, 
one of which has the same cross-section as the room temperature curve. 
The rapid decline in activity at low doses has been associated with a partial 
heat inactivation of catalase.‘ It should be pointed out, however, that 
dry catalase is not appreciably inactivated for several days at a temper- 
ature of 43°C. Since the reasons for this initial rapid decline in some of 
the curves obtained are not well understood at present, we shall neglect 
these data and concentrate on the part with the small cross-section. Other 
typical data are shown in figure 3. While cross-sections from these data 
may be obtained with high precision, it is only fair to recognize that the 
absolute accuracy is nowhere nearly as good. The major source of this 
error lies in the possible lack of uniformity in the cross-section of the deu- 
teron beam and in the calibration of the beam galvanometer. 

Figure 4 is a compilation of data obtained from runs similar to those 
already presented. Each point represents one run except for 300° and 
90°K. which represent averages of five runs each. The general trend 
is for the cross-section to increase as the temperature goes up. This in- 
crease, however, is not by any means uniform. There definitely seem to 
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FIGURE 3 


A semilog plot of the relative activity of catalase vs. the number of incident deuter- 
\ log plot of tl lat ty of catal tl I f lent deut 


ons/em.? for several irradiation temperatures. Note the similar slopes for 60.5° and 
93°C. and the big difference between 104° and 108.5°C. The points of the 104°C. run 


which were at a higher temperature are indicated. 
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FIGURE 4 


The inactivation cross-section for 3.8-m. e. v. deuterons vs. the temperature during 


irradiation, The temperature scale is not uniform. 
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be temperature intervals over which the cross-section remains constant. 
Somewhat similar results have been found for invertase.'' At the highest 
temperatures indicated the rate of heat inactivation is appreciable. At 
385°K. the specific rate constant for inactivation is about 2 X 10~* sec. ~!.4 

Discussion. —In the analysis that follows we shall assume that the molec- 
ular weight of catalase is 225,000. As mentioned previously, we can take 
other data to indicate that a cross-section of 1420 A* corresponds to the 
inactivation of a molecular weight of 110,000, or just about one-half the 
molecule, by a primary ionization. 

The only thing that we are varying in this experiment is the manner of 
degradation of the energy released within the catalase molecule by an ion- 
ization. Since the average energy released for all the temperatures in- 
volved is the same, the relative cross-sections at different temperatures 
are measures of the sensitivity of the molecule at these different tempera- 
tures. Because a molecule is a three-dimensional structure, it is more 
significant to calculate the volumes associated with these cross-sections, 
Since the volume depends upon the 3/2 power of the area, it is a simple 
matter to find the volume inactivated by a primary ionization at various 
temperatures, from the experimental observation that 1420 A? cor- 
responds to a molecular weight of 110,000 inactivated per primary ioniza- 
tion. The results have been put in terms of molecular weights since these 
are more recognizable than volumes. We thus obtain the values pre- 
sented in table 1. 


FABLE 1 
THe MoLecuLarR WEIGHT INACTIVATED PER PRIMARY IONIZATION FOR DIFFERENT 
‘TEMPERATURES 
MOLECULAR WEIGHT MOLECULAR FRACTION 
PER PRIMARY PER PRIMARY 
TEMPERATURE, °K CROSS-SECTION, A? IONIZATION IONIZATION 
300 1420 110,000 
90 925 58,000 
350 2450 250,000 
385 4200 500,000 


Qualitatively there is a simple explanation for these results. At liquid 
air temperatures all the bonds in the molecules are in their ground vibra- 
tional states. Any large amount of localized energy introduced into the 
molecule will be degraded to thermal energy within a relatively small 
volume. As the temperature is increased the energy will still fall off 
rapidly with distance from the site of introduction, but because some of the 
higher vibrational levels are now excited, it is easier for a new molecular 
configuration to arise. The data indicate that the catalase molecule can 


function in parts. If this were not so, the cross-section versus temperature 


curve would show a steady increase with temperature. Since catalase has 
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four heme groups, the correlation between the molecular fraction inactivated 
and parts of the molecule is excellent. It is of interest that there is no 
evidence for a cross-section plateau corresponding to three-quarters of the 
molecule. This may mean that a quarter catalase molecule is inactive, 
as has been suggested,'’ or that inactivation of three-quarters of the mole- 
cule alters it so drastically that its entire structure is changed, and it no 
longer is a closely knit entity. At temperatures above 380°K. the cross- 
section rises above that corresponding to the whole molecule. At these 
temperatures the heat inactivation rate is large, and any slight additional 
energy could cause inactivation of the entire structure. If this energy 
comes from a deuteron passing outside the molecule, owing, for example, to 


a secondary electron or an energetic excitation, then obviously one would 


measure a volume greater than the size of the molecule. 

Summary.—The sensitivity of catalase to irradiation by 4.8-m. e. v. 
deuterons depends on the temperature at which the bombardment is per- 
formed. From 90° to 385°K. three different regions of constant sensitiv- 
ity have been found. They correspond to inactivation of one-quarter, 
one-half and the entire enzymatic activity of the molecule. There is no 
evidence for the inactivation of three-quarters of the molecule. 
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Although invertase was one of the first enzymes discovered, many facts 
concerning its physical properties are still unknown. Conflicting figures 
are given for its molecular weight" * and data quoted on the thermal inac- 
tivation® do not take into account the partial reversibility of inactivation of 
invertase discovered by Herriott.‘. In this work we have observed the 
effect of heat on invertase in the dry and wet conditions, of fast electron 
bombardment, of deuteron bom- 
bardment and of combined heat and 





deuteron bombardment. In addi 
tion we have studied the effect of 
deuteron bombardment on solubility 
of the enzyme and on the serological 
activity as measured by combina- 
tion with specific antibody. 

The result has been to enable a 
measurement of the molecular weight 
of the unit of invertase which is 
enzymatically active, of the unit 
which is serologically active, and to 
show the relationship between ther- 
mal and ionizing radiation action 


PER CENT REMAINING 


on invertase regarded as a protein. 
Thermal Inactivation. The general 











stability of invertase to various r ‘ , ‘ . 
or 7 ‘4 26 40 ”) rr) 100 120 
forms of inactivation was measured TIME IN MINUTES 


roughly in the early days of enzy- FIGURE | 
mology.’ Very little study has been The thermal inactivation of wet inver 
made of dry inactivation. Herriott® tase at pH 4.5. There is clearly a rever 
showed that at pH 4.6 invertase is sible process occurring at low tempera 
reversibly inactivated when heated — tures 
in solution, though the reversibility 
applied apparently only to part of the enzymatically active material. 
Chase, Reppert and Ruch*® showed clearly that inactivation reached an 
equilibrium at 50°C. and pH 7. 

We have studied the effect of temperature at pH 4.5, pH 7 and in the 
dry state over a range of temperatures. The enzymatic activity was deter- 
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mined by the method of Sumner and Howell.’ The enzyme (melibiase 


free) was supplied by Nutritional Biochemical, Inc. 

Typical results of the study are shown in figure 1. This shows the wet 
inactivation at pH 4.5. In these curves the per cent activity remaining is 
plotted as a function of time at various temperatures. It can be seen that 
at low temperatures the inactivation proceeds to a definite fraction and 


thereafter remains steady. This is the feature discovered by Chase, 
Reppert and Ruch. At pH 7 and 41°C. the inactivation remained at 50% 
after seven hours of exposure. Such samples regain still more activity 
after cooling to lower temperatures but the process is not completely re- 
versible. This is in agreement with Herriott’s findings. 

If we treat the process as one of reversible inactivation so that we have 
a reaction constant ; for inactivation of native enzyme and k» for the re- 
covery from the inactive form the measurements allow values of k, and 
ky to be deduced from the curve. These are shown in table 1. 


TABLE 1 


TEMPERATURE, 
¢ ki, FRACTION/SEC ko, FRACTION /SEC Fei, CAL. /MOLE 


1% OS 7 10"! 4.14 10°4 24,500 
5D4.0 ‘ 10-¢ 6.8 XK 10~4 24,500 
59.0 x 1073 8.0 K 10°34 23 ,950 
61: 3.26 10~' 2.1 X 10 23 ,550 
64 75 X 10 23 ,350 
41: oH EE 4.26 X 10 23 ,550 
AW 9% . oo 3.90 X 10 23,150 
1): 56 10~* 3.90 XK 10 23,000 
19 4 ¢ 10 22,800 
54 9 X 107? 22,150 
59 3.20 % 107 21,900 
YO 5.8 X 10 30,250 
100 : C107 30, 100 
107 5.5: 107 29,900 
117 ( 1! 30, 100 
126 < 107" 30,200 
144 ay X 107° 30,250 


It can be seen that for both pH 4.5 and pH 7 &, the rate of the forward 
reaction increases markedly with temperature. On the other hand, k» 
first rises and then falls very sharply. For this reason studies of the re- 
activation of invertase above 65°C. fail to show any appreciable reactiva- 
tion. 

Analyzing the variation of k; with temperature in terms of the theory 
of absolute reaction rates** we can derive the free energy AF* of activa- 
tion. The resulting values are given in table 1 for the forward reaction. 
It can be seen that AF * diminishes with temperature. If we put 


AF* = AH* — T AS* 
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we deduce A//* = 53,000 cal /mole at pH 7 and AS* = 95 cal./mole per 
degree while at pH 4.5 the figures are AJ7* = 48,500 and AS* = 75 
cal./ mole. 

These are much lower than those quoted by Glasstone, Laidler and 
Eyring. We suggest that the discrepancy is due to the existence of the 
reverse recation at low temperatures. 

In the dry state the reaction constants are as shown in table 1. Invertase 
is very much more stable when dry. We have not found any conclusive 
evidence for reactivation in the dry 
state. The temperature variation 
fits a constant value of AF* at 
30,000 cal./mole. The entropy of 
activation is accordingly zero or very 
nearly so. This fits a general be- 
havior already noticed for some 


bacterial viruses and for two or three 


REMAINING 


enzymes for which it is found that 
in the dry state AS* is either zero 
or slightly negative." '! 

Inactivation of Dry Invertase by 
Tonizing Radiation. The fact that 


PER CENT 





ionizing radiation produces a number 
of separated but highly active 
agents the primary ionizations 

renders it possible to use this dis- 





tributed energy to form some esti- pF 


10 a: 
mates of the size and shape of DEUTERONS PER SQUARE CENTIMETER 
molecules which are sensitive to the FIGURE 2 
action of radiation. This was used The inactivation of invertase by deu 
by Svedberg and Brohult'? who  teron bombardment. The per cent re 
maining, both for enzymatic and serologi 


studied the splitting of hemocyanin eae en J 
by alpha articles in solution cal activity, is a logarithmic function of 
ees dane I ; ; : * the number of deuterons per square centi 
Further studies on ribonuclease and eter 
myosin by Lea, Smith, Holmes 

and Markham,'* indicated that very large fractions of the proteins 


were radiation sensitive. Deuteron and electron bombardment is more 


precise in that the location of the primary ionizations is definite. This 
was tried for pepsin and trypsin by Pollard, Buzzell, Jefferys and Forro'™ 
and reasonably accurate values for the molecular weights of these enzymes 
found. Setlow'® found a similar result to hold for urease. Catalase'® and 
hemoglobin” are apparently not wholly radiation sensitive under normal 
conditions, a fact which is linked to their multiple structure. 

Deuteron and electron bombardment as described by Pollard, Buzzell, 
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Jefferys and Forro'' was applied to invertase. It consisted in irradiating 
dry samples with known numbers of deuterons or electrons. We are in- 
debted to Mr. Marvin Slater for help with the electron irradiation. 

The results of deuteron bombardment on enzymatic action are shown in 
figure 2 (lower line). The relation In n/np = —SD (where n/n is the 


fraction remaining, S the apparent cross-section per deuteron and D the 


number of deuterons incident per square centimeter) is obeyed. S should 
be related to the molecular area if the deuteron beam is 100% efficient in 
inactivation. 

To test this, the energy of the deuteron beam was varied and S measured. 
For very slow deuterons the ionization per unit path is much greater than 
for fast deuterons and so a greater efficiency should result. The findings 
are plotted in figure 3 where the observed S is seen to increase with ionization 
density to a value of 5.1 X 
10-'* em.*. Some correc- 
tion has to be applied for 
the fact that at low veloci- 
ties deuterons more readily 
produce delta rays. The 
proportion of these can be 
estimated from theory'* and 
have been tabulated for 
deuterons by Lea;'® the esti- 
mated correction is shown 
200300400800 —“600 700 by the crosses. The area 
ENEROY LOSS PER 10K ow of the sensitive part of the 

FIGURE 3 molecule is 4.0 KX 107 em.?. 

The apparent target cross-section for deuteron The dotted curve is derived 
bombardment measured for various energy deu- * from simple theory on the 
terons with a variety of values of the energy loss per basis that the thickness of 


100 A. The crosses are corrected for delta ray ‘pee 
sae ; se the molecule is 50 A and the 
effects and the dashed line is a theoretical efficiency | 
curve assuming the invertase molecule to be 50 A fit is reasonably good. 
thick The above measurement 


gives information on the 

effective area of the molecule. If ionization which is random in volume is 
used the inactivation curve takes the form In n nm) = — VI, where J is the 
number of primary ionizations per unit volume and I’ 1s the target volume. 
Such volume-random ionization is produced by fast electrons. Accord- 
ingly dry enzyme preparations were bombarded by fast electrons from a 
linear accelerator. The appropriate logarithmic inactivation was found 
and the mean value for the inactivation volume was 1.52 & 10~'%em.'. 
We thus have three pieces of information about the inactivation of 


« 
° 











invertase. 
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(a) The target area for perfect efficiency = 4.0 X 10°" em.? 
(b) The target volume is 1.5 K 107' em. 


(c) The effective thickness is 50 A. 


If it is assumed that the invertase molecule is a cylinder of radius r and 
length /, and that the molecules are lying flat when bombarded, the first 
two lead to a radius of 24.2 A and a length of 83 A. The third piece of 
information checks reasonably well with a diameter of 48.4 A. 

The use of ionizing radiation in this way leads to the idea that invertase 
is a somewhat elongated molecule. 

The molecular weight, assuming a density of 1.35 is 123,000. 

The only available data on the molecular weight of invertase are con- 
tained in the diffusion measurements by Euler, Josephson and Myrback' 
which when analyzed by the Einstein-Sutherland relation (a correction 
suggested by Dr. R. B. 
Setlow) lead to a value of — y.o%,' 
120,000 for the dry enzyme. 2 

Serological Assay. In " 
order to investigate the 
effect of deuterons on the 
ability of the molecule to 
combine with specific anti- 
body, injections of invertase 


were made into a rabbit, ion 


“ 


and specific antiserum pre- _———— a 
pared. We were assisted eee 
in this work by Mr. Wilbur Q “ae 
J. McNulty. The loss of 

ability to precipitate anti- FIGURE 4 

body is shown in the upper The apparent cross-section per deuteron for differ 
line of figure 2 and itisclear ent temperatures of the specimens bombarded. 
The enzyme sensitivity increases to a steady value 
at the cross-section of the molecule and thereafter 


S 


INVERTASE 


APPARENT CROSS SECTION 








wall, 


eeeieemias " Canaan 
TEMPERATURE (DEGRELS KELVIN) 


that the cross-section is less 
than for enzymatic inactiva- 
tion. The target area found 
is 1.25 K 10°-"cem.*. The 
estimated molecular weight of the fraction responsible for serological 
behavior is 27,000. 

One very important factor in protein inactivation is solubility. To 
check on this the amount of material which goes into solution after irra- 
diation was measured by weighing the amount recovered after solution and 
evaporation. Within the accuracy of the experiment all the invertase was 
soluble after bombardment, unless very heavy doses were given. 

Combined Effect of Temperature and Deuteron Bombardment... The re- 


rises sharply as excitations and photons become able 
to produce inactivation. 
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markable behavior of ionizing radiation im removing enzymatic activity 
from a volume as large as the measured volume of an invertase molecule 
requires some interpretation. In order to find some basis for this, the 
effect of keeping invertase preparations at a wide variety of temperatures 
while irradiating was studied. An effect of combined temperature and 
ionizing radiation on a bacterial virus in the dry state was found by Adams 
and Pollard” and an effect of heating before and after irradiation has been 
reported by Setlow for catalase.'® The results of runs taken from liquid 
air temperature to 100°C. are shown in figure 4. 

Two bombarding energies were used, representing the extremes possible 
in ionization density. ‘The series at 7 e. v. per A show three significant 
facts. ‘The first is that at low temperatures the apparent cross-section is 
less than that at room temperature. The second is the constant cross-sec- 
tion at 4.2 * 10~'em.* from 220° to 290°K., which represents the cross- 
section for the whole molecule and the third is the considerable upward 
sweep at high temperatures. 

Two runs with the electron accelerator at 300° and 100°K. showed that 
at 100°K. the inactivation volume is about 0.67 times that at 300°K. 

Discussion. In the dry state, for thermal inactivation, invertase has an 
activation energy of 30,000 cal./mole with no entropy of activation. 
This corresponds to the passage over a potential barrier of 1.3 electron 
volts. No reversibility is observed. This barrier can actually correspond 
to the breaking of one bond, which seems rather hard to believe in a cross- 
bonded protein, or to the breaking of several. Three hydrogen bonds are 
a possibility 

Wet inactivation differs in that the molecule is hydrated and therefore 
subject to the effect of electrostatic forces due to groups which can be ion- 
ized. The dissociation of these groups at low pH accounts for the varia- 
tion of the heat of activation. 

The action of ionizing radiation alone appears to be explainable in terms 
of the requirement that a primary ionization fall within the molecular vol- 
ume. This releases 110 e. v. and undoubtedly causes the rupture of many 
bonds. It nevertheless appears as though this effect, large as it may be, 
does not have a 100° chance of inactivation because at low temperatures 
there is a clear drop in efficiency. It is possible that the cause of this lies 
in the differential internal expansion of the protein. Unlike most crys- 
talline solids, proteins are structurally dependent on three different 
classes of bonds, those in the polypeptide chain, those linking the side 
chains in a helix and those linking residues between helices. The effect of 
lowering the temperature may be to cause stronger bonding between res- 
idues holding the helices in place. Raising the temperature would cause 
a weakening. 

On this view the inactivation at low temperatures requires more than 
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the normal energy release at a primary ionization, At higher temperatures 
where excitation, as opposed to ionization, is sufliciently energetic to cause 
inactivation the apparent cross-section rapidly rises 

The serological cross-section to deuteron bombardment is definitely lower 
than the enzymatic cross-section. Two explanations of this can be given. 
The first is that the serological properties are less specific and therefore 
that damage can be sustained without removing the ability to combine 
with antibody. This may, indeed be true but we are inclined to propose a 
second explanation in terms of the presence of a second non-enzymatically 


active protein in large excess. The invertase we used did not have either 


the high activity or the thermolability of preparations described by Sumner 
and O’Kane*! and so must contain some protective protein. The serological 
cross-section would indicate that this protem has a molecular weight of 
27,000. Judging from the low effect of deuteron bombardment on solu- 
bility this proteim must be highly soluble. 

Summary.—Thermal inactivations of invertase at pH 4.5 and pH 7 shows 
that a partly reversible loss of activity occurs with the reverse reaction 
constant passing through a maximum at 59° for pH 4.5 and 44° for pH 7. 
For the forward reaction the value of A//* is 53,000 cal./ mole at pH 7 and 
48,500 at pH 4.5. Dry invertase shows a different thermal behavior with 
no observed reversibility and with a value of A//* of 30,000 cal./mole. 
There is no entropy of activation in the dry state. 

Deuteron and electron bombardments of dry invertase fit with the idea 
that one primary ionization inactivates the whole molecule. The cal 
culated size of the molecule, assumed cylindrical, is 4S A in diameter and 
83 A long. The molecular weight is 123,000. 

Combined thermal and deuteron bombardment gives a varying cross 
section which shows a constancy over 100° K. at the molecular size but is 
smaller at low temperatures and much higher at high temperatures. 
The variable sensitivity is tentatively explained in terms of differential ex- 
pansion inside the molecule. 

Serological studies indicate that the invertase used also contains a 
protective protein of high solubility and molecular weight 27,000. 
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Communicated by P. R. Burkholder; read before the Academy, November 5, 1951 


Introduction. The type specific transformation of pneumococcus R 
cells to S cells by an extract of S cells was first demonstrated by Griffith! 
for intraperitoneal injection into mice, and later was demonstrated in vitro 
by Dawson and Sia.2- Avery, MacLeod and McCarty* showed that the 
active transforming principle is in the desoxyribosenucleic acid (DNA) 
fraction of the S cells. The DNA which they extracted was similar to 
DNA prepared from other sources and its molecular weight was therefore 
estimated to be about 500,000. Direct observation of the active DNA 
in the ultracentrifuge does not provide a determination of the molecular 
weight, since the boundary diffusion is insufficient to give the necessary 
information about particle shape. 

Our purpose is to report evidence consistent with a larger particle size 
for the transforming principle. Dry preparations of the active DNA have 





Vou. 38, 1952. BIOCHEMISTRY: FLUKE, DREW AND POLLARD 181 


been bombarded with 2-m. e. v. electrons and 3.8-m. e. v. deuterons with 
indication of a molecular weight about 6,000,000. 

The first attempt to determine a molecular weight for a biological mole- 
cule by ionizing radiation was by Svedberg and Brohult* and later Lea, 
Smith, Holmes and Markham? found molecular weights for ribonuclease 
and myosin. The technique has since been applied to numerous materials 
by various authors. Pollard and Forro® and Pollard, Buzzell, Jeffreys 
and Forro’ give a complete discussion of the technique as it has been 
applied in the work reported here. 

The process by which ionizing radiations indicate particle size or structure 
may be summarized briefly. These radiations in passing through material 
lose energy to the material in many local events, whose size and spatial 
distribution may be determined. If the material is dry, the action is 
direct and can affect the structure of the material only at the locales where 
the energy is released. Where the material embodies some _ biological 
principle, easily destroyed by even minor changes in its structure, the 
relation of inactivation of the principle to distribution of energy release 
gives an indication of the size and shape of the particles of the principle. 

Deuterons of the energy used in this work distribute ionization events 
densely along straight paths through the material, inactivating in ac- 
cordance with the relation N/N) = e~>”, where N/ Np is the fraction of 
activity remaining after subjection to a dose D deuterons/em.’. The 
deuterons thereby probe for the frontal area of sensitive regions of cross- 
section S cm.*/deuteron. 

Since the electrons used ionize much more sparsely and do not follow 


straight paths through material, the incident doses are translated into 
ionizations/cm.*, The ionizations are distributed randomly in volume, 
inactivating as N/Ny = e-'’, where / is the number of ionizations per 
unit volume and V is the volume in which at least one ionization must 


occur in order for inactivation to result. A molecular weight can be 
calculated for |’, knowing the density of the material. Comparison with 
the deuteron cross-section gives an indication of the shape of the particles. 

Experimental. The data for this work were obtained in eleven runs 
with the 4-m. e. v. cyclotron at Yale, two runs with the 2-m. e. v. Van de 
Graaff machine in the chemistry department at Brookhaven National 
Laboratory, and one run with the linear accelerator at Yale, also for 2- 
m. e. v. electrons. The bombardment chambers used with the cyclotron 
and with the Van de Graaff machine are described in reference 7. The 
chamber used with the linear accelerator was fully insulated for beam 
collection, was arranged to collect secondary electrons and opened directly 
to the accelerator vacuum. 

A preparation of active DNA diluted with buffered saline was pipetted in 
amounts from 5 to 105 micrograms of DNA onto '/»-inch glass cover slips, ar- 
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ranged in lots of eight or ten on bombardment disks. ‘The DNA specimens 


were dried in a refrigerated desiccator pumped on through a liquid air trap. 


A smooth strain, pneumococcus Type III (A66) was used as a source 
of DNA. The active principle was prepared according to the improved 
ruethod of McCarty and Avery.’ Strain R36A, a rough variant derived 
criginally from a virulent ‘‘s”’ culture of pneumococcus Type I] was used 
as the test organism for transforming activities remaining for the various 
DNA samples after exposure to ionizing radiation. 

The special environmental conditions necessary for transformation of the 
susceptible R cells have been discussed in detail by Avery, MacLeod and 
McCarty® and by Taylor.’ /n vitro transformation of the unencapsulated 
Kt strain to the encapsulated $ strain requires not only the addition of the 
DNA to the nutrient materials but the further addition of an agglutinating 
agent for the R cells, an accessory protein constituent of serum and sodium 
pyrophosphate. In the present study, anti-R36A rabbit serum was used 
as the agglutinating agent and the accessory protein constituent was sup- 
plied in the form of a purified fraction bovine albumin.’ It was found 
convenient to use a solution of albumin pyrophosphate which was prepared 
just prior to use by mixing one part of a (7,15 solution of sodium pyro- 
phosphate with ten parts of a 4°, albumin solution, Charcoal-adsorbed 
heart infusion broth prepared as described by MacLeod and Mirick" 
constituted the assay medium. 

After exposure to the beam the cover slips with the DNA specimens 
were dropped into tubes containing 0.5-ml. buffered saline This was 
done as soon as possible, so that the total time dry was usually not over 
4to 6 hours. Tubes containing 8.6 ml. of broth, 0.1 ml. of a 1:10 dilution 
of antiserum, 0.4 ml. of albumin-pyrophosphate solution and 0.25 ml. of 
a 104 dilution of a 6- to 7-hour blood broth culture of the test organism 
were set up. The specimen cover slips in their tubes were crushed into 
fine pieces using the tip of a |-ml. serological pipette and the saline solution 
of DNA was thoroughly remixed by repeatedly drawing it up and down 
in the pipette. The solutions were transferred to tubes of the assay 
medium and the resuspension tubes rinsed with 0.2 ml. of broth. After 
careful and uniform mixing of the DNA solutions with the assay medium, 
each mixture was distributed in 2-ml. amounts among five small tubes. 
The tests were incubated at 37°C. for iS hours and examined for clouding 
of the broth which indicates the growth of transformed cells. In each 
series of five tubes the number in which transformation occurred was 
noted. To confirm the results a loopful of culture from each tube was 
streaked onto a section of a blood agar plate. The plate cultures were 
examined after incubation at 37°C. for 16 hours for the appearance of 
glistening, mucoid colonies typical of pneumococcus Type ITT. 

The assay procedure, giving as the result some number of positive trans- 
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formations among five observations, is a sensitive indicator of the amount 


of transforming principle present only in the neighborhood of some critical 


amount, corresponding to two or three transformations out of five. 


The 


assay does not permit simply putting some DNA on a cover slip and after 


bombardment removing it to assay the per cent activity remaining. 


Instead, the same amount of DNA was put on several cover slips (usually 


four or five), which were then bom 
The 


assay result was then interpolated to 


barded with different doses. 
find a dose corresponding to survival 
of the critical amount, which when 
divided by the amount originally 
put on the slips gave the relative 
activity remaining. 

The critical amount was deter- 
mined for each bombardment by a 
control series subject to the same 
experimental conditions throughout, 
except for actual exposure to the 
beam. Cover slips with graduated 
amounts of DNA were fastened to 
the bombardment dises along with 
the specimens to be irradiated. —In- 
terpolation of their assay results de 


termined the critical amount applicable in that experiment. 
figure | shows the result of a typical run. 


were 12, 24, 42, 60 and - 
ETS bp ' 

105 micrograms of DNA, 80 

Taking three out of five 


positive as the arbitrary 


remaining 
i 
oO 


end-point, the critical 
indicated is 50 
The 


diated points all received 


unount 


micrograms. imra- 


% activity 


105 micrograms of DNA 
and were exposed to doses 
of 1.3, 2.4, 5.7, 7.7 and 18 
(all X 


em.*. Amounts and doses 





10'") deuterons 


Pots 


ff critical amount 
50 pe * Ope 


controls, no dose 


out of 5 


6 80 100 
ve@ DNA 


responses 


critical dose 
(6 5°25)x 10° deut Am? 
all points 105 pg 


positive 





24 6 6 10 12 4 16 18 20x10" 
deut /em? 


FIGURE 1 


‘Transformations out of five as a fune 
tion of the amount of DNA and of the 
dose, illustrating the method of deter 
mining per cent activity remaining 


Asan example, 


The control series amounts 


37% Dose = 7«10'° deut /cm? 


pst 


20* 10° deut /cm’ 





FIGURE 2 


Per cent activity remaining (log seale) as a function 


were usually chosen in 


approximately geometric 


of deuteron dose 


series, since the critical amount varied considerably from run to run and 


there was risk of losing the result by attempting to pull the data closer to a 
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predicted end-point. At three out of five, the same arbitrary end-point 
as for the control series, the interpolated dose in figure | is 7 XK 10" 
deuterons/em.’, corresponding to 50/105 micrograms, or approximately 
50% activity remaining. 

The data in figure 1 also illustrate the estimation of assay uncertainty, 
carried over into the dose-effect curves presented below. In the example 
the critical amount is deemed reliable only between 40 and 60 micrograms 
(40 to 60% activity remaining) and the dose between 4 and 9 X 10'° 
deuterons/em.?. On the dose-effect curves these uncertainties are ex- 
pressed as limits of error 

Results and Interpretation... The principal data of this experiment are 
collected in the next two figures. Figure 2 shows the per cent activity 
remaining as a function of dose in deuterons/em.*. The error brackets 

refer only to the uncer- 

0 T T tainties in interpreting 





the transformation assay. 
The two points which are 
not blackened are the 


37% dose 5 « 10" tlectrons 
cm2 


s 


result of runs at approxi- 
mately doubled ion den- 
sity (by placing foils in 
the beam), and it may be 
observed that there is no 
change in the cross-sec- 


> 
=] 


remaining 


& 
> 
® 20 
td 


: tion within the dispersion 
10 « 10 elec em? of the data. If the sen- 
sitive regions had a di- 
a ; mension less than about 

Per cent activity remaining (log scale) as a function : J 
of electron dose 60 to SO A in the beam 
direction (or randomly 
oriented), such a change should be detectable. A straight line has been 











FIGURE 3 


drawn through the data, corresponding to a cross-section of 1.5 &K 107"! 
cm.*/deuteron, 

Figure 3 shows the result for bombardment with 2-m. e. v. electrons, 
plotting per cent activity remaining as a function of dose in electrons/em.?. 
The 87°) dose is 5 & 10" electrons/em.*. The point for the largest dose 
is the result of bombardment with the pulsed beam from the linear accel- 
erator. Between the other two points, obtained with the continuous beam 
from the Van de Graaff machine, the current differed by a factor of five. 
The agreement among these three runs is taken as an indication of rec- 
iprocity in intensity and time for dose as well as that the other experi- 
mental differences involved do not greatly influence the result. 

Before analyzing the data one additional series of runs, shown in figure 

















VoL. 38, 1952. BIOCHEMISTRY: FLUKE, DREW AND POLLARD 185 


4, is presented. These data were taken at liquid air temperature and show 
a very much diminished cross-section for the 3.8-m. e. v. deuterons, about 
one-third that at room temperature. When the ion density was doubled, 
however, resulting in the two open circle points, the cross-section increased 
to larger than that at room temperature. 

In these experiments resuspension efficiency was checked by examining 
the optical densities of the amounts recovered, using a Beckman ultraviolet 
spectrophotometer. A control and a point irradiated with enough deu- 
terons to have inactivated it to about 10°) activity remaining were com- 
pared, 60°) as much material being recovered from the latter as from the 
former. There was no change in the ultra-violet spectrum, which was that 
characteristic of nucleic acid. Our results have been adjusted by cor- 
recting the dose-effect 





curves at 10°) activity 100 ‘ 
remaining and extrapolat- Lew Temperetuve 


ing back to 37°), where 


a 
oOo 


a cross-section was chosen 
between the two curves, 
assuming that the same 


remaining 


~ 
o 


proportion of non-resus- 
pended as of resuspended 
DNA had been inacti- 
vated. The correction 


% activity 
r=) 


r 


at 37% amounts to about 7 











; ; 20x 10° devtAm* 

10% decrease in cross- Sese 

section. A similar cor- a 

rection has been made Bad a ; 

s Per cent activity remaining (log scale) as a function 

for the electron data. of deuteron dose at 90°K. The line represents room 
The primary data pre- temperature deuteron inactivation, and the open circle 

. 5 points are for data at doubled ion density at the low 
sented above, correcting — temperature 


for resuspension, are: 

(1) deuteron cross-section 1.4 X 10~'! em.?/deuteron; (2) electron critical 
dose 5.5 & 10" electrons/cm.?. The electron critical dose is now multt- 
plied by the energy loss per centimeter for 2-m. e. v. electrons in DNA, 
taken to have a composition approximately NisCyCoaH Ps, on the basis of 
the DNA base distribution given by Davidson.'' The density of DNA, 
1.65 g./cm.*, was taken from the same reference. From these figures for 


composition and density the energy loss calculated from the Bethe-Block 


formula!’ is 3.1 m. e. v./em. The only appreciable difference from that 
for protein is due to the increased density. The average energy per pri- 
mary ionization event is taken as 110 e. v./p. i. Multiplying the electron 
critical dose per centimeter squared by 3.1 & 10° e. v./em. and dividing 
by 110 e. v./p. i. gives 1.5 & 10" p. i./em.*. The sensitive volume in 
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which at least one ionization must occur for inactivation is then 6.5 X 
10° em.*/p. i. Multiplying this volume by the density of DNA (which 
cancels its previous use) and by Avogadro's number gives a molecular 
weight: 6,500,000. 

The fact that the DNA deuteron cross-section is larger than that of a 
sphere of the above critical volume is in agreement with much evidence!! 
that the DNA particles are fibrous, or long thin cylinders of high axial 
ratio. Dimensions of a eylindrical shape may be obtained by comparing 
the area and the volume. The most reasonable assumption about the 
orientation of such cylinders in our work is that they are randomly oriented 
in direction, since the dried DNA preparations were about 100,000 A 
thick. For a randomly oriented cylinder of length / and diameter d, the 
average cross-sectional area (assuming the ends negligible) is (m/4)/d. 
Taking (9/4) ld 14 X 10-!! em? and (w7/4)id? = 6.5 X 10738 cm.? 
leads to d 1 A and | 3800 A, giving an axial ratio of about SO. 

Discussion. Our data indicate a sensitive region of molecular weight 
about 6,000,000 associated with the transforming activity of pneumococcus 
DNA. At least one ionization must occur in a region of this size for 
inactivation to result. We may now identify these sensitive regions with 
the actual particles of DNA, since in order for the particles to have a 
smaller size it would be necessary to account for their inactivation by 
radiation-produced events not taking place within the actual structure 
itself and in the absence of effective means of energy transport. 

The molecular weight of nucleic acid has been variously reported," 
ranging from values corresponding to quite simple polynucleotides up to a 
molecular weight of 3,000,000 reported by Gulland, Jordan and Threlfall'* 
for DNA. In determining the molecular weight of nucleic acid a consider- 
able source of uncertainty has been its possible degradation by the methods 


required to separate it from protein. As these methods have been refined 
the trend has been toward higher molecular weights reported for nucleic 
acid. There has also been uncertainty as to whether aggregation may occur 


during purification and extraction procedures. 

By working with a biological activity of nucleic acid rather than with 
the physical properties of entire preparations it is possible to by-pass 
many of these uncertainties. Biologically active DNA 1s presumably 
more closely identified with DNA as it actually exists in place in living 
cells and their fragments. Use of the ionizing particles looks through 
possible degradation products and even heterogeneous materials to show 
the size of these active particles. Aggregation may be ruled out also, 
since activity is not enhanced by the fractionation procedures. 

Our results with low-temperature bombardment need further clari- 
fication. In particular the low-temperature work should be extended 
to include electron bombardment. The decrease in deuteron cross-section 
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at low temperature may miean that the DNA particles are enabled to sus 
tain one ionization within their structures, but no more. Increasing the 
ion density by a factor of two could then give the observed increase in 
cross-section. Another possible extension of the work would be the bom- 
bardment of active DNA before the separation from protein. 

Summary.— Pneumococcus transforming principle DNA has been 
bombarded with 2-m. e. v. electrons and 3.8-m. e. v. deuterons with indica- 
tion by assay of activity remaining of a molecular weight of about 6,000,000. 
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STUDIES ON THK HEMOGLOBIN OF COOLEYV’S ANEMIA 
AND COOLEY’S TRAIT* 


By ALEXANDER Ricut 


GATES AND CRELLIN LABORATORIES OF CHEMISTRY, CALIFORNIA INSTITUTE OF TECH 
NOLOGY, PASADENA, CALIFORNIA 


Communicated by Linus Pauling, January 26, 1952 


The diseases sickle cell anemia and Cooley's anemia (also known as 
thalassemia or Mediterranean anemia) have associated with them ‘“‘minor™’ 
or ‘trait’ forms of the disease in which the symptomatology is minimal, 
and the abnormalities in the red cells are correspondingly less. Extensive 


genetic studies carried out on these diseases! * have led to the conclusion 


that in each case the trait is inherited as a Mendalian dominant, the person 
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being heterozygous in the sickle cell gene or Cooley's gene. Homozygosity 
in either of these genes leads to the corresponding anemia. 

Our understanding of sickle cell anemia has improved considerably with 
the discovery of an abnormal hemoglobin and a molecular interpretation 
of the pathological behavior of the red cell.° In that disease it was pos- 
sible to demonstrate heterozygosity and homozygosity directly, by observ- 
ing varying amounts of the abnormal hemoglobin. 

Liquori® has recently reported the presence of fetal hemoglobin in cases 
of Cooley's anemia. We have carried out a number of experiments which 
support his findings, and have extended the investigation to include cases 
of Cooley's trait, for which no fetal hemoglobin was found. A hypothesis 
is advanced to explain these facts. 

Methods.-Blood from eleven individuals, representing five different 
families, was investigated (table 1). All of the five patients with Cooley's 


TABLE 1 
PATIENT AGE ix RELATION DIAGNOSIS 
Cooley’s anemia 
Cooley’s trait 
Cooley’s trait 
Cooley's anemia 


S. F 18 mos 8 
D.F 3 yrs Brother of S. F. 
M.F. Adult } Mother of S. F. 
V.P syrs 


J.P 
J. 


PM 


M.L 


5yrs 
Adult 


Adult 


19 mos 


Brother of V. P. 


Maternal Aunt 
Vf, 

Maternal Aunt 
VF. 


Father of M. L. 


Cooley’s trait 
Cooley’s trait 


Cooley’s trait 


Cooley's anemia 
Cooley's trait 


ii & Adult 
B.S 20 mos 
W.B 18 mos 


Cooley's anemia 
Cooley’s anemia 


anemia were typical in that they came from families of Italian extraction 
and displayed the characteristic hematological, roentgenographic and 
physical findings.’ 

None of the patients with Cooley's anemia had been transfused within 
three months. Patients B.S. and W. B. had never received transfusions; 
M. L. had been transfused 101 days before our blood specimen was taken, 
and the other two had not been transfused for over four months. 

Blood was collected from these patients using sodium citrate as an anti- 
coagulant, and stroma-free concentrated solutions of hemoglobin were 
prepared by the method of Drabkin.* Blood samples from new born in- 
fants were taken from the placental end of the cord and prepared by the 
same method. The carbon monoxide derivative was used for both the 
electrophoretic and spectrophotometric studies. 

Electrophoresis. Two types of electrophoresis experiments were under- 
taken. Those at 0.1 ionic strength (potassium phosphate buffer) and pH 
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range 5.9 to 7.0 were carried out in an attempt to separate components 
without introducing large boundary anomalies. Another group of ex- 
periments were performed in a 0.03 ionic strength sodium phosphate buffer 


at PH 8.2, in which the ratio of protein concentration to buffer ions is large 
and the resultant boundary anomalies permit the detection of components 
which differ only slightly in mobility." The latter experiments were modi- 
fications of a procedure described fully by Beaven, Hoch and Holiday” in 
an investigation of adult and fetal hemoglobins. 

The carbonmonoxyhemoglobin solutions were diluted with the phos- 
phate buffer preparation to a concentration of 0.5 to 1.0%), and were then 


rABLE 2 
ELECTROPHORETIC MOBILITY AS A FUNCTION OF PH 


PATIENT DIAGNOSIS pH OF BUPPER MOBILITY XK 10° 
cm? voL_t~! sec"! 


S.F Cooley's anemia 5.98 + 1.88 

Cooley’s anemia 9 —().27 

Vis Cooley's anemia >. 50 +0.45 

>. 50 +0.50 

>. 92 —().14 

1.98 —().29 

Cooley’s trait 91 —0.33 

Cooley's trait 5.98 1.81 

9 —0.20 

d.91 —().19 

Cooley's trait 5.93 1.84 

88 —().20 

>.91 —0.19 

Cooley’s trait 59.93 1.78 

4.75 0.04 

>. 88 —0.19 

r90 —0.16 

“ Average of ascending and descending limbs 
Buffer = 0.1 ionic strength phosphate; gradient = 3.8-6.7 volts cm~'; protein con- 
centration = 0.8-1.0g per 100 ml 


dialyzed against the buffer at 4°C. for 12 to 24 hours in a carbon monoxide 
atmosphere. Values of PH were read on samples of the dialyzate after 
they had come to room temperature. 

The modified Tiselius apparatus described by Swingle'! was used. The 
experiments at 0.1 ionic strength employed potential gradients of 3.8 to 6.7 
volts per cm. and those at 0.03 ionic strength employed gradients of 17.5 
to 18.7 volts percm. In the latter experiments, the boundaries were kept 
within the cell compartment by counter compensation. Before each elec- 
trophoresis the diluted carbonmonoxyhemoglobin solution was run in the 
preparative ultracentrifuge at 40,000 g for 40 minutes. It was found that 
this procedure effectively removed small amounts of colorless (presumably 
stromal) contaminants. 
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Alkaline Denaturation. Some experiments were performed to determine 
the resistance of the hemoglobin from Cooley's anemia patients to alkaline 
denaturation, A small volume (0.05 ml) of the blood sample was diluted 
with 3 ml of N/12 NaOH to produce a 0.20% carbonmonoxyhemoglobin 

solution. This was read spectro- 





photometrically at 650 my, as 
suggested by Jonxis,'? and the 
(a) formation of the brown alkaline 
‘aint aed hemochromogen was recorded as a 
function of time. 
Spectrophotometry..-The concen- 
trated = carbonmonoxyhemoglobin 
solutions were diluted with carbon 
monoxide-saturated distilled water 
and were then read immediately in a 
(c) Beckman spectrophotometer, model 
ee NEMBORN DU. The shit width was 0.25 mm, 
which corresponds to a_ spectral 
(9) width at half maximum of 0.88 muy. 
since alacicacaes The dilutions were always made to 
give a standard optical density, 
0.430 at 290 muy. 
le) For the determination of small 
COOLEY'S TRAIT differences in optical density, mul- 


~_ 


(b) 
NEWBORN 











tiple readings were made and the 
standard error was determined. 
Results..-In table 2 are listed 
the mobilities of some of the hemo- 
globin preparations at different PH 


FIGURE 1 
Electrophoretic patterns of the ascend 
ing boundary of carbonmonoxyhemo 


globins. (a) Normal adult hemoglobin 
(b) Newborn cord blood. (c) 50:50 mix- values. Special emphasis was placed 
ture of (a) and (b). (d) Cooley's anemia — upon examination of the hemoglobin 
nt ee a M “ % peed of individuals with Cooley's trait, 
trait hemoglobin, patient .. Sodium : : m if ; 
. ' mid wet Ae since it was felt that these might 
phosphate buffer at pH 8.2; ionic strength . . 
0.03; carbonmonoxyhemoglobin concen represent cases with mixtures of 
tration 1.08 g per 100 ml.; time 280 to 300 different types of molecules, in 
min. The boundary is moving to the right analogy with sickle cell trait. 
However, only one peak was con- 
sistently evident. The mobility of the trait blood between the pH 
values 6.88 and 6.92 was —0.21 + 0.05 em®* volt! sec ~', and the average 


value for the two anemia samples also was —0.21. A comparison of these 
values with the mobility-PH results cited for normal carbonmonoxyhemo- 
globin® shows that they are equivalent within experimental error. 

Because of Liquori’s report,® mobilities were measured of mixtures of 
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hemoglobin from Cooley's anemia patients with hemoglobin from normal 
adults and newborns. In no case could two peaks be resolved, and the mo- 
bility results were similar to those cited above. Fetal hemoglobin was 
found to migrate effectively with the same mobility as the adult form in the 
pH range at 0.1 ionic strength. 

Additional experiments were then carried out at the higher PH values 
and lower ionic strengths in which it has been possible to resolve fetal and 
adult hemoglobins electrophoretically."° 

In figure | are reproduced some of the patterns obtained from various 
blood samples. It can be seen that normal adult carbonmonoxyhemoglobin 
(a) produces a nearly symmetric, rapidly moving peak, in contrast to cord 
blood (6), for which two peaks are Gbserved. The pattern for a mixture 
of the two (c) permits identification of the faster component as adult hemo 
globin and of the slower as fetal hemoglobin. 


TABLE 3 
ELECTROPHORETIC ANALYSES OF COOLEY'S ANEMIA AND FETAL CARBONMONOXYHEMO- 
GLOBINS 


APPARENT ‘7, OF CORRECTED ‘, OF 
BLOOD FASTER COMPONENT FASTER COMPONENT 


B. S.--Cooley’s anemia 6.5 2.8 
W. B.-Cooley’s anemia 5.6 2.0 
M.L.—-Cooley's anemia 63.4 
Newborn 11.0 
Newborn and normal 

adult (50:50 mixture ) 83. 


62. 


Ascending boundary only; pH 8.2, ionic strength 0.1, concentration 1.08 g per 100 
ml, time 280-300 min. 


Hemoglobin from Cooley's anemia patients (d) produces a two-peak 
pattern, with mobilities equal to those for a newborn and adult mixture, 
while the hemoglobin from Cooley’s trait (e) produces one peak, indistin- 
guishable from the normal adult (a). 

This analysis has been carried out on the hemoglobin from three cases 


of Cooley’s anemia (table 3) and three cases of Cooley's trait. Each of 


the anemia patterns shows the two components described above, while the 
trait patterns show only one. These results for Cooley's anemia substan 
tiate the observations of Liquori. 

In a buffer medium of low ionic strength in which there are two compo- 
nents with only slightly differing mobility, the faster component tends to 
concentrate itself ahead of the boundary of the slower component to bal- 
ance the decrease in migration velocity. By using Svensson’s equation'* 
describing the conductivity change across a boundary, Hoch-Ligeti and 
Hoch'* * have derived an expression relating the areas in the electrophoretic 
pattern of a two-component system to the concentrations. They assume 
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absence of interaction between the two components and linear dependence 


of migration velocity on concentration. These approximations have been 
justified for the hemoglobin system in 0.03 ionic strength.’ We have uti- 
lized their equation to calculate the amounts of fetal and adult hemoglobins 
in the cases of Cooley's anemia. The apparent and calculated percentages 
of the faster moving (presumably normal adult) component are given in 


table 3. 

Patients B.S. and W. B. have a very small amount of the faster moving 
material. It is interesting to note that the specimens B. S. and W. B. 
represent unusually pure samples of fetal hemoglobin. Newborn infants 
have 15 to 20% of adult hemoglobin, 
as observed — electrophoretically, 
while even the 20-week fetus is 
found to have 6%." This nearly 
pure preparation of fetal hemoglobin 
can be used for the calibration of 
spectrophotometric readings. 

Alkaline Denaturation..-We have 
carried out experiments on alkaline 
denaturation, with results in agree- 
ment with those of Liquor. 
Cooley’s anemia hemoglobin was 
Pe found to denature more slowly than 

FIGURE 2 normal adult hemoglobin; samples 

Ultra-violet absorption spectrum of car W. B.and B.S., with only 2 to 3° /, 

bon: monoxyhemoglobins. (a) Adult; (6) of the faster component, denature 


50:50 adult newborn mixtures, (c) new more slowly than samples of hemo 
born. The variation in the difference in 
optical density between 288 and 289 my 


4 


globin from newborns, in comforin- 
ity with the electrophoretic analyses. 
all go through the same point at 290 my Spectrophotometry..Jope'® dis- 
but are displaced vertically here for clarity. covered a difference in the ultra- 
violet absorption spectra of fetal 
and adult hemoglobins. The inflection point of the tryptophan fine 
structure band is found at 2910 A in fetal and at 2898 A in normal 
adult hemoglobin. Beaven, Holiday and Jope'® have measured the 
position of the inflection point in a series of fetal and infant bloods 
and have correlated it with the ratio of fetal and adult hemoglobin. 
Beaven, Hoch and Holiday" have shown that the position of the inflection 
point for artificial fetal-adult mixtures is not a linear function of concentra 
tion, since a small amount of fetal hemoglobin has a marked effect in pro 
ducing a shift in the inflection point. Thus, they find their spectral method 
only of limited usefulness. 
In the course of checking their work, we have devised a method which 


can be readily seen. The curves actually 
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we believe is of somewhat greater utility both in instrumentation and range 
of applicability. Curves such as those illustrated in figure 2 were obtained 
on the Beckman Model DU spectrophotometer. Curve (a) is the absorp- 
tion spectrum of normal adult carbonmonoxyhemoglobin when examined 
with a slit width of O.SS my, curve (c) is the spectrum for a sample of new- 
born carbonmonoxyhemoglobin, and curve (4) is the spectrum of a 50:50 
mixture of the adult and newborn hemoglobins, the total concentration 
being constant. The increased prominence of the tryptophan shoulder can 
be seen in curve (c), although there is no resolution of a separate peak. 
Instead of attempting to determine the position of the inflection point, we 
measure the change in optical density for the interval 288 to 289 mu. 

If Beer's law is valid, the difference in optical density over this interval 
should be a linear function of the concentration of fetal or adult hemoglobin. 


rABLE 4 
SPECTROPHOTOMETRIC AND ELECTROPHORETIC ANALYSES OF CooLEYy’S ANEMIA AND 
Coo.ey’s Trarr CARBONMONOX YHEMOGLOBIN 


A288- 280 “¢ ADULT HEMOGLOBIN 
BLOOD x10 SPECTROPHOTOMETRIC ELECTROPHORETIC 


B.S.—-Cooley’s 
anemia 
W. B.—Cooley’s 
anemia 
M. L.-Cooley’s 
anemia 
L.— Cooley's 
trait 12.5 % LOO 
5. J. C.—-Cooley’s 
trait 13.0 101 100 


To test the assumption of Beer's law, artificial mixtures of adult and 
newborn hemoglobins were prepared, and the differences in optical density 
between 288 and 289 my (A288-289) were measured. All of these mix- 
tures were prepared to a constant total concentration initially, and the per- 
centage of adult hemoglobin in the newborn sample used was determined 


electrophoretically as cited above. The results, plotted in figure 3, show 
a linear dependence of A2SS-289 upon the concentration of adult hemo- 


globin. 

Determination of small differences in optical density is possible provided 
one exercises great care in the measurements. By taking multiple measure- 
ments in the region near 0.400, optical densities can be determined to 
+().0005. This corresponds to an error of +6°) in the determination of 
the amount of the adult component. 

The accuracy of the curve in figure 3 was tested by determining A288—-289 
for various samples of hemoglobin, the fetal-adult compositions of which 
had been determined electrophoretically. The percentages of adult hemo- 
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globin as obtained by the two methods are given in table 4. They are 
found to be in agreement to within 7°. 

Discussion. Hoch" has described a faster moving component present 
in normal adult hemoglobin in concentrations of 2 to 3%. Our cases 
B.S. and W. B. have 2 to 3% of a component which migrates faster than 
fetal hemoglobin and could be either normal adult hemoglobin or the small 
component described by Hoch. The mobility data are not sufficiently ac- 
curate to settle this question. 

As B.S. and W. B. had never received transfusions, they have continued 
to live beyond the fetal stage with use of fetal hemoglobin alone as oxygen 














1 4 
40 100% ADULT 
60 O % FETAL 


FIGURE 3 


The change in the difference in optical density between 288 and 289 
my as a function of adult hemoglobin concentration in fetal-adult 
mixtures 


carrier. To our knowledge, these are the first such cases reported in the 
literature. 

Because M. L. had received a transfusion 101 days before the present 
sample was taken, we do not know with certainty what percentage of her 
normal adult hemoglobin was derived endogenously. Normal adult eryth- 
rocytes continue to live out their normal life span when transfused into 
patients with Cooley’s anemia.” From a record of M. L.’s transfusion 
and hemoglobin levels, we estimate that one-half or less of her adult hemo- 
globin was obtained by transfusion. Liquori estimates that his two eight- 
year-old patients have 40 to 60°) fetal hemoglobin. 
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We must explain why the homozygous patient with Cooley's anemia has 
large amounts of what we believe to be fetal hemoglobin, nearly 100°;, 
while the heterozygous parent or sibling has no fetal hemoglobin. 

Sickle cell disease differs from Cooley's anemia in that the presence of 
one mutant allele in the heterozygote results in the production of an al- 
tered hemoglobin molecule. Further, the homozygous patient with Coo- 
ley's anemia produces a molecule which appears not to be unique to the dis- 
ease, but to be identical with a hemoglobin molecule normally produced 
at an earlier stage of development. The disease has prevented the cessa- 
tion or has reactivated a metabolic process which normally stops at birth. 

In order to bring this about, the Cooley's gene has both effectively to 
block the production of normal hemoglobin and to cause the continued 
production of fetal hemoglobin. We suggest that the primary effect of 
Cooley's gene is in some way to prevent the production of normal erythro- 
cytes, and that there then occurs the continued production of erythrocytes 
containing fetal hemoglobin as a result of the consequent anemia. Accord- 
ing to this interpretation fetal hemoglobin is not similar to sickle cell hemo- 
globin in being the aberrant product of an altered gene, but is instead the 
result of a compensatory response. 

If this interpretation is valid, there should be other cases of severe 
anemia with which is associated the presence of fetal hemoglobin beyond 
the fetal stage. That this may be the case is suggested by the studies of 
Singer, Chernoff and Singer'* in which they interpret a prolonged alkaline 
denaturation time as possibly indicating small amounts of fetal hemoglobin 
in various anemias. 

The patients with Cooley's trait have one noimal allele and one Cooley's 
allele, which together have the capacity to produce malformed, but stable 
erythrocytes. These patients have no severe anemia, and hence no strong 
hematopoetic stimulus, and it is observed that they do not produce cells 
containing fetal hemoglobin. 

The exact nature of the defect in Cooley's anemia remains uncertain. 
The mutant gene may prevent some necessary step in the production of 
adult hemoglobin, or it may prevent the synthesis of normal adult stromal 
material, or in some other way impair normal red cell production. Any 
one of these defects would produce a severe anemia which, according to our 


postulate, would activate the production of cells containing fetal hemo 


globin. 
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THREONINE INHIBITION IN A STRAIN OF NEUROSPORA 
By C. O. DouDNEY AND R. P. WAGNER 
GENETICS LABORATORY OF THE DEPARTMENT OF ZOOLOGY, UNIVERSITY OF TEXAS 
Communicated by J. T. Patterson, January 4, 1952 


It is becoming increasingly apparent that metabolites normally present 
within the cell may act as inhibitors of some phases of metabolism while 
at the same time being necessary as intermediates in others.! This fact 
is of interest not only in the analysis of growth factor requirements, but 
also in the analysis of genetic blocks which prevent the production of cell 
metabolites. For it would follow that the inherited inability of an organ- 


ism to synthesize a compound may not necessarily be due at all times to 
the absence of an enzyme catalyzing its synthesis, but rather to an inhibi- 
tion of its synthesis by some other compound produced by the organism 


or supplied from the outside. 

Work previously done in this laboratory on the pantothenicless mutants 
of Neurospora crassa has indicated that some such explanation as given 
above may apply to certain biochemical mutants of Neurospora, since the 
enzymes involved in the synthesis of pantothenic acid are present both in 
the pantothenicless mutants and the wild type strain.*~* Although the 
nature of the mechanism causing the enzyme system to be inactive in the 
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pantothenicless mutants is unknown, the presence of some type of inhibi- 
tion is definitely indicated. 

As a preliminary approach to the investigation of inhibitors as one 
of the factors involved in the formation of genetic blocks, a search was 
made for mutants of Neurospora which would grow like wild type on 
minimal medium, but be inhibited unlike wild type by specific compounds 
known to participate in the metabolism of Neurospora. One of the first 
of such mutants to be found, a strain which is strongly inhibited by threo- 
nine, is described in this communication together with experiments demon- 
strating the relief of threonine inhibition by other compounds. 
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FIGURE 1 


The effect of threonine on the growth of UT 77 and Wild Type, 


Emerson 5256 


E:xperimental,The threonine inhibited strain, UT 77 originated from 
a single ascospore derived from a cross involving conidia irradiated with 
ultra-violet light. The method of Lein, Mitchell and Houlahan® for se- 
lective screening of biochemical mutants was used with one modification 
to select for the mutant ascospores. A high concentration of threonine 
and other amino acids was added to the agar plates upon which the spores 
were collected for germination. The threonine inhibited strain was 
separated from the other biochemical mutants detected by further screen- 
ing in minimal agar tubes 

Compared to the wild type Emerson 5256 from which it was derived, 
the mutant UT 77 grows at about the same rate on minimal medium at 
35°C., but about 50% slower than wild type at 25°C. Threonine inhibition 
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is most strongly expressed at 45°C. and at pH values below 5.8. At pH 
values above 5.8 threonine inhibition is slight and of short duration. 
Accordingly, in the experiments described below the temperature was 
maintained at 35°C. and the initial pH was adjusted to 4.8. Experiments 
were carried out using 125-ml. flasks to which 25 ml. of medium was added. 
All experimental data given are based on the mycelial dry weights for three 
flasks with identical medium. The minimal medium used was that 
described by Beadle and Tatum.’ 

Figure | compares the response of UT 77 and wild type after 60 hours 
growth on minimal medium supplemented with various concentrations 
of L- and DL-threonine. Maximum inhibition by L-threonine is given 
at a concentration of | * 10°* molar, whereas 2.3 * 10°* molar DL- 
threonine is required to give maximum inhibition, At these concentrations 
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Phe inhibition of UT 77 by choline, and the relief of threo 
nine inhibition by choline 


of threonime growth of UT 77 is only 10 to 20°; of that produced on mini 
mal medium. Threonine is, on the other hand, slightly stimulatory to the 
wild type strain. 

Threonine inhibition of UT 77 is effective for a period of about SO hours 
after which the mutant grows at a slowed rate compared to wild type. 
A final level of growth is attained in 152 hours (99 mg. dry weight of 
mycelium in 25 ml, of medium) which is equal to the final level of growth 
attained in SS hours by the mutant on a medium free of added threonine. 
The final level of growth of UT 77 on minimal surpasses that of the wild 
type by 35° . 
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In addition to threonine, choline also inhibits the growth of UT 77. 
Choline at concentrations of 7.5 & 10-4 molar and higher inhibits growth 
to 65° of that obtained on minimal medium in 60 hours (Fig. 2). Choline 
does not inhibit the wild type strain. Complete relief of choline inhibition 
can be obtained by the addition of as little as 0.15 wg. of thiamin per ml, 


of medium, (table 1). No other B-vitamin tested shows this relief. 


TABLE | 
THIAMIN RELIEF OF CHOLINE INHIBITION 
DRY WRIGHT MYCELIUM, MG 


CHOLINE UT 77 WILD TYPE 


CONCENTRATION NO THLAMIN THIAMIN” NO THIAMIN THILAMIN® 


0 
56 X 10° M 
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28 


+} 
th 
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17 
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OX 10-3 M PA 1 
5 X 10-3? M 28 Wy 52 50 


* Thiamin concentration 0.15 pg. per ml. 
Choline partially antagonizes the inhibition produced by threonine 


(Fig. 2). The maximum relief of threonine inhibition obtained gives 
about the same level of growth as obtained with choline alone. 
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The relief of threonine inhibition of UT 77 by methionine and 
homocysteine-thiolactone 


Four compounds reverse threonine inhibition completely: L- or DL 
methionine, DL-homocysteine, DL-homocysteme-thiolactone and DL 
homoserine (Fig. 3 and table 2). DL-methionine is as effective as L 
methionine and in both cases the reversal of threonine inhibition is non 
competitive, siiice no increase of methionine concentration is required to 
effect the seme level of relief of inhibition when the threonine level is 
increased. DL-homocysteine and DL-homoserine are appreciably less 
effective than methionine for the reversal of threonine inhibition (table 2). 
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Relief of threonine inhibition by homocysteine-thiolactone is clearly of a 
competitive type, since an equivalent increase in the concentration of 
the thiolactone is required to antagonize an increased amount of threonine 
in the medium (Fig. 3). Homoserine, while slightly more effective than 
homocysteine in the competitive relief of threonine inhibition in the lower 
concentration range, requires a many fold increase in concentration to 
relieve threonine inhibition at higher concentrations of threonine (table 2). 
It should be noted in this respect that homocysteine is a highly unstable 
compound in the culture medium, being subject to oxidation to homo- 
cysteme, which is inactive in the support of growth of Neurospora.” 
rABLE 2 . 
MINIMAL MoLarR CONCENTRATIONS OF VARIOUS COMPOUNDS REQUIRED TO COMPLE 
REVERSE L-THREONINE INHIBITION 


MOLAR CONCENTRATION L-THREONINE 
10 xk 10 2.0 * 10 


L-methionine 5.0 X 107% 5.0 X 107 


oO 

DL-homocysteine 2.4 X 10 $9 10 

DL-homoserine 2X 107 7.5 X 10 

DL-homocysteine-thiolactone 5.0 * 10 10x 10 
e 


Sulfanilamide 5 X 107° 5 x 107° 


Sulfanilamide relieves the inhibition of threonine completely at 7.5 X 


10°-® molar (Fig. 4). This relief is non-competitive (table 2). If the 


@ L- THREONINE 
4 NO THREONINE 


an 
°o 


nN 
o 


—-e— 


> 
= 


» 
i.e) 


e LEVEL OF THREONINE INHIBITION 


MGMS. ORY WEIGHT OF MYCELIUM 
$ 


=) 


4 A 4 A 1 i L 
75x107® i107? 25 5 75 ixl07*) 25 
MOLAR CONCENTRATION SULFONAMIDE 
FIGURE 4 
The relief of threonine inhibition of UT 77 by the sulfonamide, sulfanil 
amide 


sulfanilamide concentration is increased above 7.5 & 10~° molar, growth 
of both UT 77 and the wild type is suppressed. As expected, p-amino- 
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benzoic acid antagonizes sulfanilamide reversal of threonine inhibition 
competitively; one mole of p-aminobenzoic acid antagonizes 100 moles 
of sulfanilamide. In the absence of sulfanilamide, p-aminobenzoic acid 
has no effect on the growth of the mutant even at concentrations as high 
as | X 10~° molar. 

Genetic studies of a cross of UT 77 to wild type show that threonine 
inhibition is due to a single gene mutation The difference between the 
inhibited and non-inhibited segregants obtained from the cross to wild 
type is decisive. For example, of 31 segregants derived from ascospores 


picked at random, 17 showed definite inhibition by threonine ranging from 
55 to 100°, with a mean of 77 + 120). On the other hand, the remaining 
14 segregants were inhibited from 0 to 10°; with a mean of 3.6 + 2.30; 
inhibition. It may be assumed therefore that the threonine inhibition 
character is a distinct one dependent upon a single gene change. 


Discussion. The results of the growth response experiments demonstrate 
that UT 77 grows as well as wild type controls at 85°C. on minimal medium 
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The relationship postulated for threonine, homocysteine and threonine metabolism 
in Neurospora. 


in flask cultures, but is markedly inhibited by the addition of the amino 
acid, threonine to the medium. DJL-threonine is about one-half as active 
as L-threonine, suggesting that the D-isomer is inactive in the process 
and that the observed suppression of growth is specifically produced by 
L-threonine. This 1s in agreement with the finding that the threonineless 
mutants of Neurospora cannot use )-threonine for growth.’ 

Threonine inhibition can be reversed by methionine non-competitively, 
and by homocysteine and homoserine in a competitive manner at higher 
concentrations. Horowitz* has shown that homocysteine is a precursor 
of methionine in Neurospora. (See Fig. 5.) The results reported here 
indicate that threonine causes inhibition by reacting with homocysteine 
(or a product to which it is readily converted) in an essentially irreversible 
reaction. The postulated reaction may be assumed to prevent growth 
by engaging the homocysteine synthesized by the mutant strain, thus 
preventing the formation of methionine. Therefore the strain becomes 
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“methionineless’ or almost so when threonine is in the medium. The 
addition of methionine is able to relieve threonine inhibition by satisfying 
the threonine created methionine deficiency. Since homocysteine relieves 
the inhibition of threonine competitively, it, or a product derived from it, 
would seem to be involved both in reacting with threonine and in meth- 
ionine production. When sufficient homocysteine is available in’ the 
medium to saturate the reaction with threonine, the synthesis of methionine 
can proceed at an adequate rate allowing usual growth. 

Homoserine is believed to be a common precursor of both threonine and 
methionine, the latter being formed through cystathionine and homo- 
cysteine as intermediates.’ Recently homoserine was demonstrated for 
the first time in a biological system when Fling and Horowitz’ showed it 
to be present in extracts of a methionineless mutant of Neurospora. 
Table 2 shows that in UT 77 homoserine surpasses the effectiveness of 
homocysteine in the competitive relief of inhibition by low concentrations 
of threonine, while it becomes much less effective at higher concentrations 
of threonine in producing complete relief. Presumably, therefore homo- 
serine exerts its effect by becoming converted to homocysteine at a rate 
sufficient to antagonize lower concentrations of threonine, but inadequate 
to effect complete relief of threonine inhibition when threonine is in the 
medium at higher concentrations. 

The methyl group of methionine is thought to come from choline in 
Neurospora." '® Mutant strains which cannot make their own choline 
can use methionine as a partial replacement for choline. The inhibition 
observed in the case of choline can be best explained by assuming that 
choline accelerates the production of methionine by donating methyl 
groups to homocysteine, thereby depriving the threonine-homocysteine 
reaction of homocysteine, and reducing the synthesis of an essential 
metabolite. This concept of the action of choline in the system is sup- 
ported by the fact that choline partially antagonizes the inhibition by 
threonine indicating that the production of methionine from homocysteine 
is aided by the presence of appreciable amounts of choline in the medium. 

It can be concluded that p-aminobenzoie acid has a role in the postulated 
threonine-homocysteine reaction system because of the sulfanilamide 
relief of threonine inhibition and its reversal by p-aminobenzoic acid. It 
must be assumed that the amount of sulfanilamide required to reverse 
threonine inhibition completely is sufficient to antagonize the p-amino- 
benzoic acid produced by the mutant strain to a level limiting to the 
postulated threonine-homocysteine reaction, and thereby making adequate 
amounts of homocysteine available for the production of methionine. 

According to the hypothesis presented above choline acts to deprive 


the postulated reaction between homocysteine and threonine of homo- 


cysteine. The fact that thiamin in low concentrations will reverse choline 
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inhibition suggests that the reaction between homocysteine and threonine 
may be involved in the biosynthesis of thiamin. Thus thiamin added to 
the medium would supply the deficit produced by the action of choline 
and allow normal growth. Further experiments performed in this labora- 
tory support the view that threonine is a possible precursor to the thiazole 
portion of thiamin. These experiments will be described in another 
publication. 

The results reported here with UT 77 can best be explained by assuming 
that when either of the precursors threonine or choline of the competing 
reactions is available in excess amount, that reaction is caused to go at an 
excessive rate, thereby producing a deficiency of the common precursor, 
homocysteine. Since inhibition is not produced in the wild type strain 
by either choline or threonine, a mechanism must exist in wild type which 
exerts a controlling influence over the rates of these competing reactions, 
preventing either from going to excess and depriving the other of the 
common precursor. These hypotheses together with those based on 
previous work with methionineless and threonineless mutants of Neuro- 
spora’~*: '® 16 which establish the relationship of threonine and methionine 
in biosynthesis are summarized in fig. 5. 

Emerson and Cushing'® and Emerson” report a mutation, sfo in Neuro- 
spora which causes a requirement for sulfonamide at 35°C. and at pH 
values below 6. It was postulated by Zalokar" that this strain is in- 
hibited by the p-aminobenzoie acid normally produced, and that sulfon- 
amide promotes growth by antagonizing this inhibition. Zalokar,'’ using 
the sfo character in combination with a genetic block preventing the 
synthesis of methionine, showed that low concentrations of methionine 
support growth in the absence of sulfonamide while higher concentrations 
required that sulfonamide be present in order for growth to take place. 
A similar effect was found in the case of homocysteine. It was further 
found that threonine as the sole supplement to minimal medium will 
permit the sfo strain to grow. These results have been interpreted as 
meaning that the mutation responsible for the sulfonamide requirement 
controls a deleterious reaction which, using a precursor of methionine as a 
substrate, is catalyzed by p-aminobenzoic acid and interferes with the 
normal utilization of threonine. 

It seems probable that the same reaction is responsible for the unusual 
growth responses described for the sfo mutant and UT 77. In each mutant 
the sulfonamide to p-aminobenzoic acid ratio demonstrated is approxi- 
mately 100,"° and in addition the mutant phenotype is best expressed at 
35°C. and at the lower pH ranges. Furthermore, the interacting sub- 


strates are the same for each mutant, homocysteine or its immediate 
= o 


product, and threonine. 
The differences observed between the sfo mutant and UT 77 must be 
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due to a specific genetic difference which allows the same reaction to exert 
a detrimental effect on growth, but in the case of the sfo strain interferes 
with an aspect of threonine utilization, while in UT 77 it prevents the 
formation of methionine when excessive amounts of threonine are available. 
It is to be noted that the sfo strain produces enough homocysteine to inhibit 
itself in minimal medium, but that UT 77 will grow about as well as the 
wild type strain in minimal at 35°C. Presumably, therefore, there is a 
balance in the production of the precursors to the postulated reaction in 


this case, 

From a genetic viewpoint, the behavior of UT 77 supports the concept 
that the inability to mediate a chemical synthesis by an organism due to 
a genetic block is not necessarily due to the absence or inactivation of an 
enzyme required for that synthesis, but may be due to a metabolite in- 


hibition resulting from the mutation. 

The fact that inhibited strains such as UT 77 may be induced by the 
same means as are used to induce the ordinary biochemical mutants with 
nutritional deficiencies, makes it evident that the phenomenon of inhibition 
must be regarded as an important possible causative factor of genetic 
blocks. Although UT 77 is not a biochemical mutant on minimal medium, 
it is on exogenous threonine, and it is not too far fetched to assume that a 
more drastic mutation may cause it to be inhibited by its own threonine, 
and thus become a biochemical mutant on minimal medium requiring 
methionine or its precursors. ‘The present data by no means prove that 
genetic blocks are produced in this way in biochemical mutants, but they 
do definitely indicate that further work predicated on this assumption 
is warranted. 

Summary. A new mutant strain of Neurospora is described which 
gives wild type growth on minimal medium at 35°C., but which is strongly 
inhibited by L-threonine. This inhibition is completely overcome by 
methionine, homocysteine, homoserine and sulfanilamide, and partially 
overcome by choline. Choline also inhibits the mutant, but to a lesser 
degree than threonine. The choline inhibition is relieved by thiamin. 
The significance of these observations is discussed, and their possible 
bearing on the understanding of the nature of genetic blocks considered. 
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OBSERVATIONS ON THE BEHAVIOR OF SUPPRESSORS IN 
NEUROSPORA * 
By Mary B. MITCHELL AND HERSCHEL K. MITCHELL 


KERCKHOFF LABORATORIES OF BroLoGy, CALIFORNIA INSTITUTE OF TECHNOLOGY, 
PASADENA, CALIFORNIA 


Communicated by G. W. Beadle, January 14, 1952 


A suppressor of pyrimidineless 3a (37301) and some aspects of the be 


havior of the suppressed mutant have been described earlier.! The obser- 
vation that lysine, ornithine, citrulline and arginine influence growth re 
sponses of the suppressed mutant suggested studies of the behavior of re 


combinants involving pyr Ja and s and mutants having requirements for 
these amino acids. — Effects of the pyrimidineless mutant and its suppressor 
upon certain lysine-requiring mutants have been reported.” The present 
paper deals with a somewhat greater variety of interactions observed be 
tween pyr Sa and s and mutants which utilize proline, ornithine, citrulline 
or arginine.’ These interactions include suppression of two non-alleliec 
prolineless mutants by the pyrimidineless suppressor and partial sup 
pression of pyr 3a by three non-allelic ornithineless mutants. 

I:ffects of the pyrimidineless suppressor on mutants of the proline-arginine 
series. The suppressed pyrimidineless mutant, pyr Ja-s (37301-s) was 
first described! as being characterized on minimal medium by a lower 
growth rate than that of wild type, and by stimulation by pyrimidine, ly 
sine or histidine. It has since been found that this is true only of some 
isolates. Others grow as rapidly as wild type, still others very nearly so, 
and the stimulations are very slight. ‘The reported inhibition by ornithine, 
citrulline and arginine and its release by lysine remain properties of the re 
cent isolates, however, although there may be differences in the effective- 
ness of the compounds as inhibitors. Proline is also inhibitory but very 
slightly so. 

Crosses of pyr 3a-s were made to the following strains: 21863, prolineless; 
S401, 44207 and 51506, which use proline, ornithine, citrulline or arginine ; 
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27947, 29997 and 34105, which utilize ornithine, citrulline or arginine; 
30300 and 33442, which utilize citrulline or arginine. For convenience 
these mutants have been renamed in the following fashion: 
prolineless 
21863 = proll = $5401 = prol2 = 44207 = prol3 = 51506 = prol 4 
ornithineless 


27947 = orn 1 29997 = orn ? 34105 = orn 3 


citrullineless 
30300 = cit 1 33442 = cit 2 


Growth responses and genetic relationships of these strains have been de- 
scribed by Srb and Horowitz,* Srb‘ and Srb, Fincham and Bonner,® who 
have shown that, except for prol 3 and 4, which behave as alleles, strains 
among these which use the same compounds are genetically distinct. 

Ascus dissections were done on minimal agar plates as described by Haas 
et al.,® and the spores were heat treated and allowed to germinate on the 
plates. The germinated spores were transferred to tubes containing min- 
imal agar medium supplemented with cytidine and proline or arginine and 
the resulting cultures classified by means of tests in tubes of liquid minimal 
medium with appropriate supplements. Members of spore pairs were 
isolated and cultured together as reported by Fincham? and Haskins 
and Mitchell.* 

In tables | and 2 the types of asci observed are tabulated so as to show 
the phenotypes, but not the order in the ascus, of the spore pairs. The 
ornithineless and citrullineless mutants are designated as arg in order to sim- 
plify the tabulation in table 2... The suppressor, s, is represented only in 
those spore pairs in which its presence could be detected phenotypically by 
the qualitative classification tests used. From table | it may be seen that 
from the crosses pyr 3a-s & prol 2, 3 and 4 the prolineless mutant appeared 


in only ene of the four spore pairs in asei of types 2, 5, 7 and 12, and in 
types 3, 6 and 10 it did not appear in any spore pair. Hence, it seems that 
prol 2, 3 and 4 are suppressed, and the types of segregation observed are 
consistent with the view that they are suppressed by s. This was con- 
firmed by recovering the three prolineless mutants from out-crosses of 
phenotypically pyr-s strains derived from asci of types 2 and 3, and also by 


recovering them from out-crosses of phenotypically wild strains from asci 
of types 10 and 12. For the sake of clarity it should be pointed out that 
only one of the pyr-s spore pairs from a type 2 ascus, and only one + spore 
pair from a type 12 ascus carried the prolineless mutant, since it appeared, 
unsuppressed, in one spore pair in asci of these types. 

There is no indication that prol 1] is suppressed since 1t appeared in two 
spore pairs from each ascus, although in types 4, 9 and 11 it presumably 
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combined with s in at least one spore pair, as evidenced by the appearance 
of unsuppressed pyr. As table 2 shows, the same result was obtained with 
the ornithineless and citrullineless mutants since they appear to have com- 
bined with s in asci of types 2,4, 5 and 6. The presence of s in an unsup- 
pressed arg strain from each cross was demonstrated by crossing these to pyr 
3a and recovering pyr 3a-s. 


rABLE 1 

CROSSES OF pyr 3a-s TO prol 1, 2,3 AND 4 

TYPES OF SEGREGATION prol 1 prol 2 pro 3 
pyr- pyr-s prol prol | 3 
pyr pyr-s prol + 2 
pyr pyr-s + + 1 
pyr pyr prol prol 
pyr pyr prol + 
pyr pyr + 
pyr pyr-prol + 
pyr pyr-prol prol 
pyr pyr prol 
pyr pyr { 
pyr pyr-prol prol 
pyr pyr-prol + 
pyr-prol pyr-prol + 


23 


NoTre: Only the phenotypes of spore pairs are shown, not the complete genotypes nor 


the order in the ascus 


TABLE 2 
CROSSES OF pyr 3a-s TO orn 1, 2 AND 3 AND cit 1 AND 2 
TYPES OF SEGREGATION orn | orn 2 orn 3 
pyr pyr-s arg 7 
pyr pyr arg l 
pyr pyr-arg arg l 
pyr pyr arg 
pyr arg 
pyr pyr-arg arg 


pyr-arg pyr-arg T 


Notre: Only the phenotypes of spore pairs are given, not the complete genotypes nor 
the order in the ascus. 


Growth responses of recombinants from the crosses were further tested 


by measuring dry weights of mycelium from cultures in 125-ml. flasks 


containing 20 ml. of Fries medium with appropriate supplements. The 
suppressed prolineless mutants prol 2-s, prol 3-s and prol 4-s are not dis- 
tinguishable from wild type on minimal medium and are not affected by 
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added proline or arginine. The suppressed double mutants prol 2-pyr 3a-s 
and prol 3-pyr Ja-s were not found to differ from pyr 3a-s. 

Tests of recombinants of s with the ornithineless and citrullineless mu 
tants have shown that s affects these in a way which may be considered the 


opposite of suppression. Three of these mutants, orn 3, cit 1 and cit 2, 


grow slowly on unsupplemented medium, and in about 10 days reach a dry 
weight which approaches that of wild type after a 4-day growth period. 
As would be expected, they have lower requirements than do orn I and 2, 
which do not grow on unsupplemented medium. All of the five combi 
nations with s, however, fail to grow a weighable amount on unsupplemented 


FABLE 3 


REPRESENTATIVE DRY WeiGuTs (IN MG.) PROM 125-ML. FLASK CULTURES AT 25°C 
AFTER 4 OR LO Days’ GROWTH WITH THE FOLLOWING SUPPLEMENTS IN 20) ML, OF 
Mepium: ORNITHINE = DL Orniruine HCI, 1 MoG.; CirruLtting = DL CrrrRuine, 
1 Mo.; ARGININE = L ARGININE HCI, 1 Mo.; Casein = ENzyME-HypbrRouizep, 40 MG 


ORKNITHINE CITRULLINE ARGININE CASEIN 
DAYS i 10 4 10 4 10 4 
orn 1 ‘ 21.5 30.5 nay 57 


orn 2? { 27 29.5 

orn 3 4 oe ath) 51 61 
cut 1 7 2. 5 3 16 

cu 2 §.! j 57 42.5 

orn 1-s 13.6 2) 


» 


orn 2-s 14 


> 


orn 3-s 

cu 1-s 15 ‘ 

cu 2-s 

pyr 3 

pyr Sa-orn 1 

pyr 3a-orn 2 

pyr Ba-orn Go ‘ 9.¢ ‘ o3f tr 

pyr Sa-cit 1 0 0 0 0 
Nore: tr. (trace) represents a quantity of mycelium which ts visible but too small to 


be removed and weighed 


medium in 10 days, and give about the same response to citrulline or argi 
nine as orn J and 2. Furthermore, the recombinants orn 1-s, orn 2-s and 
orn 3-s all fail to respond to ornithine even at a level which is 4 to 4 times 
higher than that required for half-maximum growth of orn I and 2. Hence 
the presence of s appears to convert these five mutants into strains with 
absolute requirements for citrulline or arginine. ‘This is shown by data 
given in tables 8 and 4. 

The effect of ornithineless mutants on the pyrimidine requirements of 
pyr 3a. In classification tests of the progeny from the crosses considered 
in the preceding section the behavior of the double mutants of pyr Sa and 
orn 1, 2 and 3 suggested that they lacked the absolute pyrimidine require 
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ment which is characteristic of pyr Ja. These double mutants were 
therefore obtained from crosses of the three ornithineless mutants to pyr 3a 
in order to make it certain that » was not present. Results from tests of 
these recombinants in flasks appear in tables 3 and 4. In pyr 3a-orn 3 
the pyrimidine requirement appears to be completely suppressed on un- 
supplemented medium, since this strain grows in the same manner as does 
orn 3. On medium supplemented with ornithine, citrulline or arginine 
growth is somewhat slower than on minimal medium, but more rapid 
growth was obtained with enzyme-hydrolyzed casein as supplement. 
Pyr 3a-orn 1 and pyr 3a-orn 2, like orn I and 2, fail to grow on unsupple- 
mented medium but grow slowly without pyrimidine when supplied orni- 
thine or arginine, until they reach the dry weight characteristic of the 


rABLE 4 
REPRESENTATIVE DRY WEIGHTS (IN MG.) FROM 125-ML. FLASK CULTURES WITH 20 ML 
oF MEDIUM, AFTER 4 Days’ GROWTH WITH THE SUPPLEMENTS GIVEN BELOW 
0 Il MG 2 MG 1 MG SMG 
DL Ornithine Hydrochloride 
13 24 
11 22 
57 D4 
0 
tr 
tr 


ite 0.5 Mg. DL Citrulline 
3.5 6.5 1S 73 


65 


pyr 3a-orn 1 


a4 


} SY 82.5 63 
pyr 3a-cit 1 2 2.5 1.5 1.é 
33 31 32 


pyr Sa-orn 3 < 


3 
pyr 3a-orn 2 } 19.5 
id 
» 
> 


pyr 3a 30.5 32 


Nore: tr. (trace) represents a quantity of mycelium which is visible but too small to 


be removed and weighed 


single ornithineless mutant on that level of the supplement. They grow 
even more slowly with citrulline as supplement, as would be expected for 
reasons which are given in the next section. After 14 days, however, with 
t mg. of DL citrulline per flask, the following dry weights were obtained: 
pyr 3a-orn 1, 109 mg.; pyr 3a-orn 2, 115 mg. 

The difference in the rate of growth of pyr 3a-orn 3 on minimal medium 
and on medium supplemented with ornithine, citrulline or arginine sug 
gests that these compounds interfere with suppression of pyr 3a by orn 3 
as they do with suppression of pyr 3a by s. This is also suggested by the 
observation that growth of pyr 3a-orn 3 is more rapid with lower than with 
higher levels of ornithine. The effect of lysine on the pyr 3a-orn double 
mutants has not been thoroughly investigated, but it appears to be com- 
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plicated by inhibition of the orn mutants by lysine. <A significant stimu- 
lation of pyr 3a-orn 3 by 6 mg. of L lysine hydrochloride was observed when 
this strain was cultured in the presence of 2 mg. of DL ornithine hydro- 
chloride. With this ratio of lysine to ornithine pyr 3a-orn 1, pyr 3a-orn 2, 
orn 1 and orn 2 were completely inhibited whereas growth of orn 3 was 
less affected. Thus it seems that lysine may exert the same effect on the 
pyr 3a-orn double mutants as it does on pyr 3a-s but that this is obscured 
by its inhibitory effect on the orm mutants. 

It is considered unlikely that a change in the pyr 3a gene, such as back- 
mutation, is responsible for the apparent suppression of the pyrimidine 
requirement in these double mutants. The absence of an absolute require- 
ment for pyrimidine was detected in pyr-orn double mutants among the 
progeny of six crosses involving two different isolates of both pyr 3a and 
pyr 3a-s. Segregants carrying pyr 3a alone, on the other hand, even those 
derived from asei which also produced double mutants, showed the absolute 
requirement. Asa further check, the isolates of the double mutants which 
were used in the above experiments were crossed to wild type and pyr 3a 
recovered with its pyrimidine requirement unchanged. 

The effect of pyr 3a on citrulline utilization.-As may be seen from data 
in table 4, the utilization of citrulline by orn 1 and ? and cit 1 appears to 
be influenced by the presence of the pyr 3a gene. This effect is more 
marked in the case of pyr 3a-cit 1, which, in the presence of cytidine, gave 
a very meager response to citrulline at any level tested, whereas pyr 3a-orn 1 
and pyr Ja-orn 2 appear to have merely increased requirements. There is 
no effect in the case of pyr 3a-orn 3 nor has there been observed any effect 
upon ornithine utilization by orn 1, 2 and 3 nor upon arginine utilization 
by any of these strains. (Pyr 3a-cit 2 has not yet been obtained because 
of rather close linkage which exists between these two loci.) Cytidine has 
not been found to affect the utilization of citrulline by the single orn and 
cit mutants, nor does citrulline interfere with cytidine utilization by pyr 3a. 

It is of interest that, as table 4 shows, dry weights from the pyr 3a-orn 
double mutants cultured on cytidine plus citrulline were in some cases as 
much as twice that given by pyr 3a alone on the same amounts of cytidine 
and citrulline. This again demonstrates the partial suppression of the 
pyrimidine requirement in these double mutants. 

Other recombinants.-There are known at present, in addition to pyr 3a, 
three non-allelic pyrimidineless mutants, pyr J (263)," pyr 2 (388502)," '° 
and pyr 4 (36601)."! It has already been reported that pyr J and 2? are not 
suppressed by s and this has been found to be true also of pyr 4. These 
mutants, and pyr 3a as well, are inhibited by arginine when uracil is the only 


exogenous source of pyrimidine, but there is no inhibition if cytidine or 
uridine is added. No greater than half inhibition has been observed at 
any level of arginine tested, and there is no release by lysine. In double 
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mutants with orn 3, pyr 1, 2? and 4 were not suppressed and, in the presence 
of added cytidine, did not affect the behavior of orn J. 

The double mutants of pyr 3a with prol 2 and 3, when tested in flasks, 
behaved as would be predicted from the properties of the single mutants. 
No effect of either mutant gene upon the requirement introduced by the 
other was observed. Some isolates of pyr 3a-prol 1 behaved is this man 
ner, but others, which grew on complete medium, failed to grow when sup 
plied cytidine and proline. This has not been investigated further. 

It has been reported’ that 36703, a mutant which uses only arginine* 
showed no interaction with pyr 3a or s. More recently a histidineless 
mutant, CS4° was tested, with the same result. In linkage tests'? pyr 3a 


or one of its possible alleles, has been crossed to mutants requiring inositol, 


pyridoxin, pantothenic acid, choline, adenine, phenylalanine, isoleucine + 
valine and tryptophan, respectively, and no suppression of the pyrimidine 
requirement was observed. 


SUMMARY OF SINGLE MUTANTS AND RECOMBINANTS 
Single Mutants 


prol 1 slight growth on minimal—uses proline, 
prol 21 
prol 3) 
orn It 


grow slowly on minimal—use proline, ornithine, citrulline or arginine. 


al no growth on minimal—use ornithine, citrulline or arginine, 
orn « 


orn 3 grows slowly on minimal —uses ornithine, citrulline or arginine, 
cut 1! 
cit 2) 
pyr I slight growth on minimal — requires pyrimidine 


grow slowly on minimal—use citrulline or arginine, 


pyr 2 no growth on minimal at 35°C.—requires pyrimidine, 
pyr Sal 
pyr 4 
s no phenotypic differences from wild type observed, 


no growth on minimal-—require pyrimidine. 


Recombinants 


prol 1-s no phenotypic differences from prol 1 observed 
prol 2-s! 
prol 3-s\ 
orn >| 


no phenotypic differences from wild type observed, 


orn 2-s . . P 
no growth, or very slight growth on minimal—use citrulline or argi 


orn 3-S] - 
nine but not ornithine. 
cut 1-s 


cit 2-s 

pyr 1 s} 

pyr 2-s> no phenotypic differences from pyr 1, 2 and 4 observed. 
pyr 4 s\ 

pyr 3 a-s 

pyr 3a-prol 2-s > 
pyr 3a-prol 3-s 


growth on minimal approximates that of wild type— inhibition by pro- 
line, ornithine, citrulline or arginine relieved by lysine. 
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pyr 3a-prol 2! no growth on minimal— require pyrimidine + proline, ornithine, 
pyr 3a-prol 3% citrulline or arginine 


pyr 3a-orn ‘) no growth on minimal grow slowly without added pyrimidine on 


medium supplemented with ornithine, citrulline or arginine—in pres- 
( ence of cytidine the citrulline requirement is higher than that of 

pyr 3a-orn 2) orm 1 or 2. 

pyr 3a-orn 3 indistinguishable from orn 3 on minimal —grows slowly without pyrim- 
idine on medium supplemented with ornithine, citrulline or argi- 
nine~-with cytidine present the response to these compounds is 
like that of orn 3. 

pyr no growth on minimal or on medium supplemented with cytidine alone. 
In the presence of cytidine the response to arginine is like that of 
cit 1-s, but there is almost no response to citrulline. 

pyr l-orn 3 slight growth on minimal—uses cytidine + ornithine, citrulline or 
arginine 

pyr 2-0rn 3 no growth on minimal at 385°C. -uses cytidine + ornithine, citrulline 
or arginine 

pyr 4-orn 3 no growth on minimal uses cytidine + ornithine, citrulline or 
arginine 


Discussion. —\n order to simplify discussion of the observations presented, 
a summary of the various interactions involving pyr 3a and s is given below, 


The Effects of the Suppressor of Pyrimidineless: 


l It suppresses pyr 3a in the absence of added proline, ornithine, citrulline or argi- 
nine, or, in the presence of these compounds if lysine also is added.! 

2. It suppresses two genetically different mutants, prol 2 and 3, which require pro- 
line, ornithine, citrulline or arginine 

3. It suppresses the double mutants, pyr 3a-prol 2 and pyr 3a-prol 3, which, in combi 
nation with it, respond to proline, ornithine, citrulline, arginine and lysine as does pyr 
JSa-s. 

4. It prevents the slow growth characteristic of orn 3, cit 1 and cit 2 on unsupple- 
mented medium, 

5. It prevents the utilization of ornithine by the ornithineless mutants. 

6. It interferes with the utilization of a-amino adipic acid by the lysineless mutant, 
ly 1 (33933) unless arginine, citrulline or ornithine is added.? 

7. Incombination with pyr 3a it prevents utilization of a-amino adipic acid by ly 1. 

8. In combination with pyr 3a it interferes with the utilization of lysine by ly 3 
(4545) unless arginine is added.? 
Effects of Pyrimidineless 3a: 

1. It interferes with the utilization of citrulline by orm J and 2 and almost completely 
stops citrulline utilization by cit 1. 

2. It prevents the slow growth characteristic of c7t 1 on medium supplemented with 
cytidine. 

3. It interferes with utilization of a-amino adipic acid by /y 1 unless a small amount 
of lysine is added.? 


Effect of orn 1 and 2 on pyr 8a: 


They partially suppress pyr 3a on medium supplemented with ornithine, citrulline or 
arginine. 


eee 
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Effect of orn 3 on pyr Sa: 

It suppresses pyr 3a, completely on minimal medium and partially on medium supple 
mented with ornithine, citrulline or arginine. 

The pyrimidineless suppressor, then, acts as a suppressor of the recog- 
nized phenotypic effects of three genes, but when it is present in combi- 
nation with any one of seven other mutant genes, it acts as an inhibitor, 
in the sense that it prevents strains carrying these genes from doing some- 
thing which they are able to do if s is replaced by its “wild type’’ allele. 
Except for dhe difference in the degree of suppression, the “wild type’’ 
allele of s might be more justly called a suppressor. In any case there ap 
pears to be no reason for considering s to be in a different category from the 
other mutant genes involved in the interactions described here, since, for 
example, pyr 3a inhibits and the ornithineless mutants suppress. It is 
possible that ‘ 
as frequent as those which do not, but their detection may be less likely. 
In cases reported by Zalokar'’ and Wagner and Haddox'' suppression of one 


‘suppressors’ which themselves introduce requirements are 


mutant by a second was predicted on the basis of growth responses of the 
first and this suppression was observed when the double mutant was pre 
pared and tested. The double mutant, involving pantothenicless and 
tyrosine or phenylalanineless, described by Wagner and Haddox, is perhaps 
unique in that, although it does not grow normally in the absence of panto 
thenic acid, it is not stimulated by addition of this growth factor. 

It now seems unlikely that s takes “over the function of the inactive 
gene at the (pyr) 3a locus” as was suggested previously.' In fact it seems 
likely that the pyrimidine requirement of pyr 3a is not due to an “inactive 
gene,’ but rather to the activity of a gene being revised in a way such that, 
in a given genetic background, a pyrimidine requirement results. The re- 
quirement, the suppression of it and the other interactions as well, may re 
sult from alterations in the balance of reaction rates under the influence 
of the various mutant genes. The alterations may be pictured as arising 
through a variety of means, such as changes in enzyme concentrations or 
activities which produce changes in substrate or inhibitor levels affecting 
many reactions. The reactions involved may be related by a common 


coupled system or intracellular organization which is disturbed or modified 


in some way by the mutant genes. The variety of interactions observed 
does not point to any one mechanism as a very probable explanation, but 
rather, appears to suggest a greater number of possible mechanisms. 

When the interactions involving the lysineless mutants were first re 
ported it appeared to the authors that these, together with the observed 
accumulation of uracil by certain lysineless mutants,’ indicated a close 
relationship between lysine and pyrimidine biosynthesis. In view of the 
further complexities reported here, a direct interrelation no longer seems 
necessarily indicated. 
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The shifting of the position of the apparent ‘“‘genetic block’? by modi- 
fiers has already been reported by Haskins and Mitchell* in a paper dealing 
with tryptophan-nicotinicless mutants. This phenomenon, which is rep 
resented here by the action of s and pyr 3a on the ornithineless, citrul- 
lineless and lysineless mutants, suggests that the relationship between the 
mutation and the nutritional requirement may not always be as direct as 
it has sometimes been considered to be. As Haskins and Mitchell® have 
suggested, the assumption that the gene affected by the mutation controls, 
in a direct fashion, the reaction which, on the basis of the results of growth 
experiments appears to be blocked, does not seem to be a safe one. 

The fact that pyr 3a and its possible alleles, 3b and 3c are suppressed by 
s, whereas pyr J and 2 are not, was at first considered by the authors to be 
evidence that pyr 3a, b and ¢ were indeed alleles. At that time suppres- 
sors which acted on non-allelic mutants were known in Drosophila!’ and 
more recently Giles'® has found that certain non-allelic methionineless mu 
tants of Neurospora have a common suppressor and Lein and Lein" have 


reported that the same is true of three non-allelic acetateless mutants, 


also of Neurospora, It is, then, abundantly clear that response to the 
same suppressor cannot be considered as proof of allelism. Suppression 
by s of mutants having different requirements is reminiscent of the cases 
known in Drosophila’ in which the same suppressor acts upon mutants 
which are not only non-allelic but which also have different and seemingly 
unrelated phenotypic effects. 
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THE EFFECTS OF GENE CHANGE ON TRYPTOPHAN 
DESMOLASE FORMATION* 


By CHARLES YANOFSKY 
OSBORN BOTANICAL LABORATORY, YALE UNIVERSITY, NEW HAVEN, CONN 
Communicated by G. W. Beadle, December 29, 1951 


Genetic and biochemical investigations have provided substantial 
evidence for the view that physiological reactions are under the control of 
nuclear genetic material. These investigations have also indicated that 
the inability to carry out cellular biochemical reactions is a phenotypic 
expression of gene change and presumably results from gene-determined 
variations in the formation andor the characteristics of particular enzyme 
systems. 

Although the existing evidence suggests a close relationship between 
gene, enzyme and physiological reaction, the manner in which enzymes 
are controlled by genes and affected by gene change is poorly understood. 
In an effort to define more clearly the relationship between gene and 
enzyme, the present investigation was initiated. 

Characteristics of Mutant Strains C-S3 and S-1952. The two Neurospora 
mutants employed in this investigation, strains C-S5!' and S-1952,? are 
characterized by a tryptophan requirement which cannot be satisfied by 
indole. The biochemical reaction presumably blocked in the two mutants, 
the coupling of indole and serine, has been demonstrated to be the principal 
and probably the sole means of tryptophan synthesis in Neurospora.’ 
Since the enzyme catalyzing tryptophan synthesis, tryptophan desmolase, 
is extractable from wild type mycelium in an active cell-free state,‘ mutants 
C-83 and $-1952 constitute material in which the relationship of gene to 
enzyme can be investigated readily. 

The tryptophan requirement characteristic of mutants C-S3 and S-1952 
was demonstrated to be due to a single gene. From crosses of each mutant 
with several standard wild type strains, invariably |: 1 segregation of the 


tryptophan requirement was observed in all asci tested. Both the C-S3 


and $-1952 mutant genes are approximately the same distance from the 
centromere (C-S3 uncorrected centromere distance = 32 based on 51 asci, 
S-1952 uncorrected centromere distance = 83 based on 57 asci) and are 
linked to the fluffy locus (fluffy locus centromere distance = 29).°  Allelism 
tests were performed with the two mutants by Dr. Dorothy Newmeyer of 
Stanford University and of 6300 ascospores obtained by crossing the 
mutants to each other, no tryptophan independent ascospores were re 
covered. Attempts in this laboratory to form tryptophan independent 
heterocaryons between the two mutants were also unsuccessful. Pre- 
liminary reversion experiments employing ultraviolet light as a mutagenic 
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agent indicate that both mutant genes revert to tryptophan independence 
(T*) at about the same frequency (C-S83, 43 T*,2.9 & 10° viable macro- 
conidia; $-1952, 44 Tt/1.8 & 10? viable macroconidia). Together these 
observations suggest that the C-S3 and $-1952 mutant genes (henceforth 
designated as fd, and fds, respectively) are either allelic or closely linked. 
In the absence of contradictory evidence, they will be considered alleles 
in this paper. However, further experiments are in progress which are 
designed to examine more fully the linkage relationships between the 
two mutant genes and the nature of the reversions to tryptophan inde- 
pendence 

Since both mutants are blocked in the conversion of indole to tryptophan 
their culture filtrates were tested for possible indole accumulation. — Indole 
(tentatively identified on the basis of color reactions with p-dimethyl- 
amimobenzaldehyde and xanthydrol) and anthranilic acid as well were 
present in culture filtrates from both strains. Mutant S-1952 was con- 
sistently found to accumulate several times as much indole as did mutant 
C-S3. Whether or not this can be used as a distinguishing characteristic 
is difficult to tell at present since tryptophan dependent cultures obtained 
by crossing the two mutants to wild type strains show many different 
levels of indole accumulation. However, strains with the fd, allele always 
accumulate more indole than do strains with the fd, allele. 

The nutritional behavior of mutants C-S3 and $-1952 suggests that 
both strains lack the enzyme present in wild type strains which catalyzes 
tryptophan synthesis from indole and serine. To obtain direct evidence 
bearing on this possibility experiments were performed to determine 
whether or not extracts from mutant mycelium contain tryptophan 
desmolase activity. These tests were carried out with extracts and 
preparations from mutant mycelium which was grown under various 
culture conditions. The methods employed in obtaining tryptophan 
desmolase preparations are described in detail elsewhere.® For greater 
assay sensitivity a modification of the tryptophan desmolase assay pro- 
cedure’ was employed in these experiments. The incubation mixture 
contained 0.2 or 0.3 4M indole, 60 1.M DL-serine, 20 wg. barium pyridoxal 
phosphate,’ 0.06 ml. of 0.5 M phosphate buffer at pH 7.8 and enzyme 
extractin al ml. volume. Assay tubes were generally incubated for two 
hours at 87°. At the end of this period 0.2 ml. of 5°) NaOH was added 
to each tube to stop the reaction. Following the addition of base, 4 ml. 
of toluene were added to each tube. The tubes were then shaken and the 
layers allowed to separate. Aliquots removed from the toluene layer 
were assayed for indole,’ and when desired, tryptophan was determined’ 
in aliquots of the lower layer. The assay conditions employed are 
capable of detecting '/ 49) of the amount of tryptophan desmolase activity 
present in extracts from wild type strains. 
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Table | presents a summary of the tests carried out with mutant extracts 
and with various fractions obtained from these extracts. In addition 
to these tests, mutant extracts were incubated in the presence of pyridox- 


amine phosphate, yeast extract, casein hydrolysate, and with various 


rABLE 1 


PREPARATIONS FROM MurtTANTS C-83 AND S-1952 EXAMINED FOR TRYPTOPHAN Des 
MOLASE ACTIVITY 


AGE OF 
GROWTH MYCELIUM, FRACTIONS 
GROWTH SUPPLEMENT TEMP DAYS METHOD OF EXTRACTION TESTED 


Optimal T 2° 2,3,4,5 Homogenized in 0.1 CE, DCE 
30° : f MI phosphate at 
pH 7.8 
Optimal T 30 ‘ Homogenized in 0.1 Cell debris CE, 
M = phosphate at 30%, 50%, 
pH 7.8 70% (NH 4)b- 
SO, ppts 
Limiting T 30° ‘ Homogenized in 0.1 CE, DCE, PP, 
Optimal T M > phosphate at DPP 
Limiting T pH 7.8 
Optimal T + CH 
Limiting T+ 5 25 Homogenized in 0.1 
shaken M_ phosphate at 
pH 7.8 
Optimal T + 5 30 7 Mycelium extracted 
with water and as- 
sayed et various 
pH’s 
Optimal T 30° ‘ Extract of lyophi- Cell debris CE, 
lized mycelium DCE 
Optimal T : ya Homogenized in 0.1 CE, DCE, PP 
Limiting T M phosphate at 
Limiting T +S pH 7.8 
Limiting T +5 25° ‘ Homogenized in 0.1 
shaken M > phosphate at 
pH 7.8 
Optimal T 30° : Homogenized in 0.1 CE, DCE, PP, 
Optimal T + $ M = phosphate at DPP 
pH 7.8 
T = Trypotophan. S = Serine. CE = Crude extract. CH = Casein hydrolysate 
DCE = Dialyzed crude extract. PP = The partially purified preparation obtained by 
fractionating a crude preparation according to the procedure used with wild type 
extracts. DPP = Dialyzed purified preparation, 


concentrations of indole, serme and pyridoxal phosphate. Assays were 
also performed at several temperatures. Tryptophan desmolase activity 
could not be detected in mutant extracts under any conditions nor could 
activity be detected in any fraction obtained from these extracts. Mutant 
extracts were also mixed with wild type tryptophan desmolase prepara- 
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tions to test for the presence of inhibitors and no appreciable stimulation 
or inhibition was observed. 

These results are in agreement with those obtained previously by 
Mitchell and Lein'’ in a study of one of these mutants, C-S3. They too 
were unable to detect tryptophan desmolase activity in extracts from 
this mutant. The only report of the presence of tryptophan desmolase 
activity in mutant C-S3'!' has been retracted'’ as a result of further experi- 
mentation, At the present time, therefore, it appears that neither extracts 
from strain C-S3 nor from strain S-1952 contain tryptophan desmolase 
activity which can be detected by the methods employed. 

Although tryptophan desmolase activity could not be demonstrated 
in mutant extracts, 1t is conceivable that there is active desmolase in the 
mutants and that the enzyme is destroyed during extraction or is inactive 
under the assay conditions employed. To obtain some data concerning 


TABLE 2 
CROSSES INVOLVING THE SUPPRESSOR GENE 


ASCUS TYPE 
CROSS T OT 2T 41 


Chilton wild type * S-1952 “reverted” 6 
Lein wild type X $-1952 “‘reverted”’ 8 
Chilton wild type carrying the suppressor & $-1952 8 
Lein wild type carrying the suppressor * S-1952 ‘ 11 
Total 
Chilton wild type carrying the suppressor & C-83 
Lein wild type carrying the suppressor X C-&83 
Total 


= Tryptophan independent. TO = Tryptophan dependent. 


these possibilities an isotope experiment was performed to determine 
whether mutant C-S3, when growing, is capable of tryptophan synthesis. 
The details of this experiment have beén published previously.'* In 
brief, mutant C-S3 was grown on N!° ring labeled tryptophan and after 
Maximum growth was attained tryptophan was isolated from the mycelium 
formed. The isolated tryptophan had the same N'° content as the labeled 
tryptophan fed, showing that mutant C-S3 did not synthesize a significant 
amount of tryptophan under the conditions of this experiment. 

These data, then, all support the view that mutants C-S83 and $-1952, 
biochemically similar in their behavior, are both unable to synthesize 
tryptophan because of the absence of a normally essential enzyme, trypto- 
phan desmolase. Thus these two mutants would appear to represent 
strains in which the result of gene mutation is the loss of a particular 
enzymatic activity. 

Suppressor Studies. During the course of this work a microconidial 
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stock of strain S-1952 reverted to tryptophan independence. Miucro- 
conidia from this reverted culture were plated out and a homocaryotic 
tryptophan independent culture was recovered. This strain was_ then 
crossed with two standard wild type strains and as can be seen in table 2 
(crosses | and 2), three types of asci were obtained. The most frequent 
ascus class had 6 tryptophan independent and 2 tryptophan dependent 
ascospores. This finding suggests that the reversion to tryptophan 
independence was due to a genetic change at a separate locus (su). This 
change appears to result in the suppression of the tryptophan requirement 
associated with the fd, allele. As a further check of this conclusion trypto- 
phan independents recovered from 4:4 aser (4T*:41T~) were crossed to 
mutant $-1952 and again three types of asci were obtained (table 2, crosses 
3 and 4). 


FIGURE 1 


Growth of various strains of 

FIGURE 2 
Neurospora on different levels 
of tryptophan,  Incubated for 72 Growth curves of various strains on mini 


hours at 30°, mal medium at 30 


The cultures obtained from these crosses were scored by inoculating 
two flasks with conidia, one flask containing minimal medium and _ the 
other flask containing minimal medium plus | mg. L-tryptophan. These 
flasks were incubated for 72 hours at 30° after which the mycelial pads 
were removed, dried and weighed. By employing this procedure it is 
possible to distinguish the four possible types of cultures: mutant, sup- 
pressed mutant, wild type and wild type carrying the suppressor gene. 
The mutant and wild type cultures can easily be identified since the former 
cultures do not grow on minimal medium while the latter cultures grow 
well on both types of media. The other two types of tryptophan inde- 
pendent cultures can be identified also, since the suppressed mutant grows 
poorly on minimal medium (in a 72-hour period) and better in the presence 
of tryptophan, while the growth of wild type strains carrying the sup- 
pressor gene 1s poor on the same two media and is not stimulated by 
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tryptophan. Comparisons of the growth of the various types of strains 
are shown in figures | and 2. In figure 2 it is seen that there is a lag in the 
growth of strains carrying the suppressor gene. Attempts to overcome 
this lag or inhibition associated with the suppressor gene by adding various 
supplements have been unsuccessful. 

To obtain some information concerning the specificity of the suppressor 
gene, wild type strains carrying this gene were crossed to mutant C-S3 
As may be seen in table 2 (crosses 5 and 6) only 4:4 asei (47T*:4T ~) were 
obtained. Thus the effect of the suppressor gene appears to be specific 
for the fd, allele. ‘That some of the tryptophan dependent cultures derived 
from these crosses actually carry the suppressor gene was demonstrated 
by performing crosses with strain S-1952. Additional crosses are being 

performed to determine if a par- 
ticular genetic background is neces- 
sary for the suppressor gene to show 
its effect on the fd, allele. Green'4 
has also reported a case in which 
a suppressor gene does not affect 
every member of an allelic series 
C-S3 strains carrying the sup- 
pressor gene can be distinguished 
from C-S35 strains without this gene 
since the presence of the suppressor 
gene results in a marked inhibition 


of growth (see Fig. 1). It can also 
hours 


; be seen in figure | that C-S3 strains 
FIGURE 3 


a eee require less tryptophan to reach 
stat nthe eutan deste axinun growth when they carr 
For the unit of enzyme activity see (6). the suppressor gene. However, 
since the growth of strains with the 
td\-su genotype is poor under all conditions, it is difficult to evaluate this 
observation. 

Examination of Strains Carrying the Suppressor Gene for Tryptophan 
Desmolase Activity. Since S-1952-su strains grow in the absence of exoge- 
nous tryptophan, enzyme studies were performed to see if this ability could 
be correlated with the presence of tryptophan desmolase. Also examined 
for tryptophan desmolase activity were two wild type strains carrying the 
suppressor gene and a C-S3 strain carrying the suppressor gene (trypto- 
phan dependent). Tryptophan desmolase activity was present in extracts 
from strain S-1952-su as well as in extracts from wild type strains carrying 


the suppressor gene. No activity could be detected in extracts from strain 


C-S3-su. Attempts to activate extracts from this latter strain were also 
unsuccessful (see table 1). Thus the suppressor gene when present with 
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the fd, allele restores both the ability to synthesize tryptophan desmolase 
and to grow in the absence of exogenous tryptophan; neither of these 
effects is observed in strains with the fd)-su genotype. 

An effort was made to compare the amount of tryptophan desmolase 
in strain 5-1952-su with the amount in wild type strains with and without 
the suppressor gene. The tryptophan desmolase content of wild type 
mycelium varies considerably with time as can be seen in figure 3 and, 
in view of the growth behavior of strains 5-1952-su and wild type carrying 
the suppressor gene (see Fig. 2), such a comparison is difficult. Accepting 
these uncertainties, the maximum amount of tryptophan desmolase found 
in a suppressed S-1952 strain was approximately |.) of the maximum 
amount of tryptophan desmolase that is formed by wild type strains. 


FABLE 4% 


COMPARISON OF THE EFFECT OF VARIOUS INHIBITORS ON TRYPTOPHAN DESMOLASE FROM 
STRAIN S-1952-su AND FROM A WILD TYPE STRAIN 


© INHIBITION 

SUPPRESSED 

FINAI MUTANT WILD TYPR 

ADDITION CONCENTRATION, MOLAR DESMOLASE DESMOLASE 

Hydroxylamine ] 10 100 Q? 

I 1 1 2 3 

Co* l 10 73 76 

CN 1 < 107? Hi) HA 

Zn** | 10 7s x] 
Mg’ l 10 i) 
NH,Cl { 29 
NH,Cl 37 44 

L-tryptophan <x 107? 23 19 

L-tryptophan 2X 107 36 29 


Thus the presence of the suppressor gene restores to strain S-1952 the 
ability to form effective concentrations of tryptophan desmolase but does 
not permit tryptophan desmolase synthesis to proceed to the same extent 
as it does in wild type strains. The stimulation of growth of strain $ 
1952-su by low concentrations of tryptophan (see Fig. 1) 1s in agreement 
with this finding. Wild type strains carrying the suppressor gene have a 
maximum tryptophan desmolase content which is about one-third that of 
wild type strains without the suppressor gene. 

Comparison of Wild Type Tryptophan Desmolase and Tryptophan Des- 
molase from Strain S-1952-su. Several characteristic properties of the 


tryptophan desmolase from a wild type strain were next compared with 


the properties of tryptophan desmolase from = strain S-19052-su. On a 
molar basis indole uptake and tryptophan formation are equivalent in the 
presence of preparations from either strain. The stability of both prepara 
tions is roughly the same, both losing approximately 50°, of their trypto- 
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phan desmolase activity when stored at 2” for 24 hours. Both preparations 
are inhibited by the same substances, and to practically the same extent 
(see table 3). As was the case with wild type preparations,’ indole was 
not taken up by suppressed-mutant tryptophan desmolase preparations in 
the presence of the following substitutes for serine: alanine, glycine, 
cysteine, cystine, methionine, threonine, aspartate, glutamate, acetate, 
pyruvate, succinate or malate. 

Before other comparisons were performed an attempt was made at 
fractionation and partial purification of suppressed-mutant tryptophan 
desmolase preparations according to the procedure developed for wild 
type preparations.® This procedure involves the use of various concentra- 

tions of ammonium sulfate at two pH values. 
Suppressed-mutant preparations and wild type 
preparations behaved identically throughout the 
steps in this procedure. The partially purified 
(about 20-fold) suppressed-mutant tryptophan 
desmolase preparations which were obtained were 
then employed in further comparisons. The pH 
activity curves for such a suppressed-mutant prep 
aration and for a wild type preparation are 
shown in figure 4. Both preparations exhibit maxt- 
rome « mum activity at pH 7.5. 

pH_-activity curves for Comparisons of the affinity of the two partially 

tryptophan — desmolase — purified enzyme preparations for their substrates 


preparations froma wild and for their coenzyme were also performed. 
type strain (@ -@) and ‘ 
from strain S-1952-su 0 — : ce ‘ 
(O-O) | hour incu. for serine, excesses of indole and pyridoxal phos 


In the determinations of the Michaelis constant 


bation at 37 2 x phate were employed and indole uptake was 
1) * M indole, 6 X 10°? measured. A Lineweaver-Burk plot'® of kinetic 
M DL-serine, 20 4g. bar- data for serine is shown in figure 5. The approxi- 
ium pyridoxal phosphate : , ; ‘ 

mate A,, value for serine for both preparations can 
be estimated from figure 5 to be 5.4 & 107% JAZ, 
This value agrees fairly well with the value of 6 & 107% A/ found for wild 
type tryptophan desmolase preparations under different conditions.* 
Pyridoxal phosphate saturation curves are shown in figure 6. The 


per ml 


pyridoxal phosphate concentration required to restore half-maximum 


activity is approximately the same for the two enzyme preparations, 3. X 
10 ® MM. It ean also be seen from figure 6 that neither enzyme prepara 
tion (both were partially purified) was active in the absence of pyridoxal 
phosphate. The velocity of indole uptake by both preparations is es 
sentially zero order at indole concentrations higher than 2 & 107! AZ, 
The energy of activation was determined for partially purified wild 
type and suppressed-mutant preparations. Excesses of indole and 
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pyridoxal phosphate and a limiting concentration of serine were employed. 


Indole uptake was measured and the corresponding decrease in serine 
concentration was calculated and is plotted in figtre 7 (a). The velocity 
constants obtained from these curves were then used to determine the 
energy of activation. The log k& against 1/7 plot is shown in figure 7 (}). 
In this figure it can be seen that the points obtained with the two enzyme 
preparations fall on the same line. The energy of activation, therefore, 
is approximately 12,000 cak./mole. in each case. 

Thus the tryptophan desmolase preparations obtained from both strains 
appear to be similar by the criteria applied. 

Discussion..-The data presented in this paper permit certain conclusions 


concerning the relationship between gene and enzyme. ‘The studies with 


FIGURE 5 

Lineweaver-Burk plot of gs 
kinetic data for serine. 1 FIGURE 6 
hour ee at P Wild Pyridoxal phosphate saturation 
Ee , . ; 
i Phaser 6 : sharers curves for tryptophan desmolase 

. »- 1952-5 ry ‘ : : ; 

: Ss — u tryptophan preparations from a wild type strain 

> » ( ) am ot 4 : . . 
: esmolase ; o Reka 10 (@--@) and from strain S-1952-su 
M indole, 20 yg barium pyri (O--O), 1 hour incubation at 37° 
doxal phosphate per mil 3 xX 10-4 M indole. 6 X 1072 M 


I) L-serine 


mutants C-83 and S-1952 have shown that genetic alteration has affected 
the ability of each strain to produce active tryptophan desmolase. The 
nature of the effect on tryptophan desmolase or tryptophan desmolase 
formation which is responsible for the absence of detectable activity 
cannot be stated. It is clear, however, that if tryptophan desmolase is 
present in the mutants it does not provide sufficient tryptophan to support 
growth. 

Although strains with the fd, allele (S-1952) lack tryptophan desmolase 
activity, they are capable of forming active tryptophan desmolase when 
they are carrying a specific form of a different gene. This second gene, 
however, although it enables strains with the fd, allele to form active 
enzyme, does not appear to be capable of mediating tryptophan desmolase 
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synthesis by itself. This conclusion is drawn from the observation that 
strains with the /d/,-su genotype do not have tryptophan desmolase activity. 
In strain S-1952-su then, two genes, both in specific forms, appear to be 
required for the synthesis of tryptophan desmolase. 

In addition to this case in which two specific genes appear to be required 
for the formation of a single enzyme, there are other cases in which muta 
tion of any one of a number of genes affects the formation of single en 
zymes.'®'* ‘Together, these findings might be interpreted as indicating 
that the relationship between genes and enzymes is exceedingly complex 
However, since so many cellular reactions are interrelated and inter 
dependent it is difficult to tell at which level the complexity exists. 


FIGURE 7 


The effect of temperature on the rate of serine uptake (indole up 
take was measured) by tryptophan desmolase preparations from a 
wild type strain (OO) and from strain S-1952-su(@--@), 2 10°! 
M indole, 4 & 10°? V7 DL-serine, 380 wg. barium pyridoxal phosphate 
per ml 


In view of the fact that two genes, each differing from its wild type allele, 
are required for iryptophan desmolase synthesis in strain S-1952-su, one 
might expect the tryptophan desmolase formed by this strain to have 
altered properties. In a preliminary comparison in which a number of 
enzyme characteristics were examined, no major difference was found 
between the tryptophan desmolase extracted from strain S-1952-sa and 
the same enzyme extracted from a wild type strain. To be sure, the 
methods employed could not reveal all types of changes that could occur 
in the suppressed mutant desmolase, nor are they sensitive enough to 
detect very slight changes. However, the catalytic portion of S-1952-su 


tryptophan desmolase, at least, does not appear to be altered appreciably. 


The question of the effect of gene change on enzyme specificity has been 


ee eae ee 
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considered in other investigations'®'’ and, as yet, there is no clear demon 
stration of a change in enzyme specificity as a result of gene mutation. 
The recent findings of Pauling, ef a/.,"’ however, do suggest that such 
changes may occur. In any case, since relatively few studies of the effects 
of gene change on enzyme specificity have been performed, we cannot 
exclude alteration of enzyme specificity as a possible result of gene mutation. 

It is difficult to account for the fact that the suppressor gene permits 
mutant $-1952 to grow in the absence of an exogenous supply of tryptophan 
and to form apparently normal tryptophan desmolase. The suppressor 
gene does not appear to be a duplication at another locus of the normal 
allele of the S-1952 locus since strains with the /d)-su genotype lack trypto 
phan desmolase activity. Nor does it seem likely that the suppressor 
gene acts by creating an alternate biosynthetic pathway of tryptophan 
formation or by reheving strain S-1952 from a tryptophan desmolase 
inhibitor. Although it is relatively easy to rule out most suggested 
explanations at the substrate level, at levels closer to the gene many 
hypotheses are consistent with the results obtained. ‘Thus, the suppressor 
gene could act by creating an alternate biosynthetic pathway of trypto- 
phan desmolase formation, by removing or inhibiting a reaction competing 
for a common enzyme precursor, by restoring the specificity of tryptophan 
desmolase, tryptophan desmolase precursors or the tryptophan desmolase 
forming system, by increasing the concentration of tryptophan desmolase 
precursors, tryptophan desmolase forming system or both or by affecting 
in any of a number of ways the rate of formation or time of appearance of 
tryptophan desmolase. It thus becomes apparent that we must have a 
better understanding of the processes involved in enzyme formation if 
we are to explain the mechanism of gene action. 

The results obtained in this study also demonstrate that the fd; and fd, 
alleles behave differently. This can be seen from the fact that the sup 
pressor gene, in contrast to its effect on strains with the fd, allele, does not 
enable strains with the fd; allele to grow in the absence of an exogenous 


supply of tryptophan, nor does it permit these strains to form detectable 


tryptophan desmolase. ‘Thus the behavior of the fd; and (dy, alleles in the 
presence of the suppressor gene suggests that these alleles may be quali 
tatively different forms of the same gene. In view of this and since there 
are several cases in which the complexity of single genetic loci has been 
demonstrated,’! ” it is tempting to speculate that each gene may be 
composed of one or more units each of which is concerned with different 
steps in the synthesis of a single enzyme. 

Summary. Two allelic tryptophanless Neurospora mutants, both 
unable to use indole for growth, were examined for tryptophan desmolase, 
the enzyme which catalyzes tryptophan synthesis from indole in wild 
type strains. Tryptophan desmolase activity could not be detected in 
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extracts or in various preparations from either strain. One of these 
mutants, S-1952, does form tryptophan desmolase and grow in the absence 
of an exogenous supply of tryptophan when it is carrying a specific sup- 
pressor gene. This suppressor gene does not similarly affect the trypto 
phanless strain (C-S3) with the allelic mutant gene. Since C-S3-su 
strains do not form tryptophan desmolase, it appears that both the S-1952 
mutant gene and the suppressor gene are required for tryptophan desmolase 


synthesis in strain S-1952-su. The tryptophan desmolase formed by a 


5-1952-su strain was compared with the tryptophan desmolase of a wild 
type strain and no major differences were detected. 
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In Euclidean /*” we will say that an operator with constant coefficients 
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of any even order is e/liptic semihomogeneous (in a generalized sense) if 
it possesses a Green's function G(x) which is defined and analytic in i" 
“?’ 


except at the origin, and if there is an ro such that for any 7, > 0, 
rr=> Or +... +4, => 7% we have at infinity 


Rt 


: — = O(|x| ~*~"), xi > o, (2) 
ON: *. .: OSy, 


We have recently demonstrated! that such a Green's function exists if the 
polynomial 


er ee (ity) 


and both its highest and lowest homogeneous parts occurring, are all 
three different from zero for any real 4, ..., ¢, not all zero. 

THEOREM |. Jn E:,, if Dts a bounded domain having a connected boundary 
B, and if U is a neighborhood of B, and if an analytic function f(x,;) in U 
satisties an elliptic semthomogeneous equation 


Af = 0 (3) 


in all its variables, and also an arbitrary (non-homogeneous, non-elliptic) 
equation 

oO! + + dm f 
“Ox,". . .dx,™ 


Af — ai } Lass i b, 


in fewer than all variables, 
m<n, 


the latter having again constant coefficients; then the function has an analytic 
continuation into all of D + LU’. 

Proof: This theorem is an outgrowth of the theorem of Hartogs for 
complex variables, and it was already previously’? given in the present 
version, and in more refined ones, but only for the operator (1) being the 
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ordinary Laplacean and the operator (4) being elliptic in its variables, 
more or less. 
The Green's function of (1) gives rise to a Green's integral 


g(x) = Sp Daas (Xa: 0%), 
—1)*~ “dg,. . dé, 1d&,41.. db, 


al ck 8) (sa ( é (6 
<2h—-1 : 7 io = Jy ») 
Y Of," dé,*” BI » 
where (i) each G® ,..() is a linear combination with constant coefficients 
of derivatives of G(é), (11) for fixed x and any solution /(£) of (3) we have 


OX, 4 4 OX, 
Of; 0, 


= () (7) 


and (1) if By bounds a domain in which /(£) is a solution of (3) then for 
x in the domain we have g(x,;) = /(x,), that 1s 


f(x;) = Sm dwXea'O. (S) 


Now, for the proof of the theorem we may assume that B is a piecewise 
differentiable simplicial closed surface, and if we take a layer around it 
which is bounded by two “parallel” surfaces, and apply (S) to the layer 
then we obtain in l’ a Laurent decomposition 


T(x,) ' h(x,) ~ fo(x i (9) 


where /; is analytic in D + Ll and fy, in (4" — D) + U, and where for x 
in D, or in ke” — D, respectively, both are represented by the “same” 
integral (5). Our next step is to show that both are again solutions of the 


“accessory? equation (4), 
= 0, Afe = 0, (10) 


and as previously® we do this in the following manner. We view the 
integral (5) as an operator, obviously distributive, from solutions of the 
elliptic equation (3) in the neighborhood of B to functions g(x) in k" — B, 
denoting it then by g(x) = L(/(); x), and due to the fact that in (6) the 
variable x occurs only in the combinations x — €& this operator is com 
mutative with partial differentiations as well, so that for any A with 
constant coeflicients we have 


A, g(x) L(A, f(E); x). (11) 


If in particular A/(é) = O then also Ag(v) = O and hence (10). It will be 
worth while to state the reason behind (11) in a comprehensive lemma, and 


af 


we again introduce the differential form o“” which for a < @ 1s defined as 
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— {err I 


dé dé LE +1 dé, és 4 dE, 


a 


for a = Bis 0, and for a> Bis —oa*". 


LemMMA |. // B is a (non-closed) surface with boundary C, and if in a 
neighborhood of B + C we are given a solution {(&) of an elliptic equation 
(3), and if we introduce the functional (11); then for any differential operator 


(4) with constant coefficients we have 


A,g(x) — L(Agf(8); x) = Sede poiYage™, 


where each Vy, ts a bilinear form in partial derivatives of G(x £), and in 
partial derivatives of {(&) of order < (2h — 1) + (M — 1). 

For the conclusion of Theorem | it suffices to show that the function 
fo in (9) is analytic everywhere. We have already obtained Af, = 0, 
and from inspection of the integral (5) it follows that at infinity /. has the 
same behavior (2) as G(x) itself, with another ro of course. Thus our 
conclusion follows from the following lemma which can be proved by 
Fourier analysis. | 

LemMa 2. [fan analytic function {(x,) in-an exterior of a sphere |x; >R 
satisfies an equation (4) in fewer than all variables, and at infinity has the 
behavior (2), then f(x,) ts a polynomial of degree <\J — 1 in all its variables 
and ts thus analytic everywhere. 

For the proof of the Laurent decomposition (9) instead of applying (8) 
to a layer around all of B we might have used as previously® the following 
“local” argument. 

LemMa 3. Jf we consider the integral (5) for any union of simplices B 
and any {(&) on it, and if ina neighborhood of point of a simplex, f(&) ts a 
solution of (3), then in the neighborhood of this point the function g(x) can be 
continued analytically across B from both sides, and the difference of g(x) on 
both sides 1s f(x) itself, 


= f(x). 


With this lemma the following sharpening of Theorem | ensues. 

THEOREM 2. Jf the piecewise differentiable surface B is a union of the 
closures of differentiable simplices By, ., B,; if ina neighborhood L’, of 
each (non-closed) B, there is a given analytic solution [, of (3) and (4); 4 
each f, together with its partial derivatives of order <(2h — 1) + (M — 1) 
has continuous boundary values on the boundary of B,; and if on each (n — 2) 
dimensional edge cemmon to two simplices all these boundary values are the 
same for the two functions; then in the domain D bounded by B there exists an 
analytic function {(x) which in each L’, coincides with the function f,(x) given. 

Another conclusion is as follows. 

THEOREM 3. /f a surface B with a boundary C does not decompose I", 


and if ina neighborhood of B there is given a solution of (3) and (4) for which 
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the partial derivatives (or only their boundary values) of order <(2h — 1) + 
(M — 1) are 0 on C; then f(x) itself is = 0. 
A full account with results in other directions will be given elsewhere. 


'“Zeta Functions and Green's Functions for Linear Partial Differential Operators 
of Elliptic Type with Constant Coefficients,”” Ann. Math. (to appear). 

* “Analytic and Meromorphic Continuation by Means of Green’s Formula,” Jhid., 
44, 652 673 (1943) 


THE DIRICHLET PROBLEM FOR LINEAR ELLIPTIC EQUATIONS 
OF ARBITRARY EVEN ORDER WITH VARIABLE COEFFICIENTS 
By Fe.ix E. BROWDER 
BOSTON UNIVERSITY 


Communicated by Marston Morse, January 22, 1952 


Let A be an arbitrary linear elliptic partial differential operator of order 
2m on a bounded domain D of Euclidean n-space. The Dirichlet problem 
for K on D asks for a solution u of the equation Au = 0 in D such that 
u and its derivatives of order less than m assume the same values on the 
boundary of D,in a suitable sense, as a prescribed function g and its corre- 
sponding derivatives. The coefficients of A being suitably differentiable, 


it is shown that if the Dirichlet problem for A has at most one solution, 
then it has exactly one for every boundary value problem. More generally 
the boundary value functions for which a solution to the Dirichlet problem 
exists are characterized in terms of the null solutions of the adjoint equation 
by consideration of the inhomogeneous equation. Finally a theory of 
eigenfunctions and eigenvalues is established for the eigenfunctions of a 
large class of self-adjoint equations. 

The proofs combine essential features of the method of orthogonal pro- 
jection with the theory of completely continuous linear operators on a 
Hilbert space. They involve no assumptions of self-adjointness or posi- 
tivity on A in contrast with the usual proofs based on the method of 
orthogonal projection or the Dirichlet principle. 

These results are a generalization of previous general existence theorems 
based on variational methods including those of Courant, Friedrichs, and 
Rellich for certain classes of self-adjoint second-order linear elliptic equa- 
tions,'! Friedrichs for the biharmonic equation,’ Soboleff and Vischik for 
the polyharmonic equations,’ Kodaira for the positive homogeneous self- 
adjoint second-order linear elliptic equation,‘ and Garding for general 
linear elliptic equations of arbitrary order with constant coefficients.® 

|. If 7 is a positive integer, C’(D) is the set of j-times continuously 
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differentiable functions on D, C?(D) the subset of C’(D) of functions vanish 
ing outside a compact subset of 7). 
For f ¢ C’(D), g « C’(D) we define a j-inner product by 


; in oO’; o’g ) 
(f, g)) = ix 
sii J ce Ox;,.. Ox, Ox,,...Ox,,/- 


and a j-norm $f; = (/, f);*. More generally the j7-norm may be defined 
for any geC’(D) whose jth derivatives are in L°()). The o-norm and 
inner product are those of L°(D). 

There exist constants ¢,; > O for each 7 < m for which (fo; < 6) filjqi 
for all fe C’”'(D). 

If // is the completion of C’”(D) with respect to the m-inner product, the 
Hilbert space // may be identified with a linear subset of L°(D). Thus 
identified, /7 will be designated as the set of functions of L?()) which vanish 
with all their derivatives of order less than m on the boundary of D. 

The Dirichlet Problem: Given geC"(D) for which (g , < ©, to find 
ueC*"(D) such that Ku = 0 in D and u-g and its derivatives of order 
less than m vanish on the boundary of D. 

K is said to be suitably differentiable if the coeflicients in A of the jth 
order derivatives belong to C’”"*’(D) and if the mth derivatives of all the 
coeflicients are continuous on the closure of D.“ The formal adjoint K 
of A is well defined for such A. 

K may be written in the forn 

- or" 
lit 2 Pi. ima. hm) OX ;,. . Og Wlp,. . . Nem re 


km 


where the indices 1, » Im, Ri, ..-, Rm range independently from 1 to n, 


the functions ~,;,  ;,,4,...%,.(%) are symmetric in the indices, and R is a 


differential operator of order less than 2m. 
K will be termed elliptic if both of the following two conditions are 
fulfilled : 
(a) The form 
Pm 5 Py 


is positive definite in the &-variables for each x in D. 
(b) If we define (f, h) for f, he C”(D) by 


Pes oO” Oh 
(f,h) = ( ee sn (X dv) 
f df ee eres Ox, OXs, ONz,. « - bbe 


then there exists a constant c > Osuch that (f, f) 2 ¢(f,f),, forall fe C"(D) 

If A is of the second-order or if the terms of highest order have constant 
coefficients then condition (>) follows from (a) which is the usual ellipticity 
condition, 
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Theorem |. Jf K is a suitable differentiable linear elliptic differential 
operator on a bounded domain D in E" and if the Dirichlet problem for zero 
boundary values has only the zero solution, then the Dirichlet problem has a 
solution for every boundary function g. 

An important part in the proof of Theorem | is played by the following 
generalization of the well-known Weyl Lemma. 

THeoreM 2. Let K bea linear elliptic operator for which the coefficients 
of the jth derivatives belong to C°"*?(D). Suppose that v is a measurable 
function on D square summable on every compact subset of D, y « C'(D). 


Suppose also that for all f ¢ C?"(D) 
Sn vKf dx = Sow dx. 


Then v is equal almost everywhere in D to u where ue C*"(D) and Ku = y 


in D. 

The proof of Theorem 2 is based on a method which has been applied 
under various hypotheses by H. Weyl,’ Kodaira,* De Rham and Kodaira,* 
L. Schwartz’ and L. Garding.'!” The existence of a fundamental solution 
in the small for linear elliptic equations with analytic coefficients as estab- 
lished by F. John'' is used in the proof. 

Proof of Theorem 1: By the Schwarz inequality there exists a constant 
d > O for which 

fi| < di film 


for all f « C"(D) where | fii = (f, Sf). Thus (f, h) may be 
extended to an inner product on // under which // is a Hilbert space with 


a norm equivalent to the m-norm. 
Let L be the self-adjoint linear elliptic differential operator defined by 


Q” a 
» OX;,.. ae? tr. «mk OX;,. - -OXkm 


t tm 


hi. kem 
Then K = L + A, A of order less than 2m. 
Consider (h, g) for A in 77 and the given boundary value function g. 
By the Schwarz inequality 
(hy g)| <a) hin! iim a >O 


and hence (h, g) is a continuous linear functional for 4 in /7. Therefore 
there exists Ap € /7 such that (h, ty) = (h, g) for all h in /7. 
Setv = g — ho. Thenv — g = —hy eH and (h, v) = O for all A in 1. 
In particular if fe C’"(D), 
0 = (f,v) = (- 1)" Sp vLf dx 
By Theorem 2 there exists vw €C°"(D) for which Loy = 0 and w% — g €/1. 


In addition ||v9\m < ©. 
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The existence of a solution uw of the equation Au = O assuming the 
boundary values of g is equivalent with the existence of a solution w of 
the equation 

Kw = K(u — ) = —Atm 
with zero boundary values. By Theorem 2, w belongs to C?"(D) 1 H 
and is a solution of the above equation if and only if 

(1) SpwKfdx = — fpArof dx for all f « C?”"(D). 

(1) can be written in the form 
(2) fwhfdx + fwd fdx = f Auf dx. 
For we H, Jpwlf dx = (—1)"(w, f). 
Secondly, there exists p, > 0 for which 
|| - SpAvof dx|| < pil! fl] m 
for all fe C?"(D). As a consequence there exists » ¢ /7 such that 


— Sf Awf dx = (n,f) 


for all fe C°”"(D). Finally it may be shown by a simple argument that 
there exists a constant p. > 0 for which 


| [wAf dx| < P2|| | m 1 f m < P2x€ m || 101 m f m 


It follows from these last inequalities and the Rellich Selection Theorem!’ 
that a subset of /7 bounded in the m-norm is conditionally compact in the 
m-\-norm and consequently that there exists a completely continuous 
linear mapping |’ of // into itself for which 
(Vw, f) = SpwAf dx 
for all fe C°”"(D). Thus (2) can be written as 
((—1)"w + Vw — 9, f) = 0 

for all fe C?"(D). 

A necessary and sufficient condition that this latter equation hold is 
that 


w+ (—1)"Vw = (—1)"n 


for some wel. (—1)”"V is a completely continuous linear mapping of 
II into itself and by a theorem of F.  Riesz'* on completely continuous 
linear transformations of a Banach space / + (—1)”V is onto if the equa- 
tion 


w+ (—-1)"Vw = 0 


has no solution in // except w = 0. 
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Any such null solution, however, would also be a solution of the equation 


0 = (—1)"(w, f) + (Vw, f) SpwLlf dx + SpwAf dx 
SpwKky dx. 


Thus 


Spwki dx == () 


for all fe C?"(D) and, by Theorem 2, w which belongs to // would be 


almost everywhere equal to a 2m-times differentiable solution wy of the 
equation Aw) = 0. But wy is a solution of the Dirichlet problem with 
null boundary values and must be identically zero. Hence w = 0 almost 
everywhere in D and is the 0 element of /7. Therefore the Dirichlet prob- 
lem for K on D with boundary value function g has a solution. Q.E.D. 

CoROLLARY TO THeoreM |: Let A be a suitably differentiable linear 
elliptic differential operator on a bounded domain D in I". Suppose that 
there exists a constant p; > O such that 


(—1)" [Kf -f dx > ps\fi\m 


for all fe C°"(M). Then the Dirichlet problem for K has one and only one 
solution for every boundary value function g. In particular the above condt- 
tion holds for a homogeneous self-adjoint differential operator like the L of 
the proof of Theorem 1. 

THeoreM 3. Let K bea suitably differentiable linear elliptic operator on 
a bounded domain D in BE", pe L?(D) n CD). The equation Ku = yp 
has a solution u « C?"(D) with zero boundary values if and only if 


Srvedx = 0 


for each null solution g of Ke = 0. The linear dimension of the set of null 
solutions for K is finite and equal to that for K. 

The methods developed in this note can be used to treat the eigenfunc- 
tions and eigenvalues of the Dirichlet problem for A. With suitable 
restrictions on the coefficients at infinity the theorems may be extended 
to unbounded domains. 


' Courant-Hilbert, v. 2, Chap. VII. 

2 Math. Ann., 98 (1927). 

§ Math. Sbornik, 2, 44 (1937); Tbid., 25 (1949). 

* Proc. Imp. Acad. Tokyo, 20 (1944). 

5 Compt. rend. Paris, 230 (1950) 

6 The second condition though not the most economical possible is introduced for 
the sake of simplicity 

7 Duke Math. J., 7 (1940), 

8 Harmonic Integrals, Notes on Lectures at the Institute for Advanced Study, 
1950, pp. 36-39 

® Theorte des Distributions, v. 1, Paris, 1950, pp. 136 ff. 
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" Proc. Roy. Physiog. Soc. Lund., 20, 250-253 (1950) 
Comm. on Pure and Applied Math., v. 3 (1950) 

2 Courant-Hilbert, v. 2, pp. 488-490 

‘8 Acta Math., 41 (1918). 


INTEGRATION OF STIFF EQUATIONS* 
By C. F. Curtiss AND J. O. HIRSCHFELDER 
THe NAvAL RESEARCH LABORATORY, DEPARTMENT OF CHEMISTRY, UNIVERSITY OF 
WISCONSIN, MADISON, WISCONSIN 
Communicated by Farrington Daniels, December 29, 1951 

In the study of chemical kinetics, electrical circuit theory, and problems 
of missile guidance a type of differential equation arises which is exceed- 
ingly difficult to solve by ordinary numerical procedures. A very satis- 
factory method of solution of these equations is obtained by making use of 
a forward interpolation process. This scheme has the unusual property 
of singling out and approximating a particular solution of the differential 
equation to the exclusion of the manifold of other solutions. This behavior 
may be explained by a simple geometrical interpretation of the significance 
of the forward interpolation process. The differential equations to which 
this method applies are called “stiff.” 

A typical example of a stiff equation is the equation representing the 
rate of formation of free radicals in a complex chemical reaction. The 
free radicals are created and destroyed so rapidly compared to the time 
scale for the over-all reaction that to a first approximation the rate of 
production is equal to the rate of depletion. This is the notion of the 


pseudo-stationary state. In some cases such as the fast reactions occurring 


in flames or detonations, this approximation is not sufficiently accurate. 
The method described in the present paper provides a means for obtaining 
solutions to equations of this type to any degree of accuracy. 

The numerical procedure described here can easily be extended to sets 
of simultaneous first-order differential equations. In any particular 
region, the differential equations can be uncoupled by introducing suitable 
linear combinations of the original dependent variables. Some of the 
uncoupled equations may be ‘‘stiff’’ in which case they can be integrated 
by the methods discussed here; whereas other uncoupled equations may 
be integrated by the more usual procedures. 

1. Concept of Stiff Equations. Consider the first-order differential 
equation, 


dy 


= [y — G(x)]/a(x, y). 
dx 
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The right-hand side of this equation represents a general function of x and 
y which for each value of x has a root, y = G(x). If Ax is the desired 
resolution of x or the interval which will be used in the numerical integra- 
tion, the equation is ‘‘stiff”’ if 


a(x, y) 


“7% 9 
Ax (2) 


“~ 


and G(x) is well behaved [i.e., varies with x considerably more slowly 
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FIGURE 1 


Slopes, dy/dx, for a Typical Stiff Equation, dy/dx = 5(y — x*) 


than does exp (x/a(x, G(x)))|. The geometrical significance of the stiffness 
is shown in figure 1. Here the family of solutions horn out as one proceeds 
in the positive x direction. This can be chosen as the positive x direction 
without loss of generality. The slope, dy/dx, is drawn at regular intervals 
of x and y. For a particular value of x: for large values of y the slope is 
very large and positive; for small values of y the slope is very large in the 
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negative sense; in the vicinity of y = G(x) the slope changes rapidly from 
large positive to large negative. 

Looking at equation (1), it appears that if a(x, y) is sufficiently small 
there is a special solution, y = ¥(x), which lies close to y = G(x). Toa 
first approximation, ¥ is given by ¥, , 


dG 


ax 


YO = G + a(x, G) (3) 
The second approximation, V, is obtained by substituting V“ for y 
in dy/dx and in a(x, y), 


oe : : dY™ 
VY «= G + a(x, Y™) (4) 
dx 


In this manner, an arbitrary order of approximation of Y may be obtained 
in terms of the lower orders of ap- 
proximation. If a(x, y) does not 
depend upon y, we can write down 
the resulting expression for JY, 

’ — DG 

Y = (5) 

k= 0 Dx 

where D/ Dx is the differential opera- 
tor, a(x)d/dx. Whereas the funce- 
tion V(x) remains close to G(x), 








.— 
FIGURE 2 
every other solution deviates ex- 
ponentially. This can be seen by subtracting the differential equation 
for VY from equation (J). In the region about ¥(«) (assuming that a(x, y) 
varies slowly with y), 

d(y — ¥) 


a(x, ¥(x)) = (y 
dx 


This integrates to give 


y : dx x 
aa = C exp a(x, V(x)) 4) 


where C is a constant of integration. The qualitative appearance of this 
family of solutions is shown in figure 2. 

2. Numerical Solution of Stiff Iquations Using Forward Interpolation. 
It is very difficult to integrate ‘‘stiff’’ equations by ordinary numerical 
methods. Small errors are rapidly magnified if the equations are inte- 
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grated in the direction such that the family of solutions horn out, whereas 
the numerical solutions oscillate violently about Y(x) if the integration is 
carried out in the opposite direction. 

We have discovered a method of numerically integrating ‘‘stiff’’ equa- 


tions which has the desirable property that regardless of the starting 


conditions and regardless of the direction of integration, a solution is 
generated which rapidly converges to V(x). Thus for example, in figure 2, 
if we use this system of integration starting at the point (x9, yo) we obtain 





FIGURE 3 


Integration of dy/dx = 5(y — x?) in the positive x direction from (x = 5, y = 50) 
Here a/Ax 1/5, G(x) = x*, Vix) = 0.08 + O.Ax + x? 


the function indicated by the dotted line rather than the solution to the 
original differential equation which passes through this point and rapidly 
approaches infinity. 

Let Xo, “1, » Xu be a set of values of x spaced a distance Ax be- 
tween successive points. Then the subscript on any quantity indicates 
that it is evaluated at the corresponding value of x. We wish to 
evaluate y, from a knowledge of y at the previous points. Suppose, for 
example, that y(v) can be approximated locally by a straight line passing 
through y, and y,.;. Then this straight line has the slope: 
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dy ie We Ss 
, : (S) 
dx/, Ax 


Now, using this expression for the slope in evaluating the terms of the 
differential equation, equation (1), at the forward point one obtains: 


(‘”) Yu ~ Vat -¥n — Gr 
dx/, Ax A(Yn» Xn) 














eit Sm 


“x1 
FIGURE 4 


Integration of dy/dx = 5(y — x?) in the negative direction from (x = 5, y = 50). Here 
a/Ax = — 1/5, G(x) = x*, ¥(x) = 0.08 + 0.4% + x? 


( A(Xns Vu) 
in = Vn-1 
Ax 
A(Xns Vn) 


Ax 
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If a(x, y) does not depend upon the value of y, equation (10) gives an 
explicit solution for y, in terms of y, 1; otherwise, equation (10) gives 
an implicit relationship between y, and y,;. Thus starting from a point 
(xo, Yo), equation (10) provides a numerical solution to the differential 
equation, equation (1). The curious fact is that quite regardless of the 
starting point, thus numerical solution seeks and approximates V(x) in 
the manner shown by the dotted curve in figure 2. The reason for this 
behavior may be seen geometrically from figure 3. Here yo lies on the 
curve y = So(x) where So(x) is a solution to the original differential equa- 
tion. Equating (dy/dx), from equations (8) and (9) is equivalent to 
asking what solution to equation (1), y = S;(x), has a tangent at the point 
x, which passes through the point (x9, yo). Similarly in passing from x, 
to X, ete. Since the slope for a stiff equation has a reasonable value 
(neither tremendously large in the positive or negative sense) only in the 
vicinity of G(x) or Y(x), it is clear that this numerical integration scheme 
limits us to approximating Y(x). 

It is interesting to notice that this method of integration may be used 
in either direction. An integration in the negative direction so that Ax 
is negative (i.e., in the direction of convergence of the manifold of solutions) 
is illustrated in figure 4. In this case the solution approaches the 
asymptote without oscillations. The integration in the two directions is 
not reversible since the definition of the forward point depends upon the 
direction of integration. This fact may be used numerically to hunt the 
special solution by first integrating in one direction and then reversing 
the direction. 

3. The Asymptotic Form.—If a(x, y) is a function only of x, it is easy to 
see how the numerical procedure leads to an approximation of VY. Since 


G, = Y, — a(x,)(dY/dx),, (11) 
equation (10) can be rewritten in the form 


a(X,) 
Ax 
a(x,) 
Ax 


[Vn - y, 1 + €| 


where 


; ; dY 
—) - + £ 1 + Ax ix 7 (13) 


Expanding ¥ and (d¥/dx) in Taylor series about the point, x,—1, it follows 
that 
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bs iaah 4 vis » (#2) ‘ ue 
= xv)* x aes +) 
, 2 *) 1x? n—l 3 dx’ Xa-l 


Looking at equation (12) it appears that there are two reasons why y, 
differs from Y,. The rate of convergence of y toward Y is measured by 
the factor 


a(x,) 
Ax 
a(Xp) 
i. Ax 


The larger the interval size Ax, the smaller the value of ) and the more 
rapidly y, converges toward Y,. The series converges for values of 4 
less than unity and approaches an asymptote given approximately by 


: E a*y l ay | 
y = Y + a(x) ] = (Ax) + — (Ax)? Pave i (16) 
2 dx? 3 dx’ 
This discrepancy between y and Y is due to the poorness of the fit of Y 
by successive linear segments. The smaller the intervals, Ax, the smaller 
this error. 

The error in the asymptote of the series may be reduced by using a 
three-point formula for (dy/dx), in which y is fit to a quadratic passing 
through y,, ¥,-1 and y,-». Since the slope is evaluated at x,, the same 
sort of forward interpolation is used as in the linear approximation. In 
place of equation (10), one obtains the following relation: 

, A(Xny Vn) 1 d(Xn, Va) 
Gy —2 Vr—l 


as ~ 2 & 


' 3 a(X_, Vu) 
2 Ax 


In this case the difference between y, and Y, is given by! 
a(x,) 
Ax 
3 a(Xx,,) 
2 Ax 


Y, = 


—2(yp Siege Y, 1) + 9 (Ve-2 — Vo 2) 


Max)s(1) L (anys (#Y) F 
3 ” Of ext et - dc’ | «2 


Hence the asymptotic form is approximately 
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v4 eink tad ee a 19 
) (Xx + op hes x)* ae ° (1 
J oe 8) . dx? 1s me dx* ; 


A geometrical argument of why y, approximates Y,, follows the same line 
as for the linear case. 
The use of a cubic to fit y leads to the expression 


' (Xn Vn) 3 
G n +- = rv n—t + Z y n—2~ , ay n~3 
Ax 2 7) 


LL G(Xpy Ve) 
6 Ax 


with a discrepancy between y, and JV, given by’ 


a(x,) 
Ax 
Ll G(x.) 
6 Ax 
; 3 l 
—3(Ya-1 — Vrs) 4+ 5» (Yn-2 — Ve-2) ~ — (Yo — 


” 


| (—) | : (=) 
(Ax)‘4 4 (Ax)5 : os 
| dx’ | 1 20 BOT o=3 


The use of a quartic leads to 


. A(Xny Vn) ‘ 
a . ly, 1 OV 
. Ax 


20 O(Xny Va) 
12 Ax 
with the discrepancy between y, and JV, given by 
- I(y,, ro va i) + 
a(x, ) 
Ax 
25 a(x, ) 
12 Ax 
(23) 


Thus by taking higher order polynomials to fit y at a large number of 
points, it is possible to obtain progressively more excellent fits of VY. In 
all cases, the use of forward interpolation in the evaluation of (dy/dx), 
in’terms of Y,, Vn—1, is the mathematical operation which forces y, to 
seek out and approximate J,,. 
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THE ELEMENTARY BASIC PRINCIPLES OF THE UNIFIED 
THEORY OF RELATIVITY* 
By VAcLAV HLAvaTy 
DEPARTMENT OF MATHEMATICS, INDIANA UNIVERSITY 
Communicated by T. Y. Thomas, January 14, 1952 


1. Introduction. —The unified theory of relativity exposed recently by 
Einstein! is based on three principles: (A) The introduction of a non- 
symmetric tensor g,, in the space X, of the relativity; (B) the introduction 
of a non-symmetric connection ['y, by means of 


0 

ox fen ™ Token + Tighe (1) 
and finally (C) the introduction of a (seemingly) overdetermined system 
of conditions imposed on the ['y,, which yields g,,,. 

In the subsequent sections we will deal with each of these basic principles. 
However, we confine ourselves to results only. The corresponding detailed 
proofs will be given in a subsequent series of three papers in the Journal 
of Rational Mechanics and Analysis. 

2. Principle A.-Denote by hy,(k,,) the symmetric (the skew sym 
metric) part of g,, and by g, h, k the corresponding determinants. Through- 
out this paper we assume h # 0. If m = 4 and if Ay, 1s of the signature 
+-+-+-— then there are in general two sets of bivectors B,, By totally 
perpendicular, which are privileged in the sense that they are polar con 


jugate with respect to the cone /,, as well as with respect to the linear 


complex k,, of bivectors. Closely connected with them are four sets of 
(imaginary) bivectors each of them being polar self conjugate with respect 
to hy, and k,,. 

Projecting this configuration from an arbitrary point P of our space 
X, into the ideal space //; of the tangent space 7(P) of XY, at P, one obtains 


a linear line complex A containing a linear line congruence C (whose axes 
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are the ideal lines of B,(P) and B.(P)). Both line manifolds K and C 
are defined intrinsically, by means of g,, only. 

In another paper’? I have shown that the projective three-dimensional 
spinor space, .S;, associated with the (centered) four-dimensional vector 
space of the (special) theory of relativity is fully characterized by a linear 
line complex and a linear line congruence contained in it. Also these 
manifolds are defined intrinsically. Identifying these manifolds with K 
and C we may identify //; and S;. This identification is not affected by 
the coérdinate transformations (which are different in //; and S;) because 
in both cases the line manifolds are defined intrinsically. Hence the 
introduction of a non-symmetric tensor gy, in X4 leads in the most natural 
way to the theory of the spinor space of the relativity theory. 

3. Principle B.—\In the subsequent sections we denote by} ie: the 


v 


Christoffel symbols of h,,, by //,;, its curvature tensor, by V, the symbol 
jv 


of covariant derivative with respect to- 2 
Aus 


Moreover we put A = 


hv,V, and °A = °h™ where °h™ = h™(P). The curvature tensor 


Ox*dx"” 

of Ty, will be denoted by P,/,. Ifn = 4 then e 

Ricci density (with components +1, —1, 0). 
The equations (1) split in two systems 


_ 
“e“* denotes the four vector 


(a) ry, = , Me } + Sy, — h(S3ken + Soko) 


(b) VA, = 2SLX Sh, (2) 


where X,°%, is a tensor, function of hy, and k,,. Computing S2, from 


(2(b)) and substituting it in (2(a)) one obtains the solution Ty, of (1). A 
necessary condition for (2(b)) to admit exactly one solution is g 4 0. A 
sufficient condition for (2(b)) to admit more than one solution Sy is g = 0. 
However, there are cases where (2(b)) does not admit any solution S., 
(and consequently (1) does not admit any solution T'y,). 

4. Almost Riemannian Space.—Our space X,, will be termed an almost 
Riemannian space R, if ky, = €Ay,, Sy, = € 5x, where €— Ois a constant 
and A,,, sx, are tensors different from zero. In R, we obtain from (2(b))* 


v 


l 
(a) Sion = 5 W(V Kay + VaKan — VuKau) (mod €) (3) 


ap Mu 


and 


(b) ry, = : ms \ + € [sx — W'(sQ. ky, + Sauk,)] (mod e?). (3) 
uM 
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Substituting from (3(a@)) in (3(5)) we obtain the solution ['Y, of (1). 
5. Principle C.-lf n = 4 the conditions 


Sf = 0 Pi ur) = () ee’ p el — ¢ (4) 


wpr 
represent 4 + 6 + 6 = 16 differential equations for 16 unknowns /y and 
kyu. These conditions reduce in R, to 
(a) V’K,, = 0 (mod e), (v* = hv) 
(b) AK, + A (Wy2,K og + 2a x’ Kuis) =O (mod e) 
(c) Veo! + Ze HTT 8 Kg — Kayllyish®) =O (mod e) (5) 


aw 


where J? = e*”°V.K,, is (up to a numerical factor) the Einstein electric 
current vector density. In order to find a solution of (5) we introduce a 
neighborhood N,(P) of a general point P(*):N,(P) is the set of all points 
whose coérdinates x” satisfy the condition (x — x")' <.«. If 


hy, = “hy, + « hy, in V.(P) (6) 


then (5) reduces to 


fa) 
°AK,, =0 (mod e), Pa ird,; =O (mod e). 


It is not difficult to find a solution A,, of this system. Thus for instance 


Ky = A(z — t), Kg = B(x + 1y), Ay, = 0 for Aw ¥ 12, 34 (mod €) 
in V(P) 


is a solution, provided B is an analytic (regular) function of x + ty (and 
hn = Ph» = hs; = —°hw = 1, the remaining °h,, = 0). 
6. The First Einstein Conditions.—If we put 


re) AY aC) eee en Hee 
Ry _ Ox? Lan i revs, e 9 (= Pha) to I io) + ryul fas) 


“ 


then Einstein's conditions imposed on the Ty, are 


SZ=0, R, =0 (7) 


“ 


and represent in this form n + n? conditions for n* unknowns Ay, and k,,. 


: .  (e-4) (n — 1) ‘ = 
However, in R, they splitin mn” 5 ar = n’ conditions namely 


(a) Hi, =0 (mod €’), (b) AKy, + h™ID2,K,,=0 (mod e). 
(8) 


This system admits a solution hy, = o°h,, (mod e*) in N,(P) where a is 
an explicitly defined function of position. Moreover we have in this 
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case /1,), = 0 (mod e) in NP). Hence a codrdinate system may be 
found for which /,, are constant in V,(P) so that (S(b)) reduces to ° AA, = 0 
(mod ¢) in NV,(P). Moreover one obtains the same condition if (6) holds, 
or if one assumes that the set of bivectors B,(B.) generates a non-holonomic 
surface which differs only mod ¢ in V,(P) from a set of holonomic surfaces. 
However, in this case (S(a)) leads to Ay, = °Ay, (mod e*) in N,(P) while 
if (6) holds the equation (S(a)) reduces to 


oO? o” ; °? 


, 


ra) / OLva ] . ] Cite 
Any - (> I t Ox*0x" tem Ox"Ox" 


/ 
hi) (mod €). (9) 


It is not difficult to find a non-trivial solution of (9). In particular if 


(mod €)* 


va 


/ 
A 
x 


then (9) reduces to 

°Ah,, = Q (mod e) (10) 
If °hy = Phe = Phys = —°hw = 1 (and the remaining °/,, = 0), then 
(10) admits a solution 


hy, = Pr, (x + iy, z — t) (mod © in N,(P) 


where ~,, = P,, are arbitrary (regular) functions of x + ty. 
7. The Second Einstein Conditions.Any solution of (8) is also a solu- 


tion of the second Einstein system 


0 
Su = 0 Rin) saa Sy oon =O (11) 


which (for m = 4) represents 4 + 10 + 4 = 1S conditions for 16 unknowns. 


This system reduces in R, to 


(a) 11), =0 (mod e*) 
(b) Al’ + &*™’ (IDLE, Kg, — 4laorV pK ,,)h =0 (mod e). (12) 


Awa awh 


In all cases mentioned in section (6) we obtain the same solution for A,,, 
while (in these cases) the equation (12(b)) reduces to °A/” = 0 (mod e). 


If °Ay = She = Phys = — Phy = 1 the remaining °h,, = 0, then /’=f’(x + 


iy,z — t) (mod e) in V,(P) is a solution of this equation, provided that 
f’ is analytic (regular) in x + vy. 

* Prepared under Army Contract DA-33-008-Ord-224 

Einstein, A., The Meaning of Relativity, Princeton University Press, 1950, 


Appendix II. 
2 Hlavaty, V., “Spinor Space and Line Geometry,”’ Canadian J. Math., 3, 442-459 
(1950), and “Spinor Space and Line Geometry II’ (to be published in the Journal of 
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Rational Mechanics and Analysis) where the Dirac equations are interpreted by means 
of the complex and the congruence meutioned above. 
3A B (mod e”) is the analytical description of the statement that the difference 


, 


A — Bis of the magaitude order of «" , m’ > m 
* It is worth while to stress that (9) and (10) are covariant equations {thanks to (6)} 
In particular (10) cannot be enforced by a proper choice of coérdinates! 


HOMOLOGY RELATIONS ON REGULAR ORIENTABLE 
MANIFOLDS 


By MARSTON MORSE 
INSTITUTE FOR ADVANCED STUDY 
Communicated January 25, 1952 


1. Introduction.-The theorems stated in 1927 by Morse in reference | 
are becoming increasingly significant in questions of homology, homotopy 
and category. Studies related to the unsolved Schoenflies problem by 
Morse and Baiada are based on these theorems. It is for this reason 
that the unpublished proofs of these theorems as written in 1927 will be 
here given with relatively obvious extensions whereby chains over the 
field of integers mod 2 are replaced by chains over any field &. The con- 
nectivities used in reference | will here be replaced by the dimensions 
R, of the respective homology groups, always using cycles over ®. 

The complexes involved will be similar to those used in reference 2. 
They are broken up into cells with the properties of “‘Block-ketten’”’ of 
reference 5, page 78. The ranks of the incidence matrices will be deter- 


mined “relative” to &, replacing the elements a of a matrix by ea, where e 


is the unit in #, and evaluating the subdeterminants of the resultant 
matrix as elements in ®. 

2. Regular Manifolds and Regions..A regular n-manifold >>, in a 
Euclidean (” + r)-space /:,4,,” > 0, r > Ois a closed, compact, connected 
Hausdorff subspace of /:,,, specially representable as follows. Let L’ be 
an arbitrary open set in a Euclidean space /¢, with rectangular coordinates 
u = (My, ..., Uy). Let y = (yy, ..., Ve4r) similarly represent a point in 
En, Let P be an arbitrary point of }°,. Some neighborhood N(P) 
of P relative to }>, shall be the homeomorph of an open set LU’ ¢ /, under a 
mapping of wel’ onto [Wilw), ..., Pay(u)|e /,4, where the functions 
y, are of class C”, m > O, and possess a functional matrix of rank n. By 
a regular subregion of >~, will be meant a closed non-empty subset of }>, 
bounded by a finite set of disjoint regular (7 |)-manifolds. Unless 
otherwise stated all manifolds including boundaries of regions are to be 

wir 


regular and given by representations of class ( Other regular repre- 


sentations of a lower class may sometimes be used. 
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By means of relatively simple homeomorphic approximations it is 
possible to establish our final theorems for manifolds and regions of class 
C’ and even for polyhedral manifolds and regions. We shall not becloud 
the principal problems with such details. Another trivial extension is 
to coordinate manifolds and regions. We have chosen the Euclidean 
case for simplicity of statement of hypotheses. 

Because of local connectedness of our manifolds and regions it 1s im- 
material whether homologies are taken in the sense of the singular theory 
or Vietoris theory. The basic determination of the R, is, however, combi- 
natorial by means of matrices of incidence as in reference 2. 

3. Non-degenerate Functions f..-We shall consider real functions / 
with values f(y) defined for ye >>,. If py: U <-> N(P) is regular and 
of class C’”’ over LU’, f(w) shall be given as of at least class C’’ over LU for 
each Pe >.,. A critical point (uw) of f(y) defines a critical point (uo) of 
f and is termed non-degenerate if the Hessian of f(y) does not vanish at 
(uo). When (uw) is a critical point of f(y), f(y) shall be given as of class 
C’”’ over some sufficiently restricted neighborhood of (wm). Given a 
critical point Py of f and critical value c of f it follows as in reference 4, 
page 173, that there exists a mapping ¢: Ul’ «> N(P») regular and of 
class C’ such that f(¢) — ¢ has the values over l’ of a form quadratic in 
the coordinates of lL’. Such a representation of {| N(Po) will be termed 
canonical. 

By methods exemplified in reference 4, page 243, it is easy to establish 
the existence of non-degenerate functions / over ~ a By the methods 
of reference 3 one can modify the values of f on the boundaries of regular 
regions on which f has no critical point so as to make the gradient of f 
constantly entrant or emergent, and / a constant on each boundary manti- 
fold. 

Let f, and f* be, respectively, the sets )>,,| f(y) S a, and ¥>,,| f(y) = a. 
If aS Blet fag be the set }>,/a 5 f(y) < B. Our ultimate theorems are 
independent of f and concern regular regions C,, their boundaries BC,, 
and their complements on >>,. By suitable choice of non-degenerate 
functions /, in accordance with the existence theorems of the preceding 
paragraph, one can reduce our theorems on regular regions and their 
complements to theorems on fag and fgy where a < min f << 8 < max f < ¥. 
When a given function is modified its class C" will in general be lowered. 
In all cases it will be of at least class C’. 

4. Variation of R,(f,) with a.— This variation is analyzed in reference 
in the case of a regular subregion of a Euclidean space. Reference 
with trivial modifications, thus affords a local analysis of ff when locally 
represented in the form /(y) in terms of the Euclidean coordinates of U’. 
Reference 2 can be simplified by the deletion of sections 10 and 11, since 
these sections become trivial when / is canonically represented. The use 


2 
2 
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of orthogonal trajectories of level manifolds of / naturally takes the differen- 
tial form appropriate to a curved space >°,, see reference 4, page 152. 
Sections 13, 14, 15 of reference 2 can be taken over here with the following 
modifications. We suppose all cells oriented; a corresponding modifica- 
tion of the matrices of incidence is required. In particular the k — | 
matrices of incidence of the (k — 1)-sphere S,_; shall have two columns, 
each consisting of a | above a —1. 

The chains in reference 2 were functions with values over the field of 
integers mod 2. Here the field ® is general and the present connectivity 
R;, defined as the dimension of the 7th homology group, accordingly differs 
from R; in reference 2. An additional difference arises in that even when 
# is the field of integers mod 2, and R,(C,,) is a connectivity in reference 2, 
our present connectivity, using the same field, is R,(C,) — | when? = 1, 

,n. For an n-sphere, for example, we here have Ry = R, = 1, Ri = 
R, = ... = R,-, = 0. Even with these changes in notation equations 
(134) of reference 2, giving the differences AR; taken in the present sense, 
are applicable, provided the ranks p, of the matrices of incidence are 
appropriately taken relative to ®. 

As in reference 2 consider a critical point Py of index k with 0 < k <n. 
We can suppose that the critical value f(P») = ¢ is assumed by no other 
critical point with /-value on the interval (c — «,¢ + ¢€). In reference 2, 
« = e. Asin reference 2, $15, let C, be the complex of unoriented cells 
decomposing f,,, 80 as to include the ‘‘new cell’ (a, ,,-,)and its boundary 
cells, and let C’, be C, with the cell (a,,_,) and the new boundary cells 
deleted. Let the cells of C, and C’, now be oriented and the resulting 
complexes denoted by I, and I’,, respectively. The special product 
cells (a,b,-,), (a’, 160), ete., of reference 2 will thereby be oriented and 


se 


denoted by [a,b,-,], {a’,—~1b’o|, ete. 

Let J, and J’,,1 = 1, , n, be the matrices of incidence of I, and [’,, 
respectively. The matrix J, will have the form of the matrix //, on 
page 383 of reference 2, provided the submatrix //’, is replaced by J’, 
and the last two columns are appropriately modified. These columns are 
to be headed by [a,b’)| and {a,6’)| and their only non-vanishing elements 
will be those in the submatrix 


| a,b’o | | a,b’ | 


(a’,-16'o] | 0 
fa’, —1b'0| —| 0 
(a’y-10'"o] 0 | 
(a’’, 10/9] 0 - | 


Let p,; and p’; be the ranks relative to ® of J, and J’,, respectively. Operat- 
ing on the matrix J, essentially as on //, in reference 2, taking account of 
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, one finds that two cases may occur. In Case 1, p, = p’,, and in Case 2, 


Pre = pe + il. 

Our present matrix J,4;, 7 = 1, 2, ,n — k — 1, similarly has the 
form of //,,; as represented on page 454 of reference 2, provided one 
replaces [1',,; by J’,,; and appropriately modifies the last two columns. 
The last two columns are to be headed by |a,b’,;||a,b’’;| and the four 1’s 
in each of these columns in //,,, here replaced by 1, —1, 1, —1, respectively. 
One similarly obtains J, from //, of reference 2 by replacing //’,, by J’, 


and modifying the last column of /7/,. 
As in reference 2, one infers that 


, 
Pres Pray tl 


and p, = p', + |. The equations at the bottom of page 384 of reference 2 
can then be applied to l,, and T'’,. Hence, as a increases through a critical 
value ¢ assumed by a critical point Py of index k, 0 < k <n, eitherR,(f,) 
increases by | or else R,_\(f,) decreases by 1. The numbers R,(f,) other- 
wise remain constant. 

5. Variation of Rf") with a.—We continue with the critical point 
P,, the critical value ¢ and the constant e«. From the orientability of 
>>, it follows that for a an ordinary value of f the (wn — 1)-manifold f* 
is orientable and hence the Poincaré law R,(f*) = R,-;-1(f") holds. The 
change in R,(f*) asa increases through cis sought. The alternatives depend 
upon the value of the index k and upon n. 

Let BY, be the (n — 1)-complex of oriented cells on the geometric 
boundary of [,. The carrier |Br,| of BY, is f°*s. Although | BI’,| # 
f°‘ it is readily seen that | BI’,,) is the homeomorph of f°‘. We introduce 
the (n | )-complex 


X,~1 = Br, V BY’,, (5.0) 


taking the cells of XY, ; with the orientation and incidence relations they 
have in [',. One obtains Y, -; from BI’, by adding the “new-cells’’ on 
Blayb,~|. Let K, and K’,; be, respectively, the 7th incidence matrices 
of X,_, and BI’, (¢ = 1, ...,m — 1). SupposeeO<k<n. Forz = 1, 
2, ...,k — 1, Kh, = KK’, The matrix A, has the form of the matrix //, on 
page 383 of reference 2, provided A’, replaces //’, therein, and the last two 
columns of //, are modified asin section 4. The matrix A,,;, 7 = 1, 2, ..., 
n — k — 1, has the form of the matrix //,,,; at the top of page 384, provided 
K',,, replaces /7',,, in H/,,; and the last two columns of //,,; are modified 
as in section 4. The equations at the bottom of page 384 then hold for 
i = 1,2, ..., — 1, where a; and p,; refer to A,, and a’; and p’; refer 
to A’, Equations (134) of reference 2 are here replaced by the equations, 


6 ee AE acc 
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R(X 1) RABY’ ,) (a; - a’;) — (py - p’;) ~ ies P's) (5-2) 


n 


, 


where po = po = ps = Pn 0. Recalling that BI’, | is the homeomorph 
of f°~* we have the result 

(a) When O < k <n the connectivities of X,—, and f° 
for the mutually exclusive changes 


€ 


are equal except 


RX. 1) = R.(f' °) + l Case la, t 
Ryear(Xeu) = Real fe) — 1 Case 2a, 


and except for an additional change, independent of k, 
Ryi(Xan) = Raalf*) + 1. (5.3) 


The complex X,,; mediates between BI’, and BI,, or if one pleases 
between f¢~' and f‘t*. In this sense (4) completes (a). The change 
AR, 1 % 1 due to (5.3) is to be compounded with the change in R,—) 
under Case la or 2a when k = n 1, making a resultant AR, ; = 2 or | 
in Case la or Case 2a, respectively. Similarly with (5.4) in (4). 


(b) When 0 < k <n the connectivities of X,, and f°** are equal except 


for the mutually exclusive changes 


Ryegate Reel) l Case Ib 
Ratha sa) 2 R,-( f°") + 1 Case 2b 


and except for an additional change, independent of k, 
Re Xn 1) = RK, i(f' r*) — il. (5.4) 


Observe that Y,_; can be obtained from BI, by adding the ‘old cells” 
on Bla,b,-,| not already in BP,. Let A’’, be the 7th matrix of incidence 
of BY,, 7 = 1, 2, ,n—1. Then K, = K”, fori . oO: n—-k—1. 
The matrix A,,, has the form of the matrix on page 383, provided /1’, 
therein is replaced by A’’,, and provided the non-zero elements in the 
last two columns are those of the submatrix 


| ‘? 
a’ by —x| 


ros ] 
a On k} 


(a'ob’ ,—x—1 | l 
fa’ oD «4-31 —| 
[a’ob’ » oa | 0 
Cat ey 0 


The matrix K,_-.4;, 7 = 1, 2, ,& — 1, similarly has the form of the 
matrix at the top of page 384 provided A’’,,,,, replaces //',,, therein, the 
last two columns are headed by |[a’,b,-,| and |a’’,b,-,|, the last six rows 
headed by [a’jb’, ,-1] [a’jb’’n-« «|, ete., while the four |’s in each of the 
last two columns are replaced by 1, —1, 1, —1, respectively. If p,; and 
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p', are, respectively, the ranks of X,_, and BI, relative to , then either 
Pn-k = p'n—-» OF Pr-k = p'»-+ + 1 while p, —k+j = p's sty + i, j = I, 2, 
...,k — 1. If a; and a’, in the equations at the bottom of page 384 
then refer to X,_; and BI’,, respectively, these equations hold for 1 = 1, 
2,...,m — 1. Statement (5) follows. 

Statements (a) and (b) imply the following theorem. 

THeoreM 5.1. Jf a value c of f is assumed at just one non-degenerate 
critical point of index k then the differences AB, = R,(f°t*) — Ri f°~*) for 
a sufficiently small ¢ > 0 have the following values: 


AB, == AB,, 1 = (k = 0) (5.5) 
AB, AB, = (k = n) (5.6) 


Either AB, = AB, _.., = k= 1,2,...,.8— 1 (3.7 
or AB,_, = AB,-, ttn iiaes il’ 


) 


Either AB, = AB, = 2 ‘ ‘ 
a cag ill (Qk=n 1) (5.8) 


Either AB, = AB, 1 
or AB, AB, 
or AB, = AB, 


Either AB, = AB.» = : ‘ (5.10) 
or SB, = ABRw = -: i 


All Other Differences AB, = 0: 

In proving Theorem 5.1 one passes from f** to X,-, to f*t* 
using (a) and (6). In deriving the values of AB, resulting from 
the changes implied by (a) and (6), equations (5.3) and (5.4) can 
be neglected since their resultant effect on AB,_; is zero. Apart 
from equations (5.3) and (5.4), four algebraically different indices appear 
in (a) and (6). The conditions that two of these indices have equal 
values give rise to side conditions in Theorem 5.1. In particular the 
equality k = n — k — 1 implies 2k = n — 1,k — 1 = n — k implies 
2k = n + 1, and the two remaining equalities, k = n — kandk — 1 = 
n — k — 1, both imply 2k = n. Let two indices 7 and 7 on the interval 
[0, nm — 1] be termed dual if i + 7 = n — 1. The conditions 2k = n — 1, 
and 2k = n + 1 imply that k and m — k are self dual, respectively, while 
2k = nimplies that k and k — | are dual. The two indices appearing in 
each equality in (5.1) are dual. According to Poincaré, 4B, = AB, if 
i and j are dual. 

The cases k = 0 or m are trivial. When 0 < k < n Cases la or 2a may 
be combined with Cases 1b or 2b giving four a priori combinations. Apart 
from the special case 2k = n, the Poincaré law, however, implies that 
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Cases la and |b occur together and Cases 2a and 2b likewise, as we shall 

Consider equations (5.7). In Case la, AR, = | as one changes from 
f°‘ to X,_4, neglecting (5.3) when k = nm — 1. Under the side condition 
2k # n — 1, nof (5.7) one hask # n — k — 1, n — k, so that on passing 
from X,_; to f*~* in accordance with (6) AR, = 0, neglecting (5.4) when 
k =n — 1. Thus (a) and (}) imply a resultant change, AB, = 1, on 
passing from f‘~‘ to f°**; by virtue of the Poincaré law one has a resultant 
AB, «1 = AB, = 1. As a consequence, Case 1b must hold with Case 
la. Similarly, in Case 2a, AB,_, = AB,_, = —1 on passing from f*~* to 
f¢**, so that Case 2b must hold. 

The remaining equations similarly follow. The first alternative in 
(5.8), (5.9), (5.10) occurs when Cases la and |b both occur, the second 
when Cases 2a and 2b both occur. When 2k = n the third alternative in 
(5.9) occurs when Cases la and 2b both occur, or Cases 2a and Ib. 

6. Critical Points of Increasing, Decreasing and Neutral Types. Let 
Z, be any one of the subsets f,, f*, Comp (fa), faa fas. AS a increases 
through a value c assumed at just one non-degenerate critical point Po 
the changes AR,(Z,) are of one sign or zero, independent of 1. If each such 
difference vanishes P, will be said to be of neutral type rel Z,. If at least 
one such difference is positive (negative) P, will be said to be of increasing 
(decreasing) type rel Z,. 

THEOREM 6.1. Jf Py is of decreasing type rel f,, it will be of decreasing 
type rel f°. 

Let k be the index of P». 

Case 1: k = Oorn. When k = 0, Po is always of increasing type rel 
fa, and when k = n, always of decreasing type rel f*. The theorem follows. 

Case I]: 2k A n,0O<k <n. As has been seen in the analysis of the 
matrix of incidence J, of I’, in section 4, Py is of decreasing type rel f, 
if and only if p, = p’, + 1, where p, and p’, are, respectively, the ranks of 
J, and J’,, or equivalently if the cycle boundary 8{a;,b’o] of [a,b’9] does 
not bound on I’,. Turning to f* recall that when 2k # n, Py is of de- 
creasing type rel f* if and only if Case 2a occurs in (a) in section 5, or 
equivalently if p, = p’, + 1 where p, and p’, are the ranks of A, and K’;, 
respectively, or equivalently if and only if B{a,b’o| not ~ 0 on BI",. 
The latter condition is implied by 6{a,b’5| not ~ 0Oon [’,. The theorem 
follows in Case IT. 

Case [I]: 2k =n, 0O<k <n. By virtue of Case II the theorem holds 
if m is odd. We take advantage of this fact when nm = 2k by replacing 
>, by a manifold }°,,: of odd dimension and f by a function F appro- 
priately defined on >, 41. 

The Definition of 3>,,: and F: We are supposing that >>, lies in a space 
E,4, of points y with coordinates (y,, ..., Yn4,). Without loss of generality 
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one can suppose that y;| 0, # 0. Let /,,,, be regarded as the subspace of 
the space /,,,,, of points y’ and coordinates (yo, Yi, -.., Ya¢r) On whieh 
yo = 0. If A is any subset of /¢,,, on which y,; > 0 let 0A be a subset of 
Fnaryi generated by revolving A around /,,, as follows. For y = (\, 


» Yn¢r) € A let Oy be the circle of points in /,4,4; with coordinates 


(v1. cos 0, yy Sin 0, Yo, ...5 Va+e) (O'S 0 << 2a); (6.1a) 


and let 0A be the union of the circles Oy for ye A. In particular let 
> EP 0>>.- One sees that » wer is regular with ee We now define 
F over 141 by setting 


Fy, cos 6, y; sin 6, (v2, ..-5 Vn¢r) = f(y) + ncosé@  (6.1b) 


for every point y (Sis... 55 3-Meaeey e : where 7 > 0 is a constant yet to 


be chosen. 

Observe that F has a critical point on the circle Oy if and only if y 1s 
a critical point of /, and if @ = Oor am. If y has the index k as a critical 
point of f then a point on the circle Oy has the index k + 1, or k, as a critical 


point of F, according as 6 = O or 7. 

Proof of Theorem 6.1 in Case [11: Suppose for simplicity that the given 
critical point Py of f has the f-value 0. Let e be so small a positive constant 
that |—e, e| is an interval of ordinary values of /, 0 excepted. Let F° = 
F\ (yn = 0). The set equalities are clear: 

yr, = Of. F, = Of.. (6.1¢) 

For e fixed an obvious use of the orthogonal trajectories of / shows that 
when 7 > 0 is sufficiently small the relations (6.le) are replaced by the 
topological equivalences: 


F_. i. = Of_. F,.<— F, = Of, (6.1d) 


including the equivalences: 
F< Of, Fee—-ofe. (6.1e) 
If Py is a critical point of / of decreasing type rel f,, then 
- Ralf.) = —1. (6.2) 
In accordance with well-known formulas for product spaces, (6.1d) and 
(6.2) imply that 
Ryal( Fe) — Rea FP.) l, R.(F.) — RF.) = —-1. (6.3) 
This is possible only if the two eritical points P; and P, of FP, with /#-values 
on the interval (—e, e), are both of decreasing type rel /,. Theorem 6.1 


has already been seen to be valid for /, as defined over jv when the 
dimension of }°,4: is odd. In the case at hand n + 1 = 2k + 1. By 
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virtue then of Theorem 6.1, P; and P2, as critical points of F, are of de- 
creasing type rel /”. Hence by definition 


RAF.) Ss RiF-*) (¢- a }...., a) (6.4) 


where the sign < prevails for some 7, Consistent with the product homol- 
ogy theory and (6.le), with (6.4), and with Theorem 6.1, one can only 
have 


RCT) Shas (¢ =O, ...,% 1), (6.5) 


where the sign < prevails for some values of 7 (in particular for 7 = k and 
k — | according to (5.9)). Thus Py is of decreasing type rel f* and the 
theorem holds in Case III. 

THEOREM 6.2. Jf the critical point Py is of decreasing type rel fq it is 
of decreasing type rel Comp fq. 

If one sets F = —/f over = Py is a critical point of F of index n — k. 
Moreover for each ordinary value a of f, R)(F_.) = R,(Comp f,), 7 = 0, 
1, ..., m. An equivalent way of stating Theorem 6.2 is accordingly as 
follows. If the critical point Py has an index k, and if for every sufficiently 
small « > 0 


R, thee.) sa R, (fe ) = —l, (6.6) 
then for every sufliciently small « > 0 
R,-x(Comp fey.) — Ra+(Comp f-_.) = —1. (6.7) 


To prove Theorem 6.2 suppose Py of decreasing type rel f,. Then Po is 


of decreasing type rel f° by Theorem 6.1, hence of increasing type rel F“ 


since F = —f, of increasing type rel Ff, by Theorem 6.1, and hence of 
decreasing type rel Comp f,. This establishes the theorem. 

THEOREM 6.3. If a < f(P») and if the critical point Py is of decreasing 
type rel faa, Po will be of decreasing type rel f*. 

The proof of Theorem 6.3 is at most trivially different from that of 
Theorem 6.1. Theorem 6.4 follows from Theorem 6.3 essentially as 
Theorem 6.2 follows from Theorem 6.1. 

THEOREM 6.4.) [f a< f( Po) < y and if the critical point Po ts of decreasing 
type rel faa, Po will be of decreasing type rel f,,. 

7. Complimentary Regions and Boundary Manifolds. Throughout this 
section the following notation will be used. The function f defined over 
a n > 1, shall be non-degenerate with distinct critical values. The 
constants a << 8 < y, min f < 6 < max f, shall equal no critical value of /f. 
Otherwise a, 8, y shall be fixed but arbitrary. The parameter a shall 
increase from ato 6. Let \/, be the number of critical points of f of index 
i with values on the interval (a, 8). Of these \/, points let P; and Q, be, 
respectively, the numbers of points of increasing and decreasing type 





256 MATHEMATICS: M. MORSE Proc. N. A. S. 


rel faa Q’,-, and P’,_, these numbers rel f,,, and let p,, qi, 7; be, respec- 
tively, the numbers of these \/, points of increasing, decreasing and 
neutral types, rel f’. Observe that 


| U Ser = San has n Se, = Te, (7.1) 

Set d; = q; — Q;. These differences are non-negative by Theorem 6.3. 

THEOREM 7.1. The algebraic increments AR, and AB, in the numbers 
Ri faa) and Ri f*) respectively, as a increases from a to 8 satisfy the relations 
(4 = O, 1, »n— 1), 
A(R, + R,, i-1 ™ B,) = 

d, + diys + d, i—l +d, iT vj + Tn i-l = 
In accordance with section 4 


M,=P,+Q, @R:° P:-Qu €=0,... 


where Qo = Q,41 = P, = 0. Using the results of Theorem 5.1 


M, =pitqatn (a4 = ©, 1,...., ”) hse 
AB, = pi — Gi4t + Pa-i-t — Gn-i (¢ = 0,1,...,#8—1),$ 


(7.4) 


where p, = go = Oand r,; = 0 for 21 An. From (7.3) and (7.4), respec- 
tively 

M, 2 M,, inf = A(R, + R,, i—1) + OF a Qin = Q, i—l + 0, i 

M, + M, i-! AB, + gi + Gist + In i 1+ qn CONE Pas 1 


and (7.2) follows on equating these two representations of MJ, + A/,,_; 1. 
Corouvary 7.1. If C, is an arbitrary regular subregion of >, 


Ri(C,) + Rr-irlC,) — R(BC,) = 0 (@=0,1,...,2—1). (7.5) 


Choose f (as is possible) so that fg = C, and choose a < min f. The 
left member of (7.2) becomes the left member of (7.5). 

COROLLARY 7.2. Necessary and sufficient conditions that >, have the 
connectivities of the n-sphere are that (a) for any non-degenerate { with distinct 
critical values each critical point with f-value < max f be of the same type 
(increasing or decreasing) rel f° as rel f,, or (b) for each regular subregion 
Gof oe the left members of (7.5) are zero. 

Choose a < min f and 8 greater than each critical value less than max /. 
Then f* is a topological (n — 1)-sphere. If }>, has the connectivities 
of the n-sphere the left members of (7.2) vanish. Hence each d, and r, is 
zero. The necessity of (a) follows. Conversely if (a) holds the right 
members of (7.2) are zero and hence the left. Choosing a and 8 as pre- 
viously one sees that fg has the connectivities of an n-cell and hence >_,, 
those of an m-sphere. Corollary 7.2 (6) follows trivially from Corollary 
7.2 (a). 
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THEOREM 7.2. Jf AR, and AR’, are, respectively, the algebraic changes in 

Ri( faa) and R,(fq,) as a increases from a to B, then fori =0,1,...,n—1, 
A(R’, i-1 ~ R;) = (Q; ada *, i) + (Qi4s —_ PP’. i-1) < 0. (7.6) 


The parentheses in the right members of (7.6) are non-positive by 
Theorem 6.4. The equalities in (7.6) follow from the relations 


Po+Q = Pon + On: CoG’, ..gf) 
AR; = P; — Qis, AR, = Q'igni — P’; (i = 0,1,...,2 —1) 


Corouiary 7.3. If C, is an arbitrary regular subregion of >, then for 
-2 O01, ..."n—-1 


R,_:-1(Comp C,) — R(C,) S Ra-is(don). (7.7) 
The left member of (7.6) has the form 


Re i i(fa,) _ Ri(f ap) _ R,, i (fay) + Rf"). (7.8) 


To obtain (7.7) from (7.6) choose f so that fg = C,. Take a < min f, 
max f < y. Then each R,(f*) = O and the remaining terms in (7.8) 
become those of (7.7). 

Corouvary 7.4. If ¥, is a regular subregion of >, and if C, and C’, 
are complimentary regular subregions of >, then fori =0,...,n —1 


R,-i(C’n) — RAC.) S Rain.) — Ril[B(C,) — (C, © C’,)]. (7.9) 


This is an obvious consequence of (7.6), using (7.8), on choosing f and 
a, B, y (as is possible) so that 


fe BH) -CGat), f-G.9t. fa~ Gh Liwe. 
(7.10) 

With >>, as in Corollary 7.2 (a), a critical point Py of f with f(Py) < 
max f is of the same type (increasing or decreasing) rel f, as rel f*, and 
hence, if f(P») > min f, of the same type (increasing or decreasing) rel F, as 
rel F*. A review of the proof of Theorems 6.2 and 6.4 then gives the 
following. 

TueoreM 7.3. If >, has the connectivities of the n-sphere a critical point 
of f of index k with O < k <n 1s of the same type (increasing or decreasing) 
rel faq as rel fa,, and (7.6), (7.7) and (7.9) become equalities fori = 1 ..., 
n— I. 

In particular if }>, has the connectivities of the n-sphere the equality 
(7.7) gives a combinatorial extension (1927) of Alexander's result 


R,-:.(Comp C,) = RA¢(C,) + 6-1 »,n— 1). (7.11) 


' Morse, M., ‘The Analysis and Analysis Situs of Regular n-Spreads in (n + r) 
Space,”” Proc. Nati. Acap., 13, 813-817 (1927). 
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2? Morse, M., “Relations between the Critical Points of a Real Function of n Inde- 
pendent Variables,” Trans. Am. Math. Soc., 27, 345-396 (1925). 

§’ Morse, M., and Van Schaack, G. B., “The Critical Point Theory under General 
Boundary Conditions,” Jhid., 35, 045-571 (1934) 

* Morse, M., ‘‘The Calculus of Variations in the Large,’’ Am. Math. Soc. Coll. Publ., 
New York, 1934. 

® Seifert-Threlfall, Lehrbuch der Topologie, Teubner, Leipzig and Berlin, 1934 


CRITERIA FOR LINEAR EQUIVALENCE 
By ANDRE WEIL 
UNIVERSITY OF CHICAGO 
Communicated by Oscar Zariski, January 5, 1952 
1. Equivalence criteria play an important part in the theory of algebraic 
surfaces; such criteria have been established by Severi! for surfaces with- 
out multiple points over the field of complex numbers. My work on 
abelian varieties? has enabled me to extend Severi’s methods and results 


to arbitrary normal varieties in projective spaces. As publication of the 
complete proofs, which I have already written up, may be somewhat 
delayed, I shall give here a precise statement of the main results which 


have been so obtained. 

Once for all, V will denote a variety of dimension n without multiple 
subvarieties of dimension » — 1, embedded in a projective space P of 
dimension .V and not contained in any hyperplane. Let ky be the smallest 
field of definition for V; all the fields to be mentioned will be assumed 
to contain ky. If // is a hyperplane in P, generic over ko, then’ the cycle 
V.// is a variety of dimension n — | without multiple subvarieties of 
dimension n — 2. If Lisa linear subvariety of P of dimension NV — n + 1, 
generic over ko, then V.Z is a curve C without multiple points which does 
not go through any multiple point of V; such a curve C = V.L will be 
called a general curve on V; if, moreover, L is generic over a field k con- 
taining ky, then the general curve C will be said to be generic over k. 

2. Linear equivalence of divisors on |’ is defined as usual, and will be 
denoted by ~. Let W be an abstract variety; let Z be a divisor on 
 & W; if Mis a simple point on IW such that the cycle Z.(1" & MM) is 
defined, we write Z.(V & AJ) = Z(M) & M. Any divisor of the form 
Z(M) — Z(N), where Z is as just stated and MW, N are such that Z(1/), 
Z(N) are defined, is said to be algebraically equivalent to 0; one of my 
lemmas states that such divisors form a group and that every such divisor 
can be written as Z(.\/J) — Z(.N) by means of a divisor Z on the product 
VX Tof V and of a curve IP. 
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Criteria for linear equivalence can be of two kinds; those of the first 
kind state that a divisor X is linearly equivalent to 0 provided its trace on 
a subvariety of |’ which depends upon X is linearly equivalent to 0 on 
that subvariety; those of the second kind state that a divisor XY is linearly 
equivalent to 0 provided it is algebraically equivalent to 0 and its trace 
on a fixed subvariety of V is equivalent to 0 on that subvariety. 

3. The criteria of the first kind are as follows: 

(A) Let X be a divisor on V, rational over a field k; put W = V.H, 
where II is a hyperplane in P, generic over k. Then, ifn 2 3,X.W~O 
on W implies X ~ O on V; the same is true for n = 2 if one assumes at the 
same time that there is an integer m # O such that mX ~ 0 on V. 

For n = 2, X.W~ 0 need not imply Y ~ 0 because of the exceptional 
surfaces occurring in the theorem of Kronecker-Castelnuovo; what is 
more, there are serious difficulties in the way of extending this theorem to 
arbitrary fields.‘ Nevertheless, we have the following general result: 

(B) There is a finite set of divisors D, on V, all algebraic over ky, with the 
following properties: (1) no linear combination of the D;, other than 0, 1s 
algebraically equivalent to 0; (it) let X bea divisor on V, rational over a field 
k; let C bea general curve on V, generic over k; then X.C ~ 0 0n C if and 


only if X is linearly equivalent on V to some linear combination of the D,. 
1. Let f be a mapping of into an abelian variety A; if X is a divisor 
on V, and C a general curve on |’ such that Y.C is defined, I define f(X) 


by putting f(Y) = S{f(X.C)], where S is as defined, loc. cit.,? pp. 28-29; 
this is independent of the choice of C; if k is a field of definition for V, 
A and f, and if X is rational over k, f(Y) is rational over k. If Z(AM) is as 
defined above in paragraph 2, there is a mapping g of W into A such that 
g(M) = f(Z(\M)) whenever Z(.\/) is defined. 

The criteria of the second kind can now be formulated: 

(C) Let X be algebraically equivalent to 0 on V. Then, if f(X) = 0 
for every mapping f of V into an abelian variety, there is an integer m # O 
such that mX ~ 0 on V. 

(D) Let W be an abstract variety; let Z be a divisor on V K W. Assume 
that, for every mapping f of V into an abelian variety, the mapping g of W 
into the same variety defined by g(M) = f(Z(M)) is constant. Then all the 
divisors of the form Z(M) are linearly equivalent to one another on V. 

(i) Let C be a general curve on V. Then, if a divisor X 1s algebraically 
equivalent to 0 en V and X .C ts defined and ~ 0 on C, X is linearly equiva 
lent to O on V. 

Only one consequence of the above results will be mentioned here: 
If k is any field of definition for V, and a divisor Xy is rational over k, 
then every divisor X satisfying a relation mY ~ Xy, where m is an integer 
other than 0, is linearly equivalent to some divisor Y, which is algebraic 
over R, 
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On the other hand, it will be clear to everyone who is familiar with this 
subject that the above results, and particularly (/), open the way to a 
theory of the so-called Picard varieties of arbitrary varieties.® 


! Severi, F., Ann. di Mat. (111), 12, 55-79 (1905); Rend. Pal., 21, 257-282 (1906); 
Rend. Acc. Lincei (V), 30, 328-332 (1° sem. 1921). Cf. also O. Zariski, Algebraic 
Surfaces, Chelsea, New York, 1948, pp. 88-89 and pp. 126-128. 

2 Weil, A., Variétés abéliennes et courbes algébriques, Act. Sc. et Ind. no. 1064, Hermann 
et Cie, Paris, 1948 

* Matsusaka, T., Kyoto Math. Mem., 26, 51-62 (1950). 

‘I am, however, informed by P. Samuel that he has been able to avoid all such diffi- 
culties by the use of a suitable birational transformation. 

® Weil, A., Colloque d’Algébre et Théorie des Nombres, Centre Nat. de la Rech. Scient., 
Paris, 1950, pp. 125-127. I have been informed by T. Matsusaka, and also by A 
Néron and P. Samuel, that they have independently developed theories of the Picard 
varieties. By means of the device referred to in footnote 4, Néron and Samuel have 
been able to derive the basic results on the existence of the Picard variety from the most 
elementary form of the criterion of the first kind. 


RAPID EFFECTS UPON THE RENAL CIRCULATION PRODUCED 
BY NEPHROTOXIC GLOBULIN ADMINISTRATION 
IN THE RAT 


By RicHARD W. LipPpMAN, HELEN U. MartI AND E. ELMo JAcoss 


INSTITUTE FOR MEDICAL RESEARCH, CEDARS OF LEBANON HOSPITAL, 
Los ANGELES, CALIFORNIA* 


Communicated by Linus Pauling, January 28, 1952 


Although rabbit anti-rat-kidmey serum has been shown to contain 
antibodies that localize in the glomeruli,' renal tubule cells also are damaged 
by the administration of this serum or of rabbit anti-rat-kidney gamma 
globulin (nephrotoxic globulin, NTG). Indeed, the damaging effect 
of NTG on renal tubule cells has been shown in tissue culture explants, 
where indirect circulatory effects are excluded.’ It has been suggested 
that the tubular damage produced in vivo by administration of NTG might 
result from either a direct effect upon the tubule cells or circulatory 
disturbances subsequent to NTG administration, or from a combination 
of these two factors.’ 

Various methods have been used to study circulatory changes in the 
kidney. The injection of foreign materials, such as neoprene or India 
ink, may produce distortions of the vascular pattern, as a result of mechan- 
ical or pharmacologic effects. Such a quasi-physiologic method as the 
measurement of renal blood flow by means of para-aminohippurate or 
diodrast clearances cannot be used after NTG administration, since the 
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animals are anuric for at least 2 hours, a period of critical importance. 
The fluorescent dye vasoflavine presents technical difficulties which we 
have not mastered with entire satisfaction. 

For these reasons, we have now devised a very simple method of examin- 
ing the renal circulation. The rat was placed under light ether anesthesia, 
the abdomen was opened widely, the abdominal viscera were retracted 
gently with pads of cotton dipped in warm 0.850) sodium chloride solution, 
and a kidney was exposed to direct observation. The solution under 
investigation was then administered by intravenous injection, in a foot 
vein, and the kidney was observed for changes in the surface appearance. 
At any desired time the renal pedicle was rapidly and tightly occluded 
by application of a clamp, thus trapping blood in the renal vascular bed, 
since the artery and vein were simultaneously occluded. The kidney, 
with clamp on pedicle, was rapidly removed en bloc and a ligature was 
placed around the pedicle before removing the clamp. The entire kidney 
was then dropped into 10°, neutral formalin in 0.85%) sodium chloride 
solution. After 24 hours the kidney was sliced in half, and fixation was 
continued for another 24-hour period. Paraffin sections were made at a 
thickness of 15 yw, and were stained by the benzidine method of Ralph.‘ 

Control globulin (GC) and NTG were prepared by the methods pre- 
viously described in detail.* ° In the control animals 1.0 ml. of GC (17.0 
mg. total protein) was administered by intravenous injection. In the 
other animals 1.0 ml. of NTG (23.0 mg. total protein) was administered. 
For this experiment 39 rats of the Slonaker-Addis strain® were used. 

In the 8 male, 150-g. animals given GC no observable change in the 
surface appearance of the kidney could be seen during 2 minutes after 
the globulin injection. The kidney was clamped and excised, after the 
2-minute period of observation, and sections were made for histologic 
examination. It was noted that the glomeruli contained erythrocytes 
in the capillary channels and, except for occasional rouleaux, their number 
was such that they could be identified as discrete bodies. Small numbers 
of erythrocytes could be seen in the peritubular capillaries and venules. 
Occasional dilated venous channels could be seen in the cortical area, and, 
adjacent to these, the peritubular capillaries appeared to contain more 
blood than vessels in other areas. Blood could also be seen in the major 
vessels and in the vasa recta, extending toward the pelvic papillae. 

After the administration of NTG a dramatic succession of events was 
observed on the surface of the kidney, 1m situ. Within 15 to 45 seconds 
after the injection was given, small livid areas, similar to petechiae, were 
seen on the exposed renal surface. ‘These were so deeply red in color that, 
at times, the intervening renal tissue appeared to be blanched. This 


appearance was only by comparison with the deeply colored areas, and no 
true blanching occurred. The ‘“‘petechiae’’ rapidly increased in number 
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until there was a gradual fusion of these livid areas, and the entire kidney 
became a deep, meaty red, but this color was not cyanotic. The complete 
change occurred within 2 minutes in most instances, with only an occasional 
animal that required 3 minutes. This experiment was performed on 19 
male, 150-g. animals and, of these, in 7 the kidney was clamped and ex- 
cised at intermediate stages, while in 12 the kidney was removed after 2 
minutes or slightly more, when the gross surface change had become maxi- 
mal. When sections were examined, it was found that the glomeruli 


contained many more erythrocytes than were contained in the glomeruli 
of animals which had received GC. Many more rouleaux were present, 
and the erythrocytes were at times so closely packed that it was impossible 


to identify them discretely. The number of dilated venous channels 
in the cortical area was increased, and there was an increase in the amount 
of blood contained in the peritubular capillaries and venules. In addition 
the vasa recta contained more blood than in the control specimens. 

For comparison, 6 specimens were obtained in which the renal vein 
was isolated and ligated for 2 minutes prior to application of the pedicle 
clamp and removal of the kidney. In addition, in 6 other animals the renal 
artery was ligated for 2 minutes prior to application of the clamp and 
removal of the kidney. In specimens in which the renal artery had been 
ligated the entire kidney appeared to be almost devoid of blood, with a 
few isolated erythrocytes in each glomerulus and with very little blood 
in the peritubular capillaries and other vessels. In specimens obtained 
from animals in which the renal vein had been ligated there was massive 
congestion of the venous channels, including the peritubular vessels, 
but there was no significant change in the glomerular blood content. 

While the technique described here does not exclude the possibility 
of intrarenal shifts in blood distribution during the period between appli- 
cation of the clamp and penetration of the fixative, the likelihood of a 
distorted picture is very much less than in most other methods that have 
been used to study the renal circulation. It is of interest to note that, 
in any given kidney section, there was some variation in the blood content 
between different radial segments. The experimental differences, however, 
far exceeded the variability observed in individual sections. 

The most remarkable observation in this experiment was the speed with 
which NTG admimstration affects the renal circulation. Since the mixing 
time in the rat is approximately 2 minutes,’ and the vascular changes ob- 
served in the kidney begin after only 15 seconds, it would seem probable 
that such changes occur almost instantaneously when NTG arrives at 
the kidney. In previous work the speed with which such major changes 
occur has not been appreciated, and observations have usually been made 
at much longer intervals, ranging from 80 minutes to | week after the 
admunistration of NTG, 
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Obstruction to the venous outflow from the kidney, as demonstrated 
here, fills the peritubular vessels but does not appreciably alter the glomer 
ular blood content. After NTG administration, the glomeruli, as well 
as the peritubular and venous channels, contain an increased amount 
of blood. Since efferent arteriolar constriction would increase the glomer- 
ular blood content but probably would diminish the amount of blood 
in venous channels, it would seem reasonable to suppose that the events 
observed are explained by an increased amount of blood reaching the 
kidney from the arterial side, rather than by obstruction to the venous 
outflow at any point. This idea is confirmed by the color change which, 
when grossly observed, does not have the appearance of cyanosis, which 
one would expect with venous congestion alone. 

Exploratory experiments performed with India ink injection and vaso- 
flavine, in spite of the technical difficulties, have given results that confirm 
the information obtained by the much simpler method described here. 

In the light of these experiments, it seems clear that the administration 
of NTG does not result in tubular ischemia. Additional preliminary 
experiments have shown that the renal hyperemia observed here lasts 
for some time, certainly more than 15 minutes, although possibly less 
than | hour. It would seem more reasonable to believe that the tubular 
damage after NTG administration is the consequence of a direct, non- 
circulatory effect upon the tubular cells. 


* This work was aided by a grant from the National Heart Institute, of the National 
Institutes of Health (H 132 C3), and by a grant from Ciba Pharmaceutical Products, 
Inc., Summit, N. J. The authors are grateful for the technical assistance of Ruth 
Taniguchi. Dr. Lippman is a Fellow of the John Simon Guggenheim Memorial Foun- 
dation. 

! Pressman, D., Hill, R. F., and Foote, R. F., Science, 109, 65 (1949). 

? Lippman, R. W., Cameron, G., and Campbell, D. H., Proc. Nati. Acan. Scr., 36, 
576 (1950). 

3 Lippman, R. W., Marti, H. U., and Campbell, D. H., Arch. Path. 53, 1 (1952). 

4 Ralph, P. H., Stain Technol., 16, 105 (1941). 

5 We are deeply grateful to Prof. Dan H. Campbell and Jay Banovitz, Gates and Crel- 
lin Laboratories of Chemistry, California Institute of Technology, for preparation of 
the NTG and GC used in these experiments. 

6 Addis, T., and Gray, H., Growth, 14, 49 (1950). 

7 Lippman, R. W., Proc. Soc. Exp. Biol. Med., 66, 188 (1947). 





264 ZOOLOGY: WEISS AND GARBER Proc. N. A. S. 


SHAPE AND MOVEMENT OF MESENCHYME CELLS AS FUNC- 
TIONS OF THE PHYSICAL STRUCTURE OF THE MEDIUM. 
CONTRIBUTIONS TO A QUANTITATIVE MORPHOLOG \* 


By PAuL WEISS AND BEATRICE GARBER 
DEPARTMENT OF ZOOLOGY, UNIVERSITY OF CHICAGO 
Communicated January 29, 1952" 


Introduction. Historically, the description and classification of cells 
has proceeded trom visual criteria. Thus shape became the predominant 
feature, and even cytological stains, which are actually cytochemical tests 
by nature, were used primarily for their value in making cellular detail 
visible or more distinct under the microscope. It is the task of contempo- 
rary biology to replace such static and purely formal descriptions of living 
systems by dynamic concepts which will define these systems objectively by 
reference to theirinherent properties expressed in measurable standard units, 
rather than in terms of subjective impressions. To quote from an earlier 
paper (Weiss, 1950), ‘‘We must form the habit, for instance, of viewing the 
shape of a given cell, tissue, or organ not as a static feature, but as a phase, 
transitory or terminal, in a continuous chain of transformations; asa cross- 
section through a stream of processes proceeding in time. Every mor- 
phological criterion is but an expression of antecedent processes that have 
contributed to its formation, Distance, in this version, becomes an index 
of past attractions, separations, intercalations, ete.; shape, other than 
spherical, the result of distortions due to elongation, secretions, shifts, 
contractions, expansions, differential growth, and the like; 

In this view, morphological appearance becomes but an incidental 


sign of prior morphogenetic transformations, and cell shape can be under- 
stood only in terms of the chains of interactions through which it has come 
about. Form thus loses its connotation of absoluteness and assumes the 
role of a mere sample taken during the course or at the end of a continuous 
formative process. ‘To translate organic form into such objective and 
‘morphostatics”’ into’ morphodynamics’— will take 


dynamic terms 
quite some time and effort. In the meantime, we shall have to continue, 
for convenience, to use such picturesque designations as fusiform, stellate, 
arborized, crenated, flame cells, goblet cells, ete. But our ultimate goal 
must remain to express them (a) m quantitative terms, and () in terms 
of their genesis. The classical case of a program of this sort has been stated 
by d'Arey Thompson (1942) in his book on Growth and Form. The present 
paper offers a small contribution to this program. It will mainly sketch 
the theoretical premises of our work and some illustrative experimental 
results, while the detailed report of the experiments will follow in a later 


paper by Garber. 
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Cell shapes can be classified according to causes in two distinct groups 
depending on whether they are due to (1) passive deformation from the 
outside, or (2) active transformation by cellular activities. Examples of 
(1) are the configuration of cells closely packed into a confined space 
(e.g., blastomeres; epithelial cells, F. T. Lewis, 1933; plant meristem, 
Matzke, 1947), following general geometrical rules and biologically of rela- 
tively minor interest. Class (2) can be subdivided into two types, namely, 
(a) deformations within a cell caused by its own internal products (e.g., 
local accumulation of secretion granules; fibrous protrusions, such as 
flagellae, sperm tails, etc.; crystalline endoskeletons, W. J. Schmidt, 1924); 
and (b) transformations of the whole cell body into elongate, discoid, 











e 
FIGURE 1 


Series of cell forms observed in fibroblast cultures in vitro 
semidiagrammiatic ) 


ramified, etc., shapes (e.g., mesenchyme cells, neurons, chromatophores, 
glia cells). It is with this last group that we shall deal in the following. 

Figure | shows a graded series of cell samples commonly encountered in 
mesenchyme and connective tissues, both in the body and in tissue cultures. 
Observation has shown that these are widely interconvertible, that is, the 
same cell can assume one or another of these guises depending on circum- 
stances. We shall subject this series to morphogenetic analysis in the 
outlined sense. This involves three steps. First, we must express the 
diverse shapes in appropriate parameters that will permit us to range them 
along a single quantitative scale. Next, we must try to correlate the 


observed morphological variations with some systematic variations either 
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within the cells or in the environment. Lastly, having succeeded in this, 
we must determine the physical basis, the “‘mechanism,”’ underlying this 
correlation; for it is one thing to ascertain, let us say, that a cell moves the 
faster, the thinner the medium (e.g., Wiilmer,1933), and another thing 
to explain just why. 

The direction in which to look for such an explanation was clearly pointed 
by earlier observations (Weiss,1929), in which the two extreme forms of 
the above series (Fig. 1) were found to be characteristically associated with 
fluid-rich and highly condensed plasma media, respectively. Subsequent 
work (Weiss, 1934, 1944a, 1945; Weiss and Wang, 1944) confirmed the 
dependence of cell shape and cell movement on the physical structure 
of the medium. This is schematically summarized in the following dia- 
gram (reproduced from Weiss, 1949), which represents crudely the fibrin 
net of a clot, the lower part of which has been stretched, and which con- 


FIGURE 2 


Effect of stretch (arrows) on the structure of the fibrin reticulum of a plasma clot 
and on the shape of cells contained therein. (From Weiss, 1949). 


tains mesenchyme cells in its different regions (Fig. 2). Where tension 
has aligned and condensed the fibrin (lower part), cells are bipolar and 
elongate in the direction of the fibrin strands, while in the relaxed portions 
(top), cells are multipolar and plump, with intergrades in the intermediate 
zones. These cell forms evidently correspond closely to the series pictured 
in figure 1. Yet, since the whole variety can also be observed, often side 
by side, in unstretched clots, one had to search for a more comprehensive 
explanation. The clue in this search was that the variety of cell forms 
was Clearly correlated with some variable property of the medium, evi- 
dently one that varied systematically with stretch, but at random in 
ordinary clots. 

Cell Parameters.—-For a quantitative formulation of cell shape, we choose 
its deviations from the simplest form, the spheroid (which most cells 
resume prior to mitosis and after death). Ignoring the third dimension, 
all forms of figure 1 can be derived from a sphere by extending certain 
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radii, varying in number and distribution: in a, there are 7 or 8 such major 
projections in relatively uniform distribution along the surface (multipolar 
cell); in 6, there are 3 major ones, about equidistant (tripolar cell); in ¢, 
3 major ones, of which two are almost antipodal to the third (pseudobi- 
polar cell); and in d, e, f, 2 major ones (bipolar cell). The number of 
major processes thus furnishes us with one quantitative scale. As the 
number of processes decreases, bipolarity becomes more marked, and it 
becomes possible then to distinguish between length and width of the cell. 
We shall use the ratio of length over width as index of cell elongation. 
The shape of the similarly elongated nucleus can likewise be expressed 
by the ratio of longest over shortest nuclear diameter. These three vari- 
ables will be designated in the following as NV for number of major proc- 
esses; F, for elongation of the cell body and £, for elongation of the 
nucleus. 

Factors of Deformation.—_Deviation from sperical shape involves shifts 
of the cell content, with some parts extending, others receding. For 
instance, when a cell attaches itself to a flat surface, its border tends to 
expand in practically all directions at the same time, thus flattening the 
mass to a disk at the expense of height (cf. Weiss, 1944a). Forms such as 
those in figure 1, on the other hand, result from the extension of the cell 
surface at relatively fewer points, where tongues of cytoplasm advance 
centrifugally, while the intermediate areas sag, assuming catenary contours. 
The active centrifugal movement in the elongating radii produces radial 
stresses, while the sagging parts of the surface come under tangential 
stress. These stress patterns can sometimes be made visible by impreg- 
nating the correspondingly aligned fibrous units in the cytoplasm with 
silver (e.g., Fig. 12 in Weiss, 1944a) or by observation under the electron 
microscope (Porter, Claude and Fullam, 1945) and thev are valuable clues 
to the mechanism of deformation; they are schematically indicated in 
figure 3. 

The protrusion of any part of the surface is opposed by the leveling 
forces of surface tension and elasticity. Therefore, in order for any protru- 
sion to last more than just momentarily, the force and rate of its thrust 
must be high encugh to overcome these opposing forces. Evidently, the 
average number of projections to be expected from any given cell will de- 
pend on the incidence of such break-throughs. 

In surveying different cell types, it is noted that those listed above as 
Class 2a can often project long arms into structureless (liquid) media by 
virtue of endocellular fibrous or crystalline axial filaments. Unlike these, 
cells of Class 26, with which we are here concerned, are devoid of solid endo- 
skeletons and seem to be unable to form long processes without the aid of 
an outside support—generally speaking, of an interface—to which they 
can cling (cf. Weiss, 1945). Surface protoplasm is drawn out along such 
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interfaces by adhesive forces which vary with the gross physical properties 
and the subtler molecular relations along the contact surfaces and which 
may be lumped here for brevity under the common term of “interfacial ten- 
sions.”’ Along a flat interface, marginal protoplasm will thus be drawn out as 
a thin veil with a broad front (W. H. Lewis, 1950; Weiss, 1944a), while along 
a cylindrical interface, e.g., the surface of a fiber, it will be drawn out into 
a pointed filament. Without going into detail, it should be mentioned 
that the forces needed for the formation of such a cell process seem to be 
provided partly by the interfacial pull and partly by intracellular meta- 
bolic energy, neither of them sufficient to do the work alone (see below). 
Perhaps the interface, by orienting the movement of the protoplasmic 
surface, just aligns crucial molecular chains, and metabolic energy then 
links them by polymerization into a framework consistent enough to 
withstand retractive forces. 

Be this as it may, the morphogenetic significance of the interface as such 
is well substantiated. Hence, any fiber of the extraneous medium that 
makes contact with, or pierces, the cell surface, will call forth a little 
pointed projection, vaguely comparable to the meniscus rising on a wire 
partly immersed in a wetting liquid. 

Now, once established, such a protrusion acts as a viscous flow, to which 
hydrodynamic considerations may be applied. That is, it exerts a drain- 
ing suction on its surroundings in proportion to its momentum. In con- 
sequence, two or more projections extending concurrently from the same 
surface cannot proceed independently, but are in competitive interaction 
with one another, as well as with the cohesive and elastic forces of the 
system. They will contend for the influx of protoplasm among themselves 
according to their relative momentums, with the larger and stronger ones 
tending to drain (hence resorb) the smaller and weaker ones. The appli- 
cation of this principle of tug-of-war among cell processes to the mechanics 
of nerve fiber growth and branching has been presented previously (e.g., 
Weiss, 1941, 1944b). A general consideration of systems of this order 
leads to the following rules: 

(1) An advancing process tends to prevent the emergence of another 
process within a certain radius around its base by competitive drainage. 
(2) This inhibitory effect declines with increasing distance from the base 
of the advancing column. (3) The effect increases with increasing mass 
and velocity of the advance. (4) Other processes already in existence 
within the drainage area can proceed if their mass and velocity are large 
enough to entertain a counter-drain. Otherwise they will be resorbed. 

To these four dynamic rules must be added a fifth one based on a micro- 
morphological consideration. If we assume that the formation of a cell 
process along an extraneous fiber involves the aligning of fibrous units 
in the cell in the direction of the advance, the probability of the emergence 
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of a new process at a given contact point along the cell surface will vary 
inversely with the difficulty with which molecular realignment at that 
point can be brought about. The difficulty will be least where orientation 
is already predominantly perpendicular to the border (e.g., at the advancing 
front of a stream) while it will reach a maximum where the orientation is par- 
allel to the border, particularly if the elements have become linked into larger 
surface-parallel chains. Rule (5) therefore states that a contact stimulus 
of given magnitude will be more likely to produce a local ‘“microlesion”’ 


¢ 


FIGURE 3 


Reaction of the cell surface (stippled) to contact with extraneous fibers of 


different sizes 


or breach of the cell surface if it is applied near the tip, than if it is near 
the flank, of a moving cell part. As reported in an earlier paper (Weiss, 
1944a), observations on phagocytosis point to the same conclusion, as 
granules are ingested only along those parts of the cell border that seem 
to possess a radial microstructure, but not across zones with tangential 
orientation (for supporting evidence, see also Hancox, 1949). 

Relation between cell processes and texture of medium.—The five cited 
rules permit us to predict the conditions making for greater or smaller 
average numbers of processes in given cell types in a medium containing net- 
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works of fibers such as the fibrin reticulum of a blood plasma clot. Figure 3 
represents schematically small samples of three cell surfaces (cytoplasm 
stippled), each in contact with three equidistant fibrin threads (black) 
of identical (a) or different (6, c) size. In line with what was said above, 
we conclude that at each contact point a protoplasmic process will be 
sucked forth by interfacial tension, the draining force being proportional 
to the size of the fiber. Larger fibers thus have a competitive advantage 
over smaller ones. In figure 3a, all three fibers being of equal size, 
equilibrium is attained with the establishment of three processes. In 
figure 3b, the larger middle fiber weakens the processes along the two 


FIGURE 4 
Morphology of mesenchyme cells in three media of (from left to right) increas- 
ing plasma concentrations. 


outer ones. In figure 3c, the largest outflow along the bottom fiber 
will drain both the middle and the top processes. We may conclude then 
that for a given density of the fibrous fabric, the number of cell processes 
that can be sustained per unit surface will be the larger, the more uniform 
the fiber diameters are, and the smaller, the greater the diversity of fiber 
sizes. Furthermore, according to rule (5), the stronger the oriented flow 
(the larger the fiber), the smoother will be the edges. In consequence, the 
number of processes should vary inversely with the proportion of large 
fibers in the medium. 

To test this assumption, we studied cell shapes in plasma clots prepared 
in a graded series of fibrin textures. The reactions to be expected are 
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illustrated in figure 4, showing in the upper row, a small segment of a cell 
border exposed to three clots of mounting fibrin concentrations, and in the 
lower row, the shapes displayed by the respective cells (for simplicity, 
cytoplasm and fibrin are not drawn to the same scale). As fibrin concentra- 
tion increases from a to ¢, more and more elementary filaments combine 
into larger fibers, resulting in an increase of both the range of the fiber 
spectrum and the proportion of large fibers. These three diagrams, 
corresponding essentially to the three phases of figure 3, summarize the 
conclusion drawn above, that the higher the incidence of large fibers in 
the medium, the greater is the probability that the cell surface will be vigor- 
ously drawn out along one of them, thus automatically suppressing the 
formation of other sprouts for some distance. The diagrams thus lead 
us to expect an inverse statistical relation between coarseness of medium 
and average number of major cell processes. 

Stated succinctly, the principle involved is as follows. Small fibers 
evoke weak protrusions with a small draining radius, hence, able to coexist 
without mutual interference even at close range. Larger fibers, evoking 
a larger outflow, establish a wider drainage field. Multiple extensions 
within such a common field, being in hydrodynamic and elastic competition, 
could only persist if all of them were of exactly identical strength, which 
would presuppose uniformity in the size of the guide fibers. The greater 
the range of guide fiber sizes, the less will be the chance for such an equili- 
bration to occur, hence, the greater will be the probability that a single 
process will assume dominance. It is for this reason that we must expect 
the number of cell processes to be a function not only of maximum fiber 
size but also of the range of fiber sizes in the medium. 

Experimental Tests. Media of controlled texture were obtained by two 
different methods. In one set of experiments, culture media were made up 
of mixtures of fowl plasma and chick embryo extract in varying ratios, 
ranging from 10 to 90° plasma concentration. With increasing con- 
centrations, the tendency of the fibrin filaments to aggregate to larger 
bundles rises steadily, thus providing us with textures of increasing coarse- 
ness. Although this series gives us the desired physical gradations, it 
also introduces chemical differentials due to the varying ratios of plasma 
and embryonic substances. Therefore, in a second set of experiments 
plasma and extract were mixed in equal amounts (40%, plasma concentra- 
tion), but coagulated at different pH over a range of from 5.6 to 8.0. 
It is known and has been crucially confirmed by the electron micrographs 
of Hawn and Porter (1947) that the lower the pH, the coarser become 
the bundles of fibrin. This then provides us with a means of grading 
physical texture while keeping chemical composition constant. Since 
both series gave identical results regarding cell shape, physical texture 
is proved to be the relevant variable. 
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The cultures were standard heart fibrocyte cultures having been passed 
through five transfers to obtain cell strains of relatively uniform character. 
Clots coagulated at unphysiological pH values were quickly buffered back 
to normal by the growing tissue. Because surface and interior of a clot 
are vastly different in texture, only cells in the depth of the clot were re- 
corded, excluding any that had spread along the coverslip or the air inter- 
face. Cells were measured in camera lucida drawings of either the living 
cultures or the fixed preparations. Samples were taken at random. The 
ranking of a given process as major (1.e., affecting cell shape significantly ) 
or minor, is, of course, somewhat arbitrary, but since large numbers of 
cells were used and the same standards were applied in the rating through- 

out, the results of the different 

pH SCALE series were comparable and con- 

a ; r —£ : sistent. In bipolar cells, length was 

measured from tip to tip; width, 

across the cell body at the level of 

the nucleus. The following results 

are based on measurements of a total 
of 4759 cells. 

General Cell Shape.-Figure 5 
expresses the most general relation 
between cell morphology and struc- 
ture of the medium, namely, an 
increasing proportion of elongate 
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(spindle) cells (types d, e, f, of 
Fig. 1) in the cell population with 
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FIGURE 5 ages of spindle cells for the five 

Proportion of spindle cells in plasma plasma concentrations tested (lower 
clots of different textures, abscissa), fitting a straight line 
which extrapolates to 100°; spindle 

cells for a theoretical plasma concentration of 100%. The proportion 
of stellate cells (types a, 6, c, of Fig. 1) decreases correspondingly. In 
media that had been coagulated at different pH (upper abscissa), the 
proportion of spindle cells likewise rises steadily with increasing coarseness 
(decreasing pH), and the resulting values (full circles) also roughly fit a 
linear function. The two curves representing the concentration and the 
pH functions have been brought to coincidence in the graph by an appro- 
priate choice of abscissa units. The fact that as good a coincidence as is 
evident could be obtained proves that the relevant variable in determining 
cell morphology is really the texture of the medium without regard as to 
whether this is accompanied by variations in chemical concentration (lower 
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scale) or not (upper scale). The validity of this conclusion is further 
confirmed by the fact reported below that if other cell parameters are 
plotted in the same coordinate systems as those used in figure 5, they also 
fit a common curve. According to these data, therefore, as far as effects 
on cell morphology are concerned, 





the texture of a 10°% plasma clot = '°° 


at normal pH is equivalent to that 
of a 40% plasma clot at pH 8, and 


an 80° @lot at normal pH corre 


sponds to a 40°, clot at pH 5. 
Number of Cell Processes.— Figure 
6 presents a more detailed break 
down according to shapes of cell 
populations (2417 cells) in clots of 
graded plasma concentrations. On 





each ordinate are added the per- 
centages of cells with 2, 3, 4 and 5 
major processes. It can readily be 
seen that as plasma concentration 
rises, the average number of cell 
processes (.V) declines steadily, 
again essentially in linear order and 
with a remarkable degree of regu- 
larity. The data thus bear out 
fully the prediction derived above 





from theoretical considerations. 
Elongation. Parallel with the re- 
duction of major cell processes goes 





. ° ° $0 

a progressive elongation of the cell, euahon “eenenvuarien 
leading from multipolarity through 
pseudobipolarity to strict bipolarity 


and then continuing, within the bi- 


FIGURE 6 
Proportions of cells with 2 (black area; 
white circles), 3 (dense stippling; crosses), 
$ (light stippling; dark circles), and 5 
ness (Fig. 1). The quantitative (white) major processes in clots of different 
expression of this transition is a textures 


polar range, to increasing slender- 


steady increase in the index of 

elongation (/2,), i.e., the ratio of length over width, of the cell body. The 
nuclear proportions (/¢,) conform to those of the cell body. Both series of 
values are plotted in figure 7 in relation to clot texture; open circles referring 
to plasma concentration (lower abscissa), full circles to pH at coagulation 
(upper abscissa). Clearly, the values of /, and those of /, define each a 
single continuous function of non-linear character. It will be noted that 
the numerical relation between the concentration and pH scales is the 
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same as in figure 5. The fact that with this setting, the values of both 
experimental series fit the same curve, justifies the assumption that lower 
and upper abscissae are truly interconvertible in point of morphological 
effects of clot texture. 

These are empirical curves. In view of the complex factors involved, 
it would be idle to try, at this early date, to interpret their detailed shape 
in physical terms. Yet, their general character can be derived from our 
theoretical concept. A bipolar spindle cell may be regarded essentially 
as two processes drawn out along a common fiber in opposite directions, 
with one tip leading, the other trailing, and the rest of the cell body 

stretched out between them. Ac- 

pH SCALE cording to our assumption, the 

: ois : pulling force is a function of inter- 
: facial tensions along the leading 

edge, which is the front line of con- 

sak tact between protoplasm and fiber. 

CELL This force increases in some direct 

proportion to the circumference, 
hence the radius, of the fiber. 
Elongation results from the lag of 
the rear process in following the 
advance. This lag is due, among 
other things, to the adhesion of the 
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rear process to the fiber, such adhe- 

sion being a function of the size of 
- L 4. | a a a 

Oo 0 30 50 70 90 the contact area, hence, propor- 





4. 4 A. 





e . . 
7 PLASMA CONCENTRATION tional to the square of the radius of 


FIGURE 7 the fiber. Therefore, since the 

Ratios of length over width for cell bodies advance increases linearly with the 

and nuclei in clots of different textures. radius, while the drag increases as 

the square of the radius, elongation 

resulting from the disparity between rate of advance at the tip and rate 

of detachment at the rear, will be the greater, the larger the diameter 

of the fiber. Consequently, media with higher plasma concentrations, 

or coagulated at lower pH values, hence containing larger proportions 

of heavy fibers, must also be expected to yield cells with higher average 
E, values, which is precisely what the curves of figure 7 demonstrate. 

An unexpected feature revealed by these curves is the remarkable simili- 
tude between cellular and nuclear proportions. The ratio /,:/, has an 
almost constant value for all observed cell shapes, averaging 3.4, with 
a range of from 3.1 to 3.8. This expresses a rather strict conformance 
of nuclear shape to the proportions of the cell body, but the mechanism 
by which this is brought about is not at all clear. It is apparent that 
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the elongation of the nucleus is directly and quantitatively related to the 
polarized structure and movement of the surrounding cytoplasm, but in 
view of our great ignorance regarding the factors determining nuclear 
shape in general, it remains uncertain whether this is a matter of passive 
deformation or of an active nuclear rearrangement. The notion of a sheer 
passive lateral compression of the nucleus might seem questionable in 
view of the common observation that other spherical cell inclusions, such 
as fat globules, remain spherical; however, this difference could be ac- 
counted for by the difference in consistency between nucleus and fat 
globules, the former being gelated (to judge from motion pictures of 
fibroblast cultures), the latter liquid. Even though the meaning of the 
peculiar constancy of the value of /:/, throughout the range of shapes 
observed remains obscure, the fact as such is noteworthy in connection 
with any future attempts to explain nuclear morphology. 

Rate of Cell Movement.—In deriving elongate shape from the disparity 
in the rates of movement of the front and rear poles of the cell, the pre- 
ceding chapter has already brought out the close relation between cell 
morphology and cell movement. As the cell elongates, the elastic strain 
within its cytoplasm (presumably in its fibrous framework) may be assumed 
to increase proportionately. The adhesive resistance of the rear end will 
be overcome when this strain has reached the proper magnitude, that is, 
when a certain critical length has been attained, and this will then be 
the shape in which we shall find the cell as a whole in locomotion. Just 
how the gliding type of locomotion characteristic of these cells (and simi- 
larly of nerve fiber tips) is achieved, is still rather obscure. It has been 
compared to the locomotor mechanism of amoebae (W. H. Lewis, 1950), 
but even the explanations offered tor the latter are not yet fully satis- 
factory. For our obiects, both the motor mechanism and the source 
of the motile energy remain to be revealed. At present, the most appro- 
priate assumption seems to be a dual source, lying partly in the interfacial 
tensions, partly in intracellular activities supported by metabolic energy. 
For the nerve fiber, such a dual ‘‘pull-push’’ mechanism of elongation has 
been suggested previously by one of us (see Weiss, 1944b), and pending 
proof to the contrary, the fibroblast may be regarded in the same light. 
That is, we may assume that interfacial tension would provide the “‘lead 
force,’ and that metabolic energy, through some unknown intracellular 
mechanism would sustain the translatory movement. 

On this assumption, the rate of gliding should be mainly a function of 
the metabolic state of the cell, rather than of the size of the fibers serving 
as substratum. Elongation and locomotion thus do not necessarily 
follow the same rules. 

Rate of migration was determined by tracing cells at the margin of a 
culture (always well inside the clot) after 24 and 48 hours of cultivation 
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and determining their average hourly progress. Figure 8 gives the average 
rates obtained for 1582 cells so measured in graded plasma concentrations, 
as well as the values for the 10 fastest cells in each series. Both curves 
are identical in character. Their extrapolated origin at the null point is 
consistent with the fact that fibroblasts inside a liquid drop (0% plasma 
concentration) cannot move at all. The rate of progress then increases 
very rapidly with increasing plasma concentration, reaching a maximum 
of about | mm. per day near a plasma concentration of 40%, and retaining 
this value rather constantly throughout the higher concentration range. 

In interpreting these curves, it is first to be noted that, as we indicated 
before, degree of elongation and rate of migration follow different rules, as 
the former continues to rise (Fig. 7), even after the latter has leveled off. 
Moreover, since the increase of fibrin fiber size (pertinently reflected in 
progressive elongation) continues steadily through the whole concentration 
range, the fact that migratory rate in the upper range remains constant 
confirms our earlier supposition (see above) that locomotor rate would 
be independent of fiber size. If this is the case, then the drop of the curves 
of migratory rate in the lower concentration range must be ascribed to 
quite different causes than that of the elongation curves, a fact which is 
further accentuated by the opposite sense of the respective curvatures. 

The slower migration in more dilute plasma can evidently be explained 
on the following grounds. The main clue is the realization that the data 
of figure 8 were obtained not by a continuous recording of the velocity of 
cell movement, but simply by measuring the total displacement of cells at 
the end of a 24-hour period. This would give us a true measure of actual 
velocity of motion only if cells were moving continuously and at a constant 
rate throughout the measured interval. However, if the movement is 
discontinuous, that is, if it is interrupted by stationary phases of un- 
known duration, slower progress would merely indicate more or longer stops 
without necessarily implying a slower basic rate. The incidence of such 
arrests is much higher in the lower plasma concentrations because of the 
greater frequency of divergent cell processes (Fig. la-c), as explained 
earlier in this paper. It is obvious that no cell can advance in two differ- 
ent directions at the same time. Therefore, as long as such processes 
stalemate each other, the cell as a whole is stalled. Sooner or later one 
or the other of the competing projections wins out in the resulting tug-of- 
war, the losing one is withdrawn and the cell can then proceed until the 
next bifurcation arises, producing another delay, and so on, Our curves, 
therefore, express cumulative delays rather than retardations of proto- 
plasmic flow. 

Another factor vitiating our ‘‘rate’’ curves is the irregularity of the course 
taken by a cell in any but a strictly oriented medium. Distances in our 
cultures (Fig. 8) were measured as straight-line connections between the 





Vou. 38, 1952 ZOOLOGY: WEISS AND GARBER 


end points of a 24-hour trek. Actually, however, all multipolar cells had 
zigzagged about this average course, hence, run a greater distance than that 
measured, in proportion to the random divergence of their lead processes. 
Therefore, on this score, too, our “rate’’ values definitely underrate loco- 
motor speed, the more so the greater the average number of divergent 
cell processes, that is, the lower the plasma concentration. 

Bipolar cells, by contrast, present a wholly different situation, if we 
remember that bipolarity itself is an expression of gliding on relatively 
heavy and unintersected path fibers. Having a straight runway, such cells 
are not subject to the periodical delays of highly intersected networks, 
for even if they do project more than one advancing process, the fact that 
all such multiple extensions lie in 
the common line of advance (Fig. 1- 
d), precludes competitive delays. 
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conclusion that the curves of figure 
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the same in both groups. The l-mm. displacement between 24 and 48 
hours of cultivation is in good accord with the maximum rate recorded in 
a preceding investigation (Weiss, 1945) for the second day of strictly 
oriented migration of Schwann cells along glass fibers embedded in plasma 
clots. But our data do not corroborate Willmer’s (1933) contention that 
migration is slower in denser plasma clots. Since neither the experimental 
conditions nor the methods of measurement nor the sizes of the samples 
used in his and our experiments were comparable, it is impossible to deter- 
mine the source of this discrepancy. On the other hand, results quite 
comparable to ours were obtained in amoebae (Pitts, 1933), in which the 
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rate of locomotor progress was found to be much slower during multipolar 
(‘‘multipodal’’) phases than when they were in the bipolar (‘“‘monopodal’’) 
form. 

Conclustons.—We have demonstrated, in the foregoing, the feasibility 
of resolving the morphology of a single cell type—the heart fibroblast— 
quantitatively into terms of its morphogenesis. Pictorial description 
has been replaced by a number of curves, representing the whole variety 
of observed cell forms as function of a single variable—the physical struc- 
ture of the medium. The trend of these curves, at the same time, has 
verified the basic assumption that the shape of cells of this kind can be 
derived from an understanding of the manner in which they react to 
fibrous substrata. It is significant to note that this deviation rests 
on statistical considerations, reemphasizing the thesis that morphogenesis, 
by its very nature, is largely statistical. Neither are the individual cells 
even of a pure cell strain all exactly identical, nor is any organic medium 
ideally uniform and homogeneous. All that is predictable, therefore, is 
the average constitution of the cells and the average structure of the 
medium, hence, the probability of observing a given morphogenetic re- 
sponse within a given range. To this way of thinking, rigid and micro- 
precise determinism of cell shape is mere illusion. We offer this demon- 
stration of how static morphological symbols can be made to yield to 
dynamic-statistical interpretations, as encouragement to apply a similar 
analysis to other cell types. 

Such comparative studies of different cell types would have another 
point in their favor. They could open a wholly new way of getting in- 
formation about certain properties of living cells uncomplicated by opera- 
tive tampering. Let us illustrate this by an example. Remember that 
our data have characterized one single cell type only—the heart fibroblast— 
in a graded series.of environments. All events recorded—the projection 
of cytoplasmic processes; competitive interaction among multiple proc- 
esses; rate of advance; elongation; surface reorientation; adhesion to, 
and detachment from, substratum; etc.—have been predicated on the 
particular physical and chemical properties of that particular cell type. 
Cell shapes are no more “caused” or ‘‘determined” by the environment 
as such than they are by any autonomous intracellular factors as such. 
Cell and environment are two interdependent variables of which we merely 
kept the one—cell constitution—practically constant throughout the tested 
range. Now, with this set of data on hand, we can reverse the procedure 
and go on to comparing different cell types, keeping the environment 
constant. We can take, for instance, other kinds of mesenchyme cells, 
nerve sheath cells, sarcoma cells, and so on, determine their parameters 
(number of processes, elongation of cell body, nuclear deformation, rate 
of locomotion, etc.) in media of given plasma concentrations and pH 
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and compare these then with the values obtained from our standard fibro- 
blast strain under corresponding conditions. Any consistent differences 
‘an then be taken as indicators of physical or chemical differences between 
the respective cell strains. For instance, an upward displacement of the 
curves of figures 5 and 6 (signifying reduced frequency of cell projections) 
would indicate greater surface stability or swifter protoplasmic movement ; 
an upward shift of the curves in figure 7 (signifying greater elongation), 
lower viscosity or a lower modulus of elasticity or greater adhesiveness; 
an increase in the ordinate ratio between the upper and the lower curve 
of figure 7 (signifying reduced relative elongation of the nucleus), greater 
turgor or rigidity of the nucleus or reduced strength of the deforming 
forces in the cytoplasm; etc.—all of these for the particular test strain with 
reference to our standard strain. Similarly, it would be possible to detect 
within a given cell strain subtle alterations of physical properties attending 
differentiation, aging, somatic mutations or pathologic changes. 

On the purely technical side, our results illustrate plainly the dangers 
inherent in the uncritical use of the spread of tissue cultures as a measure 
of growth. Considering that the increase in area is primarily a result 
of cell migration, and bearing in mind the complex dependency of the rate 
and amount of migration upon physical properties of both the cell body 
and the medium, it is plain that whenever an agent, drug, radiation, etc., 
has been observed to affect the rate of expansion, this cannot be safely 
interpreted as an effect on growth as such unless the agent in question has 
been proved to affect neither the texture of the medium nor such physical 
conditions in the cell as are relevant to its mode and rate of spreading. 
The interdependence between cell configuration and cell growth pointed 
out on a previous occasion (Weiss, 1949) calls for additional caution in the 
assaying of so-called growth stimulants and growth inhibitors by geometric 
measurements on tissue cultures. 

Summary.—The various cell shapes exhibited by a sample cell type 
the heart fibroblast in tissue culture-—have been expressed on a single 
quantitative scale and explained in terms of the reaction of the cell to the 
fibrous texture of the medium. The texture of the fibrin reticulum was 
varied systematically in two graded series, one of plasma concentrations, 
and the other of pH values at coagulation. Measurable criteria of cell 
shape (number of cell processes; frequency of bipolars in the population; 
average ratio of length over width of cell body; nuclear elongation) were 
found to vary accordingly and in the predicted sense. Rate of migration 
could likewise be related to plasma texture. It has thus been possible, 
at least for this obiect, to resolve in principle static cell morphology into 
quantitative terms of the dynamics of cell morphogenesis. 

* Research aided by the Wallace C. and Clara A. Abbott Memorial Fund of the 
University of Chicago and a graot-in-aid from the American Cancer Society upon the 
recommendation of the Committee on Growth of the National Research Council. 
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ERRATA 


In my article, ‘“A New Discussion of the Changes in the Earth’s Rate of 
Rotation,” in the January issue on page 9 in line 13 from the bottom 7 miles 
should read 2.5 miles and in line 10 from the bottom 200 miles should read 
SO miles. 


DirK BROUWER 
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